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The Weak Field Limit of Higher Order 
Gravity 



The Higher Order Theories of Gravity - f{R, Raj^R"^) - theory, where R is the Ricci scalar, R^p 
is the Ricci tensor and / is any analytic function - have recently attracted a lot of interest as alter- 
native candidates to explain the observed cosmic acceleration, the flatness of the rotation curves 
of spiral galaxies and other relevant astrophysical phenomena. It is a crucial point testing these 
alternative theories in the so called weak field and newtonian limit of a /(-R, Ra/sR"^) - theory. 
With this "perturbation technique" it is possible to find spherically symmetric solutions and com- 
pare them with the ones of General Relativity. On both approaches we found a modification of 
General Relativity: the behaviour of gravitational potential presents a modification Yukawa - like 
in the newtonian case and a massive propagation in the weak field case. When the modification of 
the theory is removed (i.e. f{R, RapR""^) = R, Hilbert - Einstein lagrangian) we find the usual 
outcomes of General Relativity. Also the Noether symmetries technique has been investigated to 
find some time independent spherically symmetric solutions. 

Arturo StabileF] 



e - mail address: arturo.stabile@gmail.com 
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Introduction: History and Motivations 



0.1 General considerations on Gravitational theories 

General Relativity (GR) is a comprehensive theory of spacetime, gravity and matter. Its formu- 
lation implies that space and time are not "absolute" entities, as in Classical Mechanics, but dy- 
namical quantities strictly related to the distribution of matter and energy. As a consequence, this 
approach gave rise to a new conception of the Universe itself which, for the first time, was consid- 
ered as a dynamical system. In other words. Cosmology has been enclosed in the realm of Science 
and not only of Philosophy, as before the Einstein work. On the other hand, the possibility of a 
scientific investigation of the Universe has led to the formulation of the Standard Cosmological 
Model [LD which, quite nicely, has matched with observations. 

Despite of these results, the study of possible modifications of Einstein's theory of gravitation 
has a long history which reaches back to the early 1920s [12l[3l|4l|51|6ll3. While the proposed early 
amendments of Einstein's theory were aimed toward the unification of gravity with the other inter- 
actions of physics, like electromagnetism and whether GR is the only fundamental theory capable 
of explaining the gravitational interaction, the recent interest in such modifications comes from 
cosmological observations (for a comprehensive review, see [8J). Such issues come, essentially, 
from Cosmology and Quantum Field Theory. In the first case, the presence of the Big Bang singu- 
larity, the flatness and horizon problems flU led to the statement that Cosmological Standard Model, 
based on GR and Standard Model of Particle Physics, is inadequate to describe the Universe at ex- 
treme regimes. These observations usually lead to the introduction of additional ad-hoc concepts 
like dark energy/matter if interpreted within Einstein's theory. On the other hand, the emergence 
of such stopgaps could be interpreted as a first signal of a breakdown of GR at astrophysical and 
cosmological scales [ITOl ITTl . and led to the proposal of several alternative modifications of the 
underlying gravity theory (see [fT2l for a review). Besides from Quantun Field Theory point view, 
GR is a classical theory which does not work as a fundamental theory, when one wants to achieve 
a full quantum description of spacetime (and then of gravity). 
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While it is very natural to extend Einstein's gravity to theories with additional geometric de- 
grees of freedom, (see for example lfl3l [T4l [TSl for general surveys on this subject as well as IIT6II 
for a list of works in a cosmological context), recent attempts focused on the old idea of modifying 
the gravitational Lagrangian in a purely metric framework, leading to higher order field equations. 
Such an approach is the so-called Extended Theories of Gravity (ETG) which have become a sort 
of paradigm in the study of gravitational interaction. They are based on corrections and enlarge- 
ments of the Einstein theory. The paradigm consists, essentially, in adding higher order curvature 
invariants and minimally or non-minimally coupled scalar fields into dynamics which come out 
from the effective action of quantum gravity ifTTl . 

The idea to extend Einstein's theory of gravitation is fruitful and economic also with respect 
to several attempts which try to solve problems by adding new and, most of times, unjustified in- 
gredients in order to give a self-consistent picture of dynamics. The today observed accelerated 
expansion of the Hubble flow and the missing matter of astrophysical large scale structures, are 
primarily enclosed in these considerations. Both the issues could be solved changing the gravi- 
tational sector, i.e. the l.h.s. of field equations. The philosophy is alternative to add new cosmic 
fluids (new components in the r.h.s. of field equations) which should give rise to clustered struc- 
tures (dark matter) or to accelerated dynamics (dark energy) thanks to exotic equations of state. 
In particular, relaxing the hypothesis that gravitational Lagrangian has to be a linear function of 
the Ricci curvature scalar R, like in the Hilbert-Einstein formulation, one can take into account an 
effective action where the gravitational Lagrangian includes other scalar invariants. 

In summary, the general features of ETGs are that the Einstein field equations result to be 
modified in two senses: i) geometry can be non-minimally coupled to some scalar field, and / 
or a) higher than second order derivative terms in the metric come out. In the former case, we 
generically deal with scalar-tensor theories of gravity; in the latter, we deal with higher order 
theories. However combinations of non-minimally coupled and higher-order terms can emerge 
as contributions in effective Lagrangians. In this case, we deal with higher-order-scalar- tensor 
theories of gravity. 

Due to the increased complexity of the field equations in this framework, the main amount 
of works dealt with some formally equivalent theories, in which a reduction of the order of the 
field equations was achieved by considering the metric and the connection as independent fields 
[fTSl [T91 1201 12TI I22I . In addition, many authors exploited the formal relationship to scalar-tensor 
theories to make some statements about the weak field regime, which was already worked out for 
scalar- tensor theories more than ten years ago [|23ll . 

Other motivations to modify GR come from the issue of a full recovering of the Mach principle 
which leads to assume a varying gravitational coupling. The principle states that the local inertial 
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frame is determined by some average of the motion of distant astronomical objects [|24ll . This 
fact implies that the gravitational coupling can be scale-dependent and related to some scalar field. 
As a consequence, the concept of "inertia" and the Equivalence Principle have to be revised. For 
example, the Brans-Dicke theory ll25l is a serious attempt to define an alternative theory to the 
Einstein gravity: it takes into account a variable Newton gravitational coupling, whose dynamics 
is governed by a scalar field non-minimally coupled to the geometry. In such a way, Mach's 
principle is better implemented Il25ll26ll27l . 

As already mentioned, corrections to the gravitational Lagrangian, leading to higher order field 
equations, were already studied by several authors |l3l[6l|7]| soon after the GR was proposed. De- 
velopments in the 1960s and 1970s (H |29l [30l [SB El, partially motivated by the quantization 
schemes proposed at that time, made clear that theories containing only a R"^ term in the La- 
grangian were not viable with respect to their weak field behavior. Buchdahl, in 1962 |28|| rejected 
pure theories because of the non-existence of asymptotically flat solutions. 

Another concern which comes with generic higher order gravity (HOG) theories is linked to 
the initial value problem. It is unclear if the prolongation of standard methods can be used in 
order to tackle this problem in every theory. Hence it is doubtful that the Cauchy problem could be 
properly addressed in the near future, for example within 1/R theories, if one takes into account the 
results already obtained in fourth order theories stemming from a quadratic Lagrangian [|33l[34l . 

Starting from the Hilbert-Einstein Lagrangian 



and combinations of them, represent a first obvious choices for an extended gravity theory with 
improved dynamics with respect to GR. Since the variational derivative of £3 can be linearly ex- 
pressed [EllSSl via the variational derivatives of Ci and L2, one can omit £3 in the final Lagrangian 
of a HOG without loss of generality. 




(1) 



the following terms 




(2) 
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Besides, every unification scheme as Superstrings, Supergravity or Grand Unified Theories, 
takes into account effective actions where non-minimal couplings to the geometry or higher order 
terms in the curvature invariants are present. Such contributions are due to one-loop or higher 
loop corrections in the high curvature regimes near the full (not yet available) quantum gravity 
regime [fTTl . Specifically, this scheme was adopted in order to deal with the quantization on curved 
spacetimes and the result was that the interactions among quantum scalar fields and background 
geometry or the gravitational self-interactions yield corrective terms in the Hilbert-Einstein La- 
grangian [36]. Moreover, it has been realized that such corrective terms are inescapable in order 
to obtain the effective action of quantum gravity at scales closed to the Planck one ll37l . All these 
approaches are not the "/m/Z quantum gravity" but are needed as working schemes toward it. 

In summary, higher order terms in curvature invariants (such as R^, RapR^^ , Rap-^sR"^'^^ , 
ROR, or Rn^R) or non-minimally coupled terms between scalar fields and geometry (such as 
4>^R) have to be added to the effective Lagrangian of gravitational field when quantum corrections 
are considered. For instance, one can notice that such terms occur in the effective Lagrangian of 
strings or in Kaluza- Klein theories, when the mechanism of dimensional reduction is used ll38l . 

On the other hand, from a conceptual viewpoint, there are no a priori reason to restrict the 
gravitational Lagrangian to a linear function of the Ricci scalar R, minimally coupled with matter 
lfT9l . More precisely, higher order terms appear always as contributions of order two in the field 
equations. For example, a term like R^ gives fourth order equations ||39ll , R OR gives sixth order 
equations B0ll4T| . R UPR gives eighth order equations [|42l| and so on. By a conformal transforma- 
tion, any 2nd order derivative term corresponds to a scalar fiel(i^: for example, fourth order gravity 
gives Einstein plus one scalar field, sixth order gravity gives Einstein plus two scalar fields and so 
on [|40l l43l . Furthermore, the idea that there are no "exact" laws of physics could be taken into 
serious account: in such a case, the effective Lagrangians of physical interactions are "stochastic" 
functions. This feature means that the local gauge invariances {i.e. conservation laws) are well 
approximated only in the low energy limit and the fundamental physical constants can vary ll44l . 



0.2 Issues from Cosmology 

Beside fundamental physics motivations, these theories have acquired a huge interest in Cosmol- 
ogy due to the fact that they "naturally" exhibit inflationary behaviors able to overcome the short- 
comings of Cosmological Standard Model (based on GR). The related cosmological models seem 

^The dynamics of such scalar fields is usually given by the corresponding Klein-Gordon equation, which is second 
order. 
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realistic and capable of matching with the Cosmic Microwave Background Radiation (CMBR) ob- 
servations |[34ll45ll46l . Furthermore, it is possible to show that, via conformal transformations, the 
higher order and non-minimally coupled terms always correspond to the Einstein gravity plus one 
or more than one minimally coupled scalar fields If33ll40l[47ll48ll49ll . 

Furthermore, it is possible to show that the /(i?) -gravity (/-gravity) is equivalent not only to 
a scalar-tensor one but also to the Einstein theory plus an ideal fluid [|50l . This feature results 
very interesting if we want to obtain multiple inflationary events since an early stage could select 
"very" large-scale structures (clusters of galaxies today), while a late stage could select "small" 
large-scale structures (galaxies today) llTO . The philosophy is that each inflationary era is related 
to the dynamics of a scalar field. Finally, these extended schemes could naturally solve the problem 
of "graceful exit" bypassing the shortcomings of former inflationary models [|46ll5TI . 

In recent years, the efforts to give a physical explanation to the today observed cosmic acceler- 
ation [221 [531 ES [55l have attracted a good amount of interest in /-gravity, considered as a viable 
mechanism to explain the cosmic acceleration by extending the geometric sector of field equations 
[ISa [571 [51 [591161 [Ml 111 [SI IM [Ml [Ml [63 There are several physical 

and mathematical motivations to enlarge GR by these theories. For comprehensive reviews, see 

Specifically, cosmological models coming from /-gravity were firstly introduced by Staro- 
binsky [|451 in the early 80'ies to build up a feasible inflationary model where geometric degrees 
of freedom had the role of the scalar field ruling the inflation and the structure formation. 

In addition to the revision of Standard Cosmology at early epochs (leading to the Inflation), 
a new approach is necessary also at late epochs. ETGs could play a fundamental role also in 
this context. In fact, the increasing bulk of data that have been accumulated in the last few years 
have paved the way to the emergence of a new cosmological model usually referred to as the 
Cosmological Concordance Model (A Cold Dark Matter. ACDM). 

The Hubble diagram of Type la Supernovae (hereafter SNela), measured by both the Supernova 
Cosmology Project [|54l UM and the High-z Team [|52l [79l up to redshift z ^ 1, the luminosity 
distance of la Type Supernovae [[521 [531 [541 [551, the Large Scale Structure [[721 and the anisotropy 
of CMBR [[731 UM are the evidence that the Universe is undergoing a phase of accelerated expan- 
sion. On the other hand, balloon born experiments, such as BOOMERanG [[80l and MAXIMA 
[|8TI . determined the location of the first and second peak in the anisotropy spectrum of the CMBR 
strongly pointing out that the geometry of the Universe is spatially flat. If combined with con- 
straints coming from galaxy clusters on the matter density parameter ^Im, these data indicate that 
the Universe is dominated by a non-clustered fluid with negative pressure, generically dubbed dark 
energy, which is able to drive the accelerated expansion. This picture has been further strength- 
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ened by the recent precise measurements of the CMBR spectrum, due to the WMAP experiment 
11731 mi HH, and by the extension of the SNela Hubble diagram to redshifts higher than 1 [l52ll . 
From this amount of data, the widely accepted ACDM is a spatially flat Universe, dominated by 
cold dark matter (CDM (~ 0.25 ^ 0.3%) which should explain the clustered structures) and dark 
energy (A ~ 0.65 ^ 0.7%), in the form of an "effective" cosmological constant, giving rise to the 
accelerated behavior. 

After these observational evidences, an overwhelming flood of papers has appeared: they 
present a great variety of models trying to explain this phenomenon. In any case, the simplest 
explanation is claiming for the well known cosmological constant A [831. Although it is the best 
fit to most of the available astrophysical data 11731 17411821 . the ACDM model fails in explaining why 
the inferred value of A is so tiny (120 orders of magnitude lower than the value of quantum gravity 
vacuum state!) if compared with the typical vacuum energy values predicted by particle physics 
and why its energy density is today comparable to the matter density (the so called coincidence 
problem). 

Although the cosmological constant [|84l [85l [86l remains the most relevant candidate to in- 
terpret the accelerated behavior, several proposals have been suggested in the last few years: 
quintessence models, where the cosmic acceleration is generated by means of a scalar field, in a 
way similar to the early time inflation [|45l . acting at large scales and recent epochs [l87l[88l : mod- 
els based on exotic fluids like the Chaply gin-gas |l89l|90l|9Tl, or non-perfect fluids [|92l ): phantom 
fields, based on scalar fields with anomalous signature in the kinetic term [|93ll94ll95ll96l , higher 
dimensional scenarios (braneworld) [|97ll98ll99lll00ll . These results can be achieved in metric and 
Palatini approaches [|20l|57l|6llMl[67llMl[Ml[69l|70lIl0^ In addition, reversing the problem, one 
can reconstruct the form of the gravity Lagrangian by observational data of cosmological relevance 
through a "back scattering" procedure. All these facts suggest that the function / should be more 
general than the linear Hilbert-Einstein one implying that HOG could be a suitable approach to 
solve GR shortcomings without introducing mysterious ingredients as dark energy and dark matter 
(see e.g. IMlIIOl). 

Actually, all of these models, are based on the peculiar characteristic of introducing new 
sources into the cosmological dynamics, while it would be preferable to develop scenarios con- 
sistent with observations without invoking further parameters or components non-testable (up to 
now) at a fundamental level. 

Moreover, it is not clear where this scalar field originates from, thus leaving a great uncertainty 
on the choice of the scalar field potential. The subtle and elusive nature of dark energy has led 
many authors to look for completely different scenarios able to give a quintessential behavior 
without the need of exotic components. To this aim, it is worth stressing that the acceleration of 



0.2 Issues from Cosmology 



xiii 



the Universe only claims for a negative pressure dominant component, but does not tell anything 
about the nature and the number of cosmic fluids filling the Universe. 

Actually, there is still a different way to face the problem of cosmic acceleration. As stressed 
in [I104L it is possible that the observed acceleration is not the manifestation of another ingredient 
in the cosmic pie, but rather the first signal of a breakdown of our understanding of the laws of 
gravitation (in the infra-red limit). 

From this point of view, it is thus tempting to modify the Friedmann equations to see whether 
it is possible to fit the astrophysical data with models comprising only the standard matter. Inter- 
esting examples of this kind are the Cardassian expansion [11051 and the Dvali-Gabadaze-Porrati 
gravity [|97l . Moving in this same framework, it is possible to find alternative schemes where a 
quintessential behavior is obtained by taking into account effective models coming from funda- 
mental physics giving rise to generalized or HOG actions |I56l [60l [59l |20l (for a comprehensive 
review see [l76ll '). 

For instance, a cosmological constant term may be recovered as a consequence of a non- 
vanishing torsion field thus leading to a model which is consistent with both SNela Hubble diagram 
and Sunyaev - Zel'dovich data coming from clusters of galaxies [I106II . SNela data could also be 
efficiently fitted including higher-order curvature invariants in the gravity Lagrangian ll58l 11071 
11081 11091 . It is worth noticing that these alternative models provide naturally a cosmological 
component with negative pressure whose origin is related to the geometry of the Universe thus 
overcoming the problems linked to the physical significance of the scalar field. 

It is evident, from this short overview, the large number of cosmological models which are 
viable candidates to explain the observed accelerated expansion. This abundance of models is, 
from one hand, the signal of the fact that we have a limited number of cosmological tests to 
discriminate among rival theories, and, from the other hand, that a urgent degeneracy problem has 
to be faced. To this aim, it is useful to remark that both the SNela Hubble diagram and the angular 
size - redshift relation of compact radio sources UllOi llll| are distance based methods to probe 
cosmological models so then systematic errors and biases could be iterated. From this point of 
view, it is interesting to search for tests based on time-dependent observables. 

For example, one can take into account the lookback time to distant objects since this quan- 
tity can discriminate among different cosmological models. The lookback time is observationally 
estimated as the difference between the present day age of the Universe and the age of a given 
object at redshift z. Such an estimate is possible if the object is a galaxy observed in more than one 
photometric band since its color is determined by its age as a consequence of stellar evolution. It 
is thus possible to get an estimate of the galaxy age by measuring its magnitude in different bands 
and then using stellar evolutionary codes to choose the model that reproduces the observed colors 
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at best. 

The resort to modified gravity theories, which extend in some way the GR, allows to pursue 
this different approach (no further unknown sources) giving rise to suitable cosmological models 
where a late time accelerated expansion naturally arises. 

The idea that the Einstein gravity should be extended or corrected at large scales (infrared limit) 
or at high energies (ultraviolet limit) is suggested by several theoretical and observational issues. 
Quantum field theories in curved spacetimes, as well as the low energy limit of string theory, both 
imply semi - classical effective Lagrangians containing higher-order curvature invariants or scalar- 
tensor terms. In addition, GR has been tested only at solar system scales while it shows several 
shortcomings if checked at higher energies or larger scales. 

Of course modifying the gravitational action asks for several fundamental challenges. These 
models can exhibit instabilities HI 12[ II 13] II 1411 or ghost -like behaviors ||32ll , while, on the other 
side, they should be matched with the low energy limit observations and experiments (solar system 
tests). Despite of all these issues, in the last years, several interesting results have been achieved 
in the framework of /- gravity at cosmological, galactic and solar system scales. 

For example, models based on generic functions of the Ricci scalar R show cosmological 
solution with late time accelerating dynamics ll20l[59l[60l[6Tl[62ll63l . in addition, it has been shown 
that some of them could agree with CMBR observational prescriptions II115II1161 . nevertheless this 
matter is still argument of debate 11117111181 . For a review of the models and their cosmological 
applications see, e.g.. ir76iril8ll . 

Moreover, considering /-gravity in the low energy limit, it is possible to obtain corrected grav- 
itational potentials capable of explaining the flat rotation curves of spiral galaxies without consid- 
ering huge amounts of dark matter IfTOl 11191 1120[ 11211 11221 and, furthermore, this seems the only 
self-consistent way to reproduce the universal rotation curve of spiral galaxies II123L On the other 
hand, several anomalies in Solar System experiments could be framed and addressed in this picture 
[fT24l[T25]l. 

Summarizing, almost 95% of matter-energy content of the universe is unknown in the frame- 
work of Standard Cosmological Model while we can experimentally probe only gravity and ordi- 
nary (baryonic and radiation) matter. Considering another point of view, anomalous acceleration 
(Solar System), dark matter (galaxies and galaxy clusters), dark energy (cosmology) could be 
nothing else but the indications that shortcomings are present in GR and gravity is an interaction 
depending on the scale. The assumption of a linear Lagrangian density in the Ricci scalar R for 
the Hilbert-Einstein action could be too simple to describe gravity at any scale and more general 
approaches should be pursued to match observations. Among these schemes, several motivations 
suggest to generalize GR by considering gravitational actions where generic functions of curvature 
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invariants are present. 

0.3 The Weak Field Limit of Higher Order Gravity 

It is well known that GR is the cornerstone theory among several attempts proposed to describe 
gravity. It represents an elegant approach giving several phenomenological predictions, but its 
validity, in the Newtonian limit regime, is experimentally probed only at Solar System scales. 
However, also at these scales, some conundrums come out as an apparent, anomalous, long-range 
acceleration revealed from the data analysis of Pioneer 10/11, Galileo, and Ulysses spacecrafts. 
Such a feature is difficult to be framed in the standard scheme of GR and its low energy limit 
II126[ 11271 . while it could be framed in a general theoretical scheme by taking corrections to the 
Newtonian potential into account [|128ll . Furthermore, at galactic distances, huge bulks of dark 
matter are needed to provide realistic models matching with observations. In this case, retaining 
GR and its low energy limit, implies the introduction of an actually unknown ingredient. We face 
a similar situation even at larger scales: clusters of galaxies are gravitationally stable and bounded 
only if large amounts of dark matter are supposed in their potential wells. 

Taking into account the weak field limit approximation, ETGs are expected to reproduce GR 
ni29l . This fact is matter of debate since several relativistic theories do not reproduce exactly the 
Einstein results in the Newtonian approximation but, in some sense, generalize them. As it was 
firstly noticed by Stelle f32], a i?^ -theory gives rise to Yukawa-like corrections in the Newtonian 
potential. Such a feature could have interesting physical consequences. For example, some authors 
claim to explain the flat rotation curves of galaxies by using such terms [11301 . Others [11311 have 
shown that a conformal theory of gravity is nothing else but a HOG model containing such terms 
in the Newtonian limit. 

In general, any relativistic theory of gravitation yields corrections to the Newton potential (see 
for example [|132|[ ') which, in the post-Newtonian formalism, could be a test for the same theory 
[|129|[ . Furthermore the newborn gravitational lensing astronomy [|133|[ is giving rise to additional 
tests of gravity over small, large, and very large scales which soon will provide direct measure- 
ments for the variation of the Newton coupling [|l34l . the potential of galaxies, clusters of galaxies 
and several other features of self-gravitating systems. Such data could be, very likely, capable of 
confirming or ruling out the physical consistency of GR or of any ETG. 

In recent papers, some authors have confronted this kind of theories even with the Post Pa- 
rameterized Newtonian (PPN) prescriptions in metric and Palatini approaches. The results seem 
controversial since in some cases [[1351 1136[[ it is argued that GR is always valid at Solar System 
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scales and there is no room for other theories; nevertheless, some other studies [11371 1138II find 
that recent experiments as Cassini and Lunar Laser Ranging allow the possibility that ETGs could 
be seriously taken into account. In particular, it is possible to define PPN-parameters in term of 
/-gravity functions and several classes of fourth order theories result compatible with experiments 
in Solar System ifTTTll . 

In a recent paper [I139II . spherically symmetric solutions for /-gravity in vacuum have been 
found considering relations among functions defining the spherical metric or imposing a constant 
Ricci curvature scalar. The authors have been able to reconstruct, at the end, the form of some 
/-theories, discussing their physical relevance. In H 14011 , the same authors have discussed static 
spherically symmetric perfect fluid solutions for /-gravity in metric formalism. They showed that 
a given matter distribution is not capable of determining the functional form of /. 

The discussion about the short scale behavior of HOG has been quite vivacious in the last years 
since OR shows is best predictions just at the Solar System level. As matter of fact, measurements 
coming from weak field limit tests like the bending of light, the perihelion shift of planets, frame 
dragging experiments represent inescapable tests for whatever theory of gravity. Actually, in our 
opinion, there are sufficient theoretical predictions to state that HOG can be compatible with New- 
tonian and post-Newtonian prescriptions [A]. In other papers II137II . it has been that this result can 
be achieved by means of the analogy of / - models with scalar - tensor gravity. 

Nevertheless, up to now, the discussion on the weak field limit of / - theories is far to be defini- 
tive and there are several papers claiming for opposite results [11251 11381 11421 11431 11441 11451 11481 
I149[|150l l. or stating that no progress has been reached in the last years due to the several common 
misconceptions in the various theories of gravity [|141[| . 

In the last few years, several authors have dealt with this matter with contrasting conclusions, 
in particular with respect to the PPN limit lfT35l [T37l [T38l [T43l [TSTTl . On the other hand, the 
investigation of spherically symmetric solutions for such kind of models has been developed in 
several papers 11391 11401 11421 1152II . Such an analysis deserves particular attention since it can 
allow to draw interesting conclusions on the effective modification of the gravitational potential 
induced by HOG at low energies and, in addition, it could shed new light on the PPN limit of such 
theories at least in a preliminary way. For example, theories like / = + which fairly address 
the cosmic acceleration issue ll56l[57l[58l . suffer a ill-defined PPN-limit since a theory containing 
terms like R^^ is singular in i? = and does not admit any Minkowski limit and then any other 
background solution which is Ricci flat. On the other hand, affirming that the unique /-gravity 
spherically symmetric "static" solution, corresponding to a realistic mass source and matching the 
present cosmic background at infinity has the PPN parameter 7 = 1/2, in conflict with experiments 
which give 7 ~ 1 could be misleading since by assuming, for example, / = R^^'' with e — > has 
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to give results compatible with GR (i.e. 7 ~ 1). 

As a matter of fact, defining a PPN-limit in such a case is a quite delicate issue, since in order 
to develop an analytical study of the deviation from the Newtonian approximation requires that 
the spacetime should be, at least, asymptotically Minkowski. Finally, approaching the problem 
considering a curvature constant metric, as in the case of the Schwarzschild - de Sitter solution, 
could induce to flawed conclusions. As a consequence, understanding the properties of HOG with 
respect to spherically symmetric solutions in the weak field limit turns out to be a key issue from 
several points of view. 

This is only one example about the debate on the weak limit field limit: authors approached the 
weak limit issue following different schemes and developing different parameterizations which, in 
some cases, turn out to be not necessarily correct. 

The purpose of this thesis (also referring to the published papers [C], [D]) is to take part to 
the debate, building up a rigorous formalism which deals with the formal definition of weak field 
and small velocities limit applied to HOG gravity. In a series of papers [C], [D], [E], the aim is 
to pursue a systematic discussion involving: i) the Newtonian limit of /-gravity, ii) spherically 
symmetric solutions toward the weak filed limit of /-gravity; and, finally. Hi) considering gen- 
eral HOG theories where also invariants as R^yR^" or Ra^^^R"'^^'^ are considered. Besides the 
Birkhoff theorem is not a general result in HOG models [C] also if it holds for several interesting 
classes of these theories as discussed, for example, in Ill53[|154l . 

The analysis is based on the metric approach, developed in the Jordan frame, assuming that 
the observations are performed in it, without resorting to any conformal transformation as done in 
several cases [I135II . This point of view is adopted in order to avoid dangerous variable changes 
which could compromise the correct physical interpretation of the results. 

A relevant aspect of HOG thoeries, in the post-Minkowskian limit is the propagation of the 
gravitational fields. It turns out that wave signals can be characterized with both tensorial and 
scalar mode [F]. This issue represents a quite striking difference between GR-like models and ex- 
tended gravity models since in the standard Einstein scheme only tensorial degrees of freedom are 
allowed. As matter of facts, the gravitational wave limit of these models can represent an interest- 
ing framework to study, in order to discuss the physical observable footprints which discriminate 
between GR and HOG experimental predictions. 

In this PhD thesis, we are going to analyze and discuss, in a general way and without specifying 
a priori the form of the Lagrangian, the relation between the spherical symmetry and the weak 
field limit, pointing out the differences and the relations with respect to the post-Newtonian and the 
post-Minkowskian limits of f -gravity. Our aim is to develop a systematic approach considering 
the theoretical prescriptions to obtain a correct weak field limit in order to point out the analogies 
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and the differences with respect to GR. A fundamental issue is to recover the asymptotically flat 
solution in absence of gravity and the well-known results related to the speciflc case f = R, i.e. 
GR. Only in this situation a correct comparison between GR and any ETG is possible from an 
experimental and a theoretical viewpoint. For example, 

In literature, there are several definitions and several claims in this direction but clear state- 
ments and discussion on these approaches urge in order to find out definite results to be tested by 
experiments II141L 

0.4 Plan of thesis 

The layout of the PhD thesis is organized as follows. In the first chapter we report a general 
review of ETGs and the fundamental aspects of GR. In particular we display all fundamental tools: 
Einstein Equation, Bianchi Identity, Conformal transformations. Metric and Palatini formalism, 
ETG theories (Scalar-tensor, HOG theories and so on). Coordinates system transformations and 
the relations between them (for example standard, isotropic coordinates etc). 

In the second chapter some "exact" spherically symmetric solution of GR is shown (Schwar- 
zschild, Schwarzschild-de Sitter, Reissner-Nordstrom solution). On the other hand, we show the 
technicality of development of field equation [C] with respect to Newtonian and Post-Newtonian 
approach: in such case we also introduce the Eddington parameters. Finally we perform the post- 
Minkowskian limit: the gravitational waves. The developments are computed in generic coordi- 
nates systems and in the gauge harmonic. 

The third chapter is devoted to some general remarks on spherical symmetry in /-gravity [D]. 
In particular, the expression of the Ricci scalar and the general form of metric components are 
derived in spherical symmetry discussing how recovering the correct Minkowski flat limit. We dis- 
cuss the spherically symmetric background solutions with constant scalar curvature considering, in 
particular Schwarzschild-like and Schwarzschild-de Sitter-like solutions with constant curvature; 
we discuss also the cases in which the spherical symmetry is present also for the Ricci scalar de- 
pending on the radial coordinate r. This is an interesting situation, not present in GR. In fact, as 
it is well known, in the Einstein theory, the Birkhoff theorem states that a spherically symmetric 
solution is always stationary and static [11571 and the Ricci scalar is constant. In / the situation 
is more general and then the Ricci scalar, in principle, can evolve with radial and time coordi- 
nates. Finally, the last part of chapter is devoted to the study of a perturbation approach starting 
from a spherically symmetric background considering the general case in which the Ricci scalar is 
R = R{r). The motivation is due to the fact that, in GR, the Schwarzschild solution and the weak 
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field limit coincide under suitable conditions. 

In the fourth one we want to seek for a general method to find out spherically symmetric 
solutions in /-gravity and, eventually, in generic ETG [B]. Asking for a certain symmetry of the 
metric, we would like to investigate if such a symmetry holds for a generic theory of gravity. In 
particular for the /-theories. Specifically, we want to apply the Noether Symmetry Approach [|26l 
in order to search for spherically symmetric solutions in generic /-theories of gravity. This means 
that we consider the spherical symmetry for the metric as a Noether symmetry and search for / 
Lagrangians compatible with it. 

In the fifth chapter, we follow a different approach. Starting from the definitions of PPN- 
parameters in term of a generic analytic function / and its derivatives, we deduce a class HOG 
theories, compatible with data, by means of an inverse procedure which allows to compare PPN- 
conditions with data [A]. As a matter of fact, it is possible to show that a third order polynomial, 
in the Ricci scalar, is compatible with observational constraints on PPN-parameters. The degree of 
deviation from OR depends on the experimental estimate of PPN-parameters. The second part of 
chapter is dedicated to very strong debate about the analogy or not between /- and Scalar-tensor 
gravity [H]. In fact for some authors the Newtonian limit of /-gravity is equivalent to the one of 
Brans-Dicke gravity with lobd = 0, so that the PPN parameters of these models turn out to be 
ill defined. We don't agree with this claim. We show that this is indeed not true. We discuss 
that HOG models are dynamically equivalent to a O'Hanlon Lagrangian which is a special case 
of Brans-Dicke theory characterized by a self-interaction potential and that, in the low energy and 
small velocity limit, this will imply a non-standard behaviour. This result turns out to be completely 
different from the one of a pure Brans-Dicke model and in particular suggests that it is completely 
misleading to consider the PPN parameters of this theory with ubd = in order to characterize 
the homologous quantities of /-gravity. 

In the sixth one we analyze the Newtonian limit of HOG theory. We are going to focus exclu- 
sively on the weak field limit within the metric approach [C], [E], [G]. At this point we remind 
the readers that it was already shown in II158II that different variational procedures do not lead to 
equivalent results in the case of quadratic order Lagrangians, casting a shadow on several newer 
works in which this equivalence was implicitly assumed. In the first part we use the development 
shown in the second chapter for a generic analytic function / and find the solution in the vacuum. 
For the sake of completeness we will treat the problem also by imposing the harmonic gauge on 
the field equations. Besides, we show that the Birkhoff theorem is not a general result for /-gravity 
since time-dependent evolution for spherically symmetric solutions can be achieved depending on 
the order of perturbations. In the second part we find again the Newtonian limit but for a so-called 
quadratic gravity lagrangian with the Green function method. We find the internal and external 
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potential generated by an extended spherically symmetric matter source. In the last part we out- 
line the general approach to find the expression of metric tensor at fourth order perturbative for a 
generic /-theory. 

In the seventh chapter we develop a formal description of the gravitational waves propagation 
in HOG models focusing on the scalar degrees of freedom and the characteristic of such scalar 
candidate in the gravity sector of gauge bosons [F]. As in the previously chapter we performed 
the Newtonian limit in vacuum with a spherically symmetric solution in standard coordinates, now 
we repeat the development but in the post-Minkowskian limit. In addition we discuss, in such 
a framework, the definition of the energy-momentum tensor of gravity which is a fundamental 
quantity in order to calculate the gravitational time delay in Pulsar timing. Some considerations on 
the differences between GR and HOG in the post-Minkowskian limit are addressed. 

Finally in the last chapter we report the discussions and conclusions. 



Chapter 1 



Extended theories of gravity: a review 



1.1 What a good theory of Gravity has to do: General Relativ- 
ity and its extensions 

From a phenomenological point of view, there are some minimal requirements that any relativistic 
theory of gravity has to match. First of all, it has to explain the astrophysical observations (e.g. the 
orbits of planets, the potential of self-gravitating structures). 

This means that it has to reproduce the Newtonian dynamics in the weak-energy limit. Besides, 
it has to pass the classical Solar System tests which are all experimentally well founded ||129I . 

As second step, it should reproduce galactic dynamics considering the observed baryonic con- 
stituents (e.g. luminous components as stars, sub-luminous components as planets, dust and gas), 
radiation and Newtonian potential which is, by assumption, extrapolated to galactic scales. 

Thirdly, it should address the problem of large scale structure (e.g. clustering of galaxies) and 
finally cosmological dynamics, which means to reproduce, in a self-consistent way, the cosmo- 
logical parameters as the expansion rate, the Hubble constant, the density parameter and so on. 
Observations and experiments, essentially, probe the standard baryonic matter, the radiation and 
an attractive overall interaction, acting at all scales and depending on distance: the gravity. 

The simplest theory which try to satisfies the above requirements was formulated by Albert 
Einstein in the years 1915 - 1916 HI 591 and it is known as the Theory of General Relativity. It is 
firstly based on the assumption that space and time have to be entangled into a single spacetime 
structure, which, in the limit of no gravitational forces, has to reproduce the Minkowski spacetime 
structure. Einstein profitted also of ideas earlier put forward by Riemann, who stated that the 
Universe should be a curved manifold and that its curvature should be established on the basis of 
astronomical observations [I160L 
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In other words, the distribution of matter has to influence point by point the local curvature of 
the spacetime structure. The theory, eventually formulated by Einstein in 1915, was strongly based 
on three assumptions that the physics of gravitation has to satisfy. 

The "Principle of Relativity" , which states that amounts to require all frames to be good frames 
for Physics, so that no preferred inertial frame should be chosen a priori (if any exist). 

The "Principle of Equivalence", that amounts to require inertial effects to be locally indis- 
tinguishable from gravitational effects (in a sense, the equivalence between the inertial and the 
gravitational mass). 

The "Principle of General Covariance" , that requires field equations to be "generally covari- 
ant" (today, we would better say to be invariant under the action of the group of all spacetime 
diffeomorphisms) [11611 . 

And - on the top of these three principles - the requirement that causality has to be preserved 
(the "Principle of Causality", i.e. that each point of spacetime should admit a universally valid 
notion of past, present and future). 

Let us also recall that the older Newtonian theory of spacetime and gravitation - that Einstein 
wanted to reproduce at least in the limit of small gravitational forces (what is called today the "post- 
Newtonian approximation") - required space and time to be absolute entities, particles moving in a 
preferred inertial frame following curved trajectories, the curvature of which (i.e., the acceleration) 
had to be determined as a function of the sources (i.e., the "forces"). 

On these bases, Einstein was led to postulate that the gravitational forces have to be expressed 
by the curvature of a metric tensor field ds'^ = gai3dx°'dx^ on a four-dimensional spacetime man- 
ifold, having the same signature of Minkowski metric, i.e., the so-called "Loren-tzian signature", 
herewith assumed to be (+, — , — , — ). He also postulated that spacetime is curved in itself and 
that its curvature is locally determined by the distribution of the sources, i.e. - being spacetime a 
continuum - by the four-dimensional generalization of what in Continuum Mechanics is called the 
"matter stress-energy tensor", i.e. a rank-two (symmetric) tensor T^^. 

Once a metric g^i, is given, the inverse g'^^ satisfies the condition 




6^ . 



(1.1) 



Its curvature is expressed by the Riemann tensor (curvature) 
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where the comas are partial derivatives. The F^^^ are the Christoffel symbols given by 



(1.3) 



if the Levi-Civita connection is assumed. The contraction of the Riemann tensor (11.21 



■per "per I per pp pp per 



(1.4) 



is the Ricci tensor and the scalar 



R = g^^R^r = R\ = g^^^i,. - g'^Ka,^ + g^^^^Kp - g^^^^L^p (i-5) 



is called the scalar curvature of g^^- The Riemann tensor (11.21) satisfies the so-called Bianchi 
identities: 



Rap.l3u;& + RapSP]U + Rap.i/S;(3 — 



Rafil3u' ~^ Rp.P\y RlJ.v\t 



2R^^'" - Rj^ -0 



(1.6) 
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where the covariant derivative is A"^-"^.^ = V^A"^-'^ = A"'^-^^ + T'^^A'''^-^ + T^^A"'''-^ + 
■■■ + ri^A'^f^-" and V„V° = □ = ^"(v^^'^aa) d'Alembert operator with respect to the 
metric g^^. 

Einstein was led to postulate the following equations for the dynamics of gravitational forces 



RpU 



(1.7) 



where X = SttG is a coupling constant (we will use the convention c = 1). These equations 
turned out to be physically and mathematically unsatisfactory. 
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As Hilbert pointed out II161II . they have not a variational origin, i.e. there was no Lagrangian 
able to reproduce them exactly (this is slightly wrong, but this remark is unessential here). Einstein 
replied that he knew that the equations were physically unsatisfactory, since they were contrasting 
with the continuity equation of any reasonable kind of matter. Assuming that matter is given as a 
perfect fluid, that is 



T^,u = {p + p)u^,u^-pgf,u (1-8) 

where u^u^ define a comoving observer, p is the pressure and p the density of the fluid, then the 
continuity equation requires T^jj, to be covariantly constant, i.e. to satisfy the conservation law 



T'^'".^ = 0. (1.9) 

In fact, it is not true that R^" vanishes (unless R = 0). Einstein and Hilbert reached in- 
dependently the conclusion that the wrong field equations (11.71) had to be replaced by the correct 
ones 



G^, = XT^, (1.10) 

where 



R 

G^y = R^y — —g^U (1-11) 

that is currently called the "Einstein tensor" of g^^,. These equations are both variational and satisfy 
the conservation laws (11.91) since the following relation holds 



G^"" = 0, (1.12) 

as a byproduct of the so-called Bianchi identities that the curvature tensor of g^^y has to satisfy 
nil [1621. 

The Lagrangian that allows to obtain the field equations (11.101 ) is the sum of a matter La- 
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grangian Cm and of the Ricci scalar: 



Che = + XCr^ 



(1.13) 



where denotes the square root of the value of the determinant of the metric g^y. The action 
of GR is 



From the action principle, we get the field equations (11.101 ) by the variation: 



5A 



-XT, 



where T^,^ is energy momentum tensor of matter: 



T, 



1 6{y/^Cr. 



The last term in (|1.15l) is a 4-divergence 



(1.14) 



Sg^'" + / d^'x^g^^SR., = , (1.15) 



(1.16) 



j dS^g^'^SR,, = j d'x^[{-5gn-,,,-U{g^''5g,,)] (1.17) 



then we can neglect it and we get the field equation (11.101) . For the variational calculus (11.151) we 
used the following relations 
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< 6R = R^pSg^^ + g^HRo,^ 



(1.18) 



{5R^p = Wg\.^,^ + 5g'^ 
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The choice of Hilbert and Einstein was completely arbitrary (as it became clear a few years 
later), but it was certainly the simplest one both from the mathematical and the physical viewpoint. 
As it was later clarified by Levi-Civita in 1919, curvature is not a "purely metric notion" but, rather, 
a notion related to the "linear connection" to which "parallel transport" and "covariant derivation" 



In a sense, this is the precursor idea of what, in the sequel, would be called a "gauge theoretical 
framework" [I164II . after the pioneering work by Cartan in 1925 II165II . But at the time of Einstein, 
only metric concepts were at hands and his solution was the only viable. 

It was later clarified that the three principles of relativity, equivalence and covariance, together 
with causality, just require that the spacetime structure has to be determined by either one or both 
of two fields, a Lorentzian metric g and a linear connection F, assumed to be torsionless for the 
sake of simplicity. 

The metric g fixes the causal structure of spacetime (the light cones) as well as its metric rela- 
tions (clocks and rods); the connection V fixes the free-fall, i.e. the locally inertial observers. They 
have, of course, to satisfy a number of compatibility relations which amount to require that pho- 
tons follow the null geodesies of F, so that F and g can be independent, a priori, but constrained, 
a posteriori, by some physical restrictions. These, however, do not impose that F has necessarily 
to be the Levi Civita connection of g II166II . 

This justifies - at least on a purely theoretical basis - the fact that one can envisage the so-called 
"alternative theories of gravitation", that we prefer to call "Extended Theories of Gravitation" 
since their starting points are exactly those considered by Einstein and Hilbert: theories in which 
gravitation is described by either a metric (the so-called "purely metric theories"), or by a linear 
connection (the so-called "purely affine theories") or by both fields (the so-called "metric-affine 
theories", also known as "first order formalism theories"). In these theories, the Lagrangian is a 
scalar density of the curvature invariants constructed out of both g and F. 

The choice (11.131) is by no means unique and it turns out that the Hilbert-Einstein Lagrangian 
is in fact the only choice that produces an invariant that is linear in second derivatives of the metric 
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(or first derivatives of the connection). A Lagrangian that, unfortunately, is rather singular from 
the Hamiltonian viewpoint, in much than same way as Lagrangians, linear in canonical momenta, 
are rather singular in Classical Mechanics (see e.g. 1116711 ). 

A number of attempts to generalize GR (and unify it to Electromagnetism) along these lines 
were followed by Einstein himself and many others (Eddington, Weyl, Schrodinger, just to quote 
the main contributors; see, e.g., II168II ) but they were eventually given up in the fifties of XX 
Century, mainly because of a number of difficulties related to the definitely more complicated 
structure of a non-linear theory (where by "non-linear" we mean here a theory that is based on non- 
linear invariants of the curvature tensor), and also because of the new understanding of physics that 
is currently based on four fundamental forces and requires the more general "gauge framework" to 
be adopted (see 1116911 ). 

Still a number of sporadic investigations about "alternative theories" continued even after 1960 
(see [11291 and refs. quoted therein for a short history). The search of a coherent quantum theory 
of gravitation or the belief that gravity has to be considered as a sort of low-energy limit of string 
theories (see, e.g., 1117011 ) - something that we are not willing to enter here in detail - has more or 
less recently revitalized the idea that there is no reason to follow the simple prescription of Einstein 
and Hilbert and to assume that gravity should be classically governed by a Lagrangian linear in the 
curvature. 

Further curvature invariants or non-linear functions of them should be also considered, espe- 
cially in view of the fact that they have to be included in both the semi-classical expansion of a 
quantum Lagrangian or in the low-energy limit of a string Lagrangian. 

Moreover, it is clear from the recent astrophysical observations and from the current cosmo- 
logical hypotheses that Einstein equations are no longer a good test for gravitation at Solar Sys- 
tem, galactic, extra-galactic and cosmic scale, unless one does not admit that the matter side of 
Eqs. (|1.10l) contains some kind of exotic matter-energy which is the "dark matter" and "dark en- 
ergy" side of the Universe. 

The idea which we propose here is much simpler. Instead of changing the matter side of Ein- 
stein Equations (11.101) in order to fit the "missing matter-energy" content of the currently observed 
Universe (up to the 95% of the total amount!), by adding any sort of inexplicable and strangely 
behaving matter and energy, we claim that it is simpler and more convenient to change the gravita- 
tional side of the equations, admitting corrections coming from non-linearities in the Lagrangian. 
However, this is nothing else but a matter of taste and, since it is possible, such an approach should 
be explored. Of course, provided that the Lagrangian can be conveniently tuned up (i.e., chosen 
in a huge family of allowed Lagrangians) on the basis of its best fit with all possible observational 
tests, at all scales (solar, galactic, extragalactic and cosmic). 
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Something that - in spite of some commonly accepted but disguised opinion - can and should 
be done before rejecting a priori a non-linear theory of gravitation (based on a non-singular La- 
grangian) and insisting that the Universe has to be necessarily described by a rather singular grav- 
itational Lagrangian (one that does not allow a coherent perturbation theory from a good Hamil- 
tonian viewpoint) accompanied by matter that does not follow the behavior that standard baryonic 
matter, probed in our laboratories, usually satisfies. 



1.2 The Extended Theories of Gravity: F{R, BR,..., B^R, 0) 



With the above considerations in mind, let us start with a general class of higher-order-scalar-tensor 
theories in four dimensions given by the action 



where F is an unspecified function of curvature invariants and of a scalar field (f). The term Cm, as 
above, is the minimally coupled ordinary matter contribution and is a generic function of the 
scalar field (j). For example its values could be uj{(p) = ±1,0 fixing the nature and the dynamics 
of the scalar field which can be a standard scalar field, a phantom field or a field without dynamics 
(see [El [ml for details). 

In the metric approach, the field equations are obtained by varying (|1.19l) with respect to g^u- 
We get 




(1.19) 
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(1.20) 



where G^u is the above Einstein tensor (11.111 ) and 
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The differential Equations (11.201) are of order (2k + 4). The (eventual) contribution of a potential 
V{(f)) is contained in the definition of By varying with respect to the scalar field 0, we obtain 
the Klein - Gordon equation 



^ + ■ - ^' . (1.22, 

2 (30 2c<;((;i)) ()0 

Several approaches can be used to deal with such equations. For example, as we said, by a 
conformal transformation, it is possible to reduce an ETG to a (multi) scalar - tensor theory of 
gravity [|23l HOI HI [ml [ml . 

The simplest extension of GR is achieved assuming 



r = f{R), a;(0) = O, (1.23) 
in the action (11.191 ); / is an arbitrary (analytic) function of the Ricci curvature scalar R. Then 



(1.24) 



where the standard Hilbert-Einstein action is, of course, recovered for f = R. Varying the (|1.24l) 
with respect to gfj,,^, we get the field equations 



H,. = f'R,. - \fg,. - /;^, + g,.Uf = XT,,, (1.25) 

which are fourth-order equations due to the terms f'^^ and □/'; the prime indicates the derivative 
with respect to R. The trace of (11.251) is 



H = g'^^H^^ = 3nf + fR-2f = XT. (1.26) 
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The peculiar behavior of / = i? is due to the particular form of the Lagrangian itself which, 
even though it is a second order Lagrangian, can be non-covariantly rewritten as the sum of a first 
order Lagrangian plus a pure divergence term. The Hilbert-Einstein Lagrangian can be in fact 
recast as follows: 



i^uE = V^R = v^^7"^(rLr^^ - r^.r^^) + "^AV^g'^^u^^) ; (1.27) 

r is the Levi - Civita connection of g and is a quantity constructed out with the variation of 
r [HI I162L Since is not a tensor, the above expression is not covariant; however a standard 
procedure has been studied to recast covariance in the first order theories [I173II . This clearly shows 
that the field equations should consequently be second order and the Hilbert-Einstein Lagrangian 
is thus degenerate. 

From the action (11.191 ). it is possible to obtain another interesting case by choosing 



r = + (1.28) 

where V{(f)) and are generic functions describing respectively the potential and the coupling 
of a scalar field 0. In this case, we get 



A^'^ = j d^x^f^[F{(t))R + + 1/(0) + XCr 



(1.29) 



The Brans-Dicke theory of gravity is a particular case of the action (11.291) in which we have Vi^cj)) 
and a;(0) = — In fact we have 



A 



BD 



d'^x^J—g 



(f)R — ujbd^—: 1- '^C,r 



(1.30) 



The variation of (11.291) with respect to g^y and gives the second-order field equations 
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(1.31) 



The third equation in (11.311) is the trace of field equation for g^y and the last one is a combination 
of the trace and of the one for cf). This last equation is equivalent to the Bianchi contracted identity 
II174II . Standard fluid matter can be treated as above. 



1.3 Conformal transformations 

Let us now introduce conformal transformations to show that any higher-order or scalar-tensor 
theory, in absence of ordinary matter, e.g. a perfect fluid, is conformally equivalent to an Einstein 
theory plus minimally coupled scalar fields. If standard matter is present, conformal transfor- 
mations allow to transfer non-minimal coupling to the matter component ni75B . The conformal 
transformation on the metric g^^ is 



with A{x^) > 0. A is the conformal factor. Obviously the transformation rule for the contravariant 
metric tensor is g^^ = A^^g^^ . The various mathematical quantities in the so-called Einstein 
frame (EE) (quantities referred to g^u) are linked to the ones in the so-called Jordan Erame (JE) 
(quantities referred to g^^) as follows 



9i^u = A{x )g^ 



(1.32) 
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where = In A^^'^. But we can have also the inverse relations 

' = - 0,^5° - 0,,5^ + 0~«^^, 

= + ^^(0;m. + 0;m0;-) " '^"(0;m/3 + Mp) + 

< R,^u = Ri,v + 20~^^ + 20;^0;^ + ~g^y 00 - 2^^^0;-0;, (1.34) 
R = e'^^{R + 600 - 60''^0;^) 

^ G^,y = G^y + 20.^.^ + 20.^0.^, - 2g^,y □0 + g^^cp'^^cp-^ 

where □ and □ are the d'Alembert operators with respect to the metric g^i, and (jf^u- The transfor- 
mation between the operators is □ = e^'^D — 20''^(9i.. 

Under these transformations, the action in (|1.29l) can be reformulated as follows 



1.3 Conformal transformations 13 



AST 



AR + n{ip)ip,aip''' + W{ip) + Xlrr 



(1.35) 



in which R is the Ricci scalar relative to the metric g and A is a generic constant. The relations 
between the quantities in two frames are 



F{4>) 2 \ 64 



2n 



A2 



r 



F(0(¥>))2^™1 F(0(^)) 



F{(^)A{x^y^ = A 
The field equations for the new fields and if are 



(1.36) 



(1.37) 



Therefore, every non-minimally coupled scalar-tensor theory, in absence of ordinary matter, 
e.g. perfect fluid, is conformally equivalent to an Einstein theory, being the conformal transforma- 
tion and the potential suitably defined by (11.361) . The converse is also true: for a given -F(0), such 
that is valid the relations (|1.36l) . we can transform a standard Einstein theory into a non-minimally 
coupled scalar-tensor theory. This means that, in principle, if we are able to solve the field equa- 
tions in the framework of the Einstein theory in presence of a scalar field with a given potential, we 
should be able to get the solutions for the scalar-tensor theories, assigned by the coupling -F(0), 
via the conformal transformation (11.321 ) with the constraints given by (11.361 ). Following the stan- 
dard terminology, the "Einstein frame" is the framework of the Einstein theory with the minimal 
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coupling and the "Jordan frame" is the framework of the non-minimally coupled theory [I175II . 

This procedure can be extended to more general theories. If the theory is assumed to be higher 
than fourth order, we may have Lagrangian densities of the form [|40lll58ll . 



£ = (1.38) 

Every □ operator introduces two further terms of derivation into the field equations. For example, 
a theory like 



-g RUR , 



(1.39) 



is a sixth-order theory and the above approach can be pursued by considering a conformal factor 
of the form 



A 



d{RUR) ^^d{RUR) 



OR 



dDR 



(1.40) 



In general, increasing two orders of derivation in the field equations (i.e. for every term OR), 
corresponds to adding a scalar field in the conformally transformed frame [40]. A sixth-order 
theory can be reduced to an Einstein theory with two minimally coupled scalar fields; a 2n-order 
theory can be, in principle, reduced to an Einstein theory plus {n — 1) scalar fields. On the other 
hand, these considerations can be directly generalized to higher - order - scalar - tensor theories in 
any number of dimensions as shown in HtI . 

The analogy between scalar-tensor gravity and HOG, although mathematically straightforward, 
requires a careful physical analysis. Recasting fourth - order gravity as a scalar - tensor theory, 
often the following steps, in terms of a generic scalar field 0, are considered 



f+c^^ nm + - nm + ^(0) + Cm , (1.41) 

where, by analogy, (p R and the "potential" is V{(f)) = f{(f)) — f'{(f))(f). Clearly the kinetic term 
is not present so that (11.411 ) is usually referred as a scalar-tensor description of / - gravity where 
co{(p) = 0. This is the so-called O'Hanlon Lagrangian [1331 : 
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^HOG 

X 

HOGEqs. 

X 



< — > 



< — )■ 



CsT < > ^HE + 

X X 

STEqs. i — Einstein Eqs. +v? 

X X 



HOG Solutions < — > ST Solutions i — > Einstein Solutions 



Table 1.1: Summary of the three approaches: Scalar-Tensor (ST), Einstein +ip, and / and their relations at 
Lagrangians, field equations and solutions levels. The solutions are in the Einstein frame for the minimally 
coupled case while they are in Jordan frame for / and ST - gravity. Clearly, / and ST theories can be 
rigorously compared only recasting them in the Einstein frame. 



A 



OH 



(PR + V{4>) + A 



(1.42) 



As concluding remarks, we can say that conformal transformations work at three levels: i) 
on the Lagrangian of the given theory; ii) on the field equations; Hi) on the solutions. The table 
11.11 summarizes the situation for HOG, non-minimally coupled scalar-tensor theories (ST) and 
standard Hilbert-Einstein (HE) theory. Clearly, direct and inverse transformations correlate all the 
steps of the table but no absolute criterion, at this point of the discussion, is able to select which 
is the "physical" framework since, at least from a mathematical point of view, all the frames are 
equivalent [11751 . 



However, the typical Brans-Dicke action is the (ll.30h where no scalar field potential is present 
and cobd is a constant, while the O'Hanlon Lagrangian (11.421) has a potential but has no kinetic 
term. The most general situation is in (11.291) where we have non-minimal coupling, kinetic term, 
and scalar field potential. This means that fourth-order gravity and scalar tensor gravity can be 
"compared" only by means of conformal transformations where kinetic and potential terms are 
preserved. In particular, it is misleading to state that PPN - limit of HOG is not working since 
these models provide lobd = and this is in contrast with observations [I1351I136I . 



Scalar-tensor theories and /-theories can be rigorously compared, after conformal transforma- 
tions, in the Einstein frame where both kinetic and potential terms are present. 
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1.4 The Palatini Approach and the Intrinsic Conformal Struc- 



As we said, the Palatini approach, considering g and T as independent fields, is "intrinsically" 
bi-metric and capable of disentangling the geodesic structure from the chronological structure of 
a given manifold. Starting from these considerations, conformal transformations assume a funda- 
mental role in defining the affine connection which is merely "Levi - Civita" only for the Hilbert- 
Einstein theory. 

In this section, we work out examples showing how conformal transformations assume a fun- 
damental physical role in relation to the Palatini approach in ETGs [I176L 

Let us start from the case of fourth-order gravity where Palatini variational principle is straight- 
forward in showing the differences with Hilbert- Einstein variational principle, involving only met- 
ric. Besides, cosmological applications of /-gravity have shown the relevance of Palatini formal- 
ism, giving physically interesting results with singularity - free solutions |l20l [64l [Ml [Ml [671 [Ml 
|69|| . This last nice feature is not present in the standard metric approach. 

An important remark is in order at this point. The Ricci scalar entering in / is i? = R{g , F) = 
g°'^Rai3{T) that is a generalized Ricci scalar and R^ui^) is the Ricci tensor of a torsion-less con- 
nection r, which, a priori, has no relations with the metric g of spacetime. The gravitational part of 
the Lagrangian is controlled by a given real analytical function of one real variable /, while 
denotes a related scalar density of weight 1. Field equations, deriving from the Palatini variational 
principle are: 



where V'^ is the covariant derivative with respect to F. We shall use the standard notation denoting 
by the symmetric part of R^^, i-e. R(^u) = liR^lu + Rufi)- 

In order to get the first one of (11.431) . one has to additionally assume that is functionally 

g 

independent of F; however it may contain metric covariant derivatives V of fields. This means that 
the matter stress-energy tensor T^^j = T^„{g, ^) depends on the metric g and some matter fields 
denoted here by ^, together with their derivatives (covariant derivatives with respect to the Levi 
- Civita connection of g). From the second one of (11.431 ) one sees that ^/—gf'g^" is a symmetric 




twice contravariant tensor density of weight L As previously discussed in 11176111771 , this naturally 



ture 




(1.43) 
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leads to define a new metric h^^, such that the following relation holds: 



v/^/y^ = v^/^^^ (1.44) 

This ansatz is suitably made in order to impose V to be the Levi - Civita connection of h and 
the only restriction is that \/^f'g^'^ should be non-degenerate. In the case of Hilbert-Einstein 
Lagrangian, it is /' = 1 and the statement is trivial. 

Eg. (| 1.441) imposes that the two metrics h and g are conformally equivalent. The corresponding 
conformal factor can be easily found to be /' (in dim M. = A) and the conformal transformation 
results to be ruled by: 



V = f'gt^v (1.45) 

Therefore, as it is well known, Eg. (1 1.431) implies that V = TLc{h) and -R(^,^)(r) = R^^{h) = 
Field equations can be supplemented by the scalar-valued equation obtained by taking the trace of 
(fL43]) 



f'R-2f = XT (1.46) 

which controls solutions of (11.431) . 

We shall refer to this scalar- valued equation as the structural equation of the spacetime. In the 
vacuum case (and spacetimes filled with radiation, such that T = 0) this scalar-valued equation 
admits constant solutions, which are different from zero only if one add a cosmological constant. 
This means that the universality of Einstein field equations holds HIVVL corresponding to a theory 
with cosmological constant [|83l . 

In the case of interaction with matter fields, the structural equation (11.451) . if explicitly solvable, 
provides an expression of i? = R{T) and consequently both / and /' can be expressed in terms 
of T. The matter content of spacetime thus rules the bi-metric structure of spacetime and, con- 
sequently, both the geodesic and metric structures which are intrinsically different. This behavior 
generalizes the vacuum case and corresponds to the case of a time- varying cosmological constant. 
In other words, due to these features, conformal transformations, which allow to pass from a metric 
structure to another one, acquire an intrinsic physical meaning since "select" metric and geodesic 
structures which, for a given ETG, in principle, do not coincide. 
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Let us now try to extend the above formalism to the case of non-minimally coupled scalar- 
tensor theories. The effort is to understand if and how the bi-metric structure of spacetime behaves 
in this cases and which could be its geometric and physical interpretation. 

We start by considering scalar-tensor theories in the Palatini formalism, calling Ai the action 
functional. After, we take into account the case of decoupled non-minimal interaction between a 
scalar-tensor theory and a /-theory, calling A2 this action functional. We finally consider the case 
of non-minimal-coupled interaction between the scalar field and the gravitational fields {g, T), 
calling ^3 the corresponding action functional. Particularly significant is, in this case, the limit of 
low curvature R. This resembles the physical relevant case of present values of curvatures of the 
Universe and it is important for cosmological applications. 

The action (11.291 ) for scalar-tensor gravity can be generalized, in order to better develop the 
Palatini approach, as: 

^1 = J d^Xy^[F{(l))R + v% + ^(0) + A'An(^, V ^)] (1.47) 

As above, the values of t<j(0) = ±1 selects between standard scalar field theories and quintessence 
(phantom) field theories. The relative "signature" can be selected by conformal transformations. 
Field equations for the gravitational part of the action are, respectively for the metric g and the 
connection F: 

F(0)[i?(^,) - f (?^,] = XT^, + iu;(0) v'' 0(7;.. + ^Vi(t))g^, 

(1.48) 

Vl{V^Fi^)gn = 
R(fj.u) is the same defined in (|1.43l) . For matter fields we have the following field equations: 



2u;(0) + cj,<^(0) Va V + V^{^) + F^{(j))R = 



5* ~ ^ 



(1.49) 



In this case, the structural equation of spacetime implies that: 
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P A'T + 2a;(0) V„ V° + 2^(0) 

^ = m ^^-^^^ 

which expresses the value of the Ricci scalar curvature in terms of the traces of the stress-energy 
tensors of standard matter and scalar field (we have to require 7^ 0). The bi-metric structure 
of spacetime is thus defined by the ansatz: 



-gFi(P)g 



(1.51) 



such that g and h result to be conformally related 



V = F{(t))9i^^ (1-52) 

The conformal factor is exactly the interaction factor. From (11.501) . it follows that in the vacuum 
case T = and a;((^) v 0+^(0) = 0: this theory is equivalent to the standard Einstein one 
without matter. On the other hand, for = Fq we recover the Einstein theory plus a minimally 
coupled scalar field: this means that the Palatini approach intrinsically gives rise to the conformal 
structure (11.521 ) of the theory which is trivial in the Einstein, minimally coupled case. 

As a further step, let us generalize the previous results considering the case of a non-minimal 
coupling in the framework of /-theories. The action functional can be written as: 



A2 = j d'x^[F{^)f{R) + a;(0) V + V{(P) + 2A'£„(^, V ^)] (1.53) 

where / is, as usual, any analytical function of the Ricci scalar R. Field equations (in the Palatini 
formalism) for the gravitational part of the action are: 



9 , 9 



a 



111/ 



(1.54) 

For scalar and matter fields we have, otherwise, the following field equations: 
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(1.55) 



where the non-minimal interaction term enters into the modified Klein-Gordon equations. In this 
case the structural equation of spacetime implies that: 



fR-2f = j;^^ (1.56) 

We remark again that this equation, if solved, expresses the value of the Ricci scalar curvature 
in terms of traces of the stress-energy tensors of standard matter and scalar field (we have to require 
again that F((f)) ^ 0). The bi-metric structure of spacetime is thus defined by the ansatz: 



^F{<p)fg^^ = v^/i^^ (1.57) 
such that g and h result to be conformally related by: 

V = F{<P)fg^, (1.58) 

Once the structural equation is solved, the conformal factor depends on the values of the matter 
fields (0, \E') or, more precisely, on the traces of the stress-energy tensors and the value of (p. From 

g g » 

equation (11.561) , it follows that in the vacuum case, i.e. both T = and iu{(j)) Va V + ^(0) = 
0, the universality of Einstein field equations still holds as in the case of minimally interacting /- 
theories II177I . The validity of this property is related to the decoupling of the scalar field and the 
gravitational field. 

Let us finally consider the case where the gravitational Lagrangian is a general function of 
and R. The action functional can thus be written as: 

A^ = J d^x^[K{<j) , R) + io{<j)) v. v" + V{<j)) + XCmC^, V ^)] (1-59) 
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Field equations for the gravitational part of the action are: 



' dK{4>,R) 



dR ^if^") 



K{4>,R) 



, S 9 " 1 / X 



^ a 



(1.60) 



For matter fields, we have: 



(1.61) 



SCm fl 

5* ^ 



The structural equation of spacetime can be expressed as: 



dK(6 R) 9 9°" 

^^-^i?-2ir(0,i?) = A'r + 2u;(0) Va0V + 2V(0) (1.62) 

This equation, if solved, expresses again the form of the Ricci scalar curvature in terms of traces of 
the stress-energy tensors of matter and scalar field (we have to impose regularity conditions and, 
for example, K{(j), R) ^ 0). The bi-metric structure of spacetime is thus defined by the ansatz: 



= ^/.- (1.63) 

such that g and h result to be conformally related by 



dK{(P,R) ^^^^^ 
V= 9,u. (1.64) 

Again, once the structural equation is solved, the conformal factor depends just on the values of 
the matter fields and (the trace of) their stress energy tensors. In other words, the evolution, the 
definition of the conformal factor and the bi-metric structure is ruled by the values of traces of 
the stress-energy tensors and by the value of the scalar field </>. In this case, the universality of 
Einstein field equations does not hold anymore in general. This is evident from (11.621) where the 
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strong coupling between R and (f) avoids the possibility, also in the vacuum case, to achieve simple 
constant solutions. 

We consider, furthermore, the case of small values of R, corresponding to small curvature 
spacetimes. This limit represents, as a good approximation, the present epoch of the observed Uni- 
verse under suitably regularity conditions. A Taylor expansion of the analytical function K^cp, R) 
can be performed: 



K{(f),R) = Ko{(P) + K^{(P)R + o{R^) (1.65) 
where only the first leading term in R is considered and we have defined: 



' Ko{<f)) = K{(f),R)n=o 

< (1.66) 

KM = (™) 

I. \ / _R=0 

Substituting this expression in (11.621) and (| 1.641) we get (neglecting higher order approximations 
in R) the structural equation and the bi-metric structure in this particular case. From the structural 
equation, we get: 



P T + 2u{<f)) V, v" + 2V{<f)) + 2Ko{<f)) 

= im ^ -''^ 

such that the value of the Ricci scalar is always determined, in this first order approximation, in 
terms of io{(f)) Va V + V{(f)), T e (f). The bi-metric structure is, otherwise, simply defined by 
means of the first term of the Taylor expansion, which is 



h^u = Ki{(P)g^, . (1.68) 

It reproduces, as expected, the scalar-tensor case (11.521) . In other words, scalar-tensor theories 
can be recovered in a first order approximation of a general theory where gravity and non-minimal 
couplings are any (compare (11.671 ) with (11.501 )). This fact agrees with the above considerations 
where Lagrangians of physical interactions can be considered as stochastic functions with local 
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gauge invariance properties [|44ll . 

Finally we have to say that there are also bi-metric theories which cannot be conformally related 
(see for example the summary of alternative theories given in [I129II ) and torsion field should be 
taken into account, if one wants to consider the most general viewpoint [fT3lll78ll . We will not take 
into account these general theories in this review. 

After this short review of ETGs in metric and Palatini approach, we are going to face some 
remarkable applications to cosmology and astrophysics. In particular, we deal with the straightfor- 
ward generalization of GR, the /-gravity, showing that, in principle, no further ingredient, a part a 
generalized gravity, could be necessary to address issues as missing matter (dark matter) and cos- 
mic acceleration (dark energy). However what we are going to consider here are nothing else but 
toy models which are not able to fit the whole expansion history, the structure growth law and the 
CMB anisotropy and polarization. These issues require more detailed theories which, up to now, 
are not available but what we are discussing could be a useful working paradigm as soon as refined 
experimental tests to probe such theories will be proposed and pursued. In particular, we will out- 
line an independent test, based on the stochastic background of gravitational waves, which could 
be extremely useful to discriminate between ETGs and GR or among the ETGs themselves. In this 
latter case, the data delivered from ground-based interferometers, like VIRGO and LIGO, or the 
forthcoming space interferometer LISA, could be of extreme relevance in such a discrimination. 

Finally, we do not take into account the well known inflationary models based on ETGs (e.g. 
Il45l ) since we want to show that also the last cosmological epochs, directly related to the so called 
Precision Cosmology, can be framed in such a new "economic" scheme. 



1.5 The general / — theory 



Let / be an analytic function of Ricci scalar R. We can formulate a HOG starting from the action 
principle (11.241) . By varying the action (| 1.241) and by using the properties (11.181) we get the field 
equations: 
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5A 



+n-{5gn;,.-^{9"'59,.)] 

J d^x^{H^, -XT^,)6g^^ = 



^9 



(1.69) 



where the symbol ~ means that we neglected a pure divergence; then we obtain the field equation 
(11.251) . Eq. (11.251) satisfies the condition if"^.^ = X T"''.^ = 0. In fact it is easy to check that 



fR'^'R-^a- f'"Ra^ = 

fR^'^R-^a - fR^Ra^ = ; (1.70) 

where we used the properties G°'^ = and [V^, Va]/''° = —f'''°'RJ^. If we develop the 
covariant derivatives in (|1.25l) and in (11.261) we obtain the complete expression for a generic /- 
theory 



(1.71) 

E = f'R-2j \n = XT 



where the two quantities 'K^v and 'K read : 
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U 



fii/ 



H = 3/" 



+ 3/'" 



(1.72) 



are the standard Christoffel's symbols defined by (11.31) . We conclude, then, this paragraph 



having shown the most general expression of field equations of /-gravity in metric formalism. 



1.6 The field equations for the i^^/^i?"^ and Rap^sR""^^^ - in- 
variants 



The technicality is ever the same one. We start from the action principle for a Lagrangian densities 
RaisR"^ and R^p^sR"^'^^ and we get their field equations: 
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RapR""^ 



2R^ Rau 2 



+2R^"'6R 



d'^x^J—g 



Rau 



+R'^''{2gP'^Sg,^^,,y -DSg^, - g^^^Sg,^., 



j d'x^g 



R R'^^ 
2R^°'Rau ^ giiu — XTf 



+ 



d^x\/—g 



9 p « p 



RapR"^ 



+^Rij.u + diivR 



5^^'^ = 0. (1.73) 



Then, the field equations are 



^iil"" — 2-R^"-Rai. ^ gii.v - '^R" (f_,-u)a + + gfiuR''^ -err " (1-V4) 

and the trace is 



H^^^ = 2nR = XT, (1.75) 
where we used the Bianchi identity contracted (11.61) . 

Let us calculate the field equations for the Ra/s-ysR"'^"'^ - invariant: 



5 A 
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5 j d^Xy/^lRap^sR"^^^ + XCm] = 

6 j d^xy/^lR^^.sg^'g^^g^^g'^Rpar^ + XC„ 



AT} p a/37 _ -"'"^7^ 



rCnB'vfi-tl 



g>.u-XT^.]Sg^'' + 



d'^Xyf—g 



^R^ajB'yRi 



D DO/975 



gf.u-XT^u]Sg^'' + 



d^x\/—g 



a7 



^gaP'-yS - Sga-y^pS + ^gp 7;ao / 



OP p a/37 -"-a/37<5-ft- /I P "/^ i 



-XT, 



111/ 



We used the expressions 



(1.76) 



S RaPyS = ^{gaaR'^ R-ys) = R-yS^ gaa + gaa^ 



13-yS 



dR" - ^' 



isjs - 2i^g%sj + ^g%p^ - ^gps'"^ - ^gy^s - ^g\,ps + ^gpi'^s) 

Then, the field equations, from (|1.76l) . are 



(1.77) 



H. 



Rie 



p pa/375 
OR p a/37 -^a/375-n. 



fl l/'jCip 



XT, 



flU 



(1.78) 



and the trace is 



H 



Rie 



-AR, 



XT. 



(1.79) 
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1.7 Generalities on spherical symmetry 

Since we are interesting to understand the modifications of predictions of GR when one considers 
a concentration of matter in the space, it is fundamental requiring particular properties of metric 
g^,y. The first step, and also the easiest, we will study, in the next chapters, the gravitational 
potential generated by spherically symmetric matter distribution (point-like and not) and the choice 
of mathematical form of metric becomes very important. Starting from the matter spherically 
symmetric we expect also the metric has the same symmetries. 

We conclude this chapter showing the principal relations between some coordinates systems 
we will use in this PhD thesis. 

The most general spherically symmetric metricS can be written as follows : 

ds'^ = g^it, |x|) dt^ + g2{t, |x|) dtx-dx + g^{t, |x|)(x ■ duf + gi{t, \s.\)d\s.\^ (1.80) 

where gi are functions of the distance |x| and of the time t. The set of coordinates is = 
{t, x^, x^, x^). The scalar product is defined as usual form: x ■ dx = x^dx^ + x'^dx'^ + x^dx^. By 
spherically symmetric form of (11.801) it is convenient to replace x with spherical polar coordinates 
r, 9, (J) defined as usual by 

x^ = rsin^cosc/), = rsin^sin*^, x^ = rcoscp. (1-81) 
The proper time interval (|1.80l) then becomes 

ds^ = gi{t,r) dt^ + rg2{t,r) dtdr + r'^g3{t,r)dr^ + g4{t,r){dr^ + r^dn) , (1.82) 

where d^l = dO"^ + sin^ Odcjy^ is the solid angle. We are free to reset our clocks by defining the time 
coordinate 

t = t' + C(t',r), (1.83) 

'The metric is spherically symmetric if it depends only on x and dx only through the rotational invariants ds?, 
X • dx and x^ . 
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with ({t', r) an arbitrary function of t' and r. This allows us to eliminate the off-diagonal element 
gtr in the metric (11.821) by setting 



dCit',r) ^ rg2{t',r) 
dr 2gi{t',r) 



(1.84) 



the metric (11.821) becomes 



ds'^ = gi{t',r) 



1 + 



dat',r) 
dt' 



dt'^ + 



dr^ + 



+gi{t' ,r)r^dVL, (1.85) 

where if we introduce a new metric coefficients gtt{t' r), grr{t' r) and gnnit', r) we can recast the 
(11.851) as follows 



ds = gu(t',r)dt' - grrit' ,r)dr - g^nit' ,r)dQ : 



(1.86) 



if we introduce a new radial coordinate (r') by considering a further transformation 



r' = {const)e 



(1.87) 



it is possible to recast Eq. (|1.86l) into the isotropic form (isotropic coordinates) 



ds^ = gtt{t',r')dr - gij{t',r')dx'dx^ 



(1.88) 



and then it is possible also to choose gw{t\r) = r"^ (this condition allows us to obtain the 
standard definition of the circumference with radius r") and to have the metric (|1.86l) in the standard 
form (standard coordinates) 



ds' = gu{t', r") dt' - grr{t', r")dr"^ - r"^dn . 



(1.89) 
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Obviously the functions (7ft (t', r") and grrit' , r") are not the same of (1 1 .861) . If we suppose (7ij(t', r') 
r')Sij we note that it is possible pass from (11.881) to (11.891) by the coordinate transformations: 



r' = r'{r") = (const) e^^ 



(1.90) 



We can, then, affirm that the expressions (|1.86l) . (11.881) and (11.891) are equivalent to the metric (11.801) 
and we can consider them without loss of generality as the most general definitions of a spherically 
symmetric metric compatible with a pseudo - Riemannian manifold without torsion. The choice of 
the form of the metric is only a practical issue. With this hypothesis, by inserting these metrics 
into the field equations (11.251) . one obtains: 



H, 



+g'' In Rt + (^g^\, + g'^ In ^ R^r + g''R,tt+ 



+g"'R, 



-f" 



R,tiR,u — g^u [ 9 R,t + g^^R,r 



H = f'R-2f + 3f" 
+g''R,u + g''^R 



(g'\, + g'' In ^ R^, + (^g^^^^ + In ^^^^ R, 



+ 3f" 



g''R/ + g"R/ 



(1.91) 



Eqs. (11.911 ) are the starting-point for the next chapter. All our studies in the next chapters are 
referred ever to field equations (11.911) . except the second part of sixth chapter where we have to 
insert in the field equations also the contribution of -invariant (11.741) . 

We conclude having shown the most general spherically symmetric metric tensor for our aim 
and before starting from third chapter with a systematic study of /-gravity we want to stop, in the 
second chapter, to consider the principal spherically symmetric solutions in GR. Some of these 
solutions will be the starting-point to find, with perturbative methods, the corrections induced by 
/-theory. 



Chapter 2 



Exact and perturbative solutions in General 
Relativity 

In this chapter we show, starting from knowledge of outcomes of GR, the mathematical tools 
needed for aims of present thesis. First of all we present the particular spherically symmetric 
solutions in GR and consequent Birkhoff theorem (§ 12.11) . In § 12.31 we show the technicality of 
development of field equations with respect to Newtonian and Post-Newtonian approach [C]. Fi- 
nally, in ^ |2.4| we perform the post-Minkowskian limit: the gravitational waves. The developments 
are computed in generic coordinates systems and in the gauge harmonic. 

The /-gravity theory, from mathematical point-view, is more complicated than GR. Giving the 
exact solutions of Eqs. (11.251) is vary hard challenge. Nevertheless, known the basic solutions 
of GR, we can try to find new solutions by requiring the Newtonian and post-Newtonian limit 
approach. This approach is very useful when we consider the astrophysical problems or the study 
of planet motion in the Solar System. An another field of comparison between GR and /-gravity is 
possible in the post-Minkowskian regime. In this case we can study the propagation of gravitational 
field induced by /-gravity. 

We dedicate, then, this chapter to understanding the outcomes of GR and to showing all math- 
ematical tools needed for the next chapters. 
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2.1 The Schwarzschild, Schwarzschild — de Sitter and Reiss- 
ner — Nordstrom solutions: the Birkhoff theorem in Gen- 
eral Relativity 

We can rewrite the metric (11.891 ) as follows 

ds^ = e'^^*'^) dt^ - e''^''''Ur^ - r^dQ , (2.1) 

where we recalled the radial coordinate. The only nonvanishing components of metric tensor g^^^, 
are 

9u = e^^*''-) , grr = -e'^^*''-) , gee = -r' , g^^ = -r^ sin^ 6 (2.2) 

with functions /^(t, r) and z/(t, r) that are to be determined by solving the field equations in GR 
(11.101) . Since g^i^ is diagonal, it is easy to write down all the nonvanishing components of its 
inverse: 

Furthermore, the determinant of the metric tensor is 



g = _eMt,0+K*,r)^4 gj^2 ^ (2.4) 
so the invariant volume element is 



-g dr dO ( 



Q ;^(t,r) + t^(t,r) 

r e 2 sin 9 drdOt 



(2.5) 



The only nonvanishing components of symbols Christoffel (11.31 ) are 
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tt — 2 ' 



p4 _ fi'{t,r) n(t,r)-u(t,r) 



-r e 



2 ' 



■pt _ t/'(f,r) 
^ tr 2 ' 



r-e ^ r(f> r ' 



[ = -sin^cos^, 

and the field equations (|1.10|) become 



rL = cot 6^ . 



^ tt — 2 ^ ' 



-r e 



(2.6) 



_^ _ M'(f,r) 



AT* 



7^2 



L.(t,r) 



/i(t,r) + 

p-(i(t,r) 



^'^{t,r) fi{t,T) L'{t,r) 



+ 



u"{t,r) 



2 r 



2 



(2.7) 



and if we suppose a tensor of matter like T^j,^ = pu^Uy with p = M 5{x) the density of matter 
time-independent we obtain the socalled Schwarzschild solution in standard coordinates: 



dr""^ 
1-^ 



(2.8) 



where = 2GM is the so called Schwarzschild radius. 



Metric (|2.8I) determines completely the gravitational field in the vacuum generated by a spher- 
ically matter density distribution. Furthermore the Schwarzschild solution is valid also when we 
consider a moving source with a spherical distribution. The spatial metric is determined by expres- 
sion of spatial distance element 
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dr""^ 

dl'' = -3^ + r"Mn . (2.9) 

We have to note that, while the length of circumference with "radius" r" is the usual one 2nr" , the 
distance between two points on the same radius is given by the integral 




(2.10) 



then the space is curved. Besides we note that gu < 1, then, by the relation between the time 
coordinate t and the proper time r (dr = ^/g^dt), we get the condition 



dT<dt. (2.11) 

At infinity, the time coordinate coincides with physical time. We can state that when we are at a 
finite distance from the the mass, there is a slowdown of the time with respect to the time measured 
at infinity. 

In presence of matter the situation is the following. In fact from the first equation in the (12.71) . 
when r — 7- 0, r) has to vanish as r^; otherwise T*^ could have a singular point in the origin. 
By integrating formally the equation with the condition r)|r = o = 0, one get 



fi{t,r) 



In 



r 



T\ dr 



(2.12) 



It is easy to demonstrate also in the matter with spherical symmetry that the proprieties (|2.10l) . 
(12.1 II) and n(t, r) + u{t, r) < are verified II162L If the gravitational field is created by spherical 
body with "radius" ^, we have T*^ = outside the body (r > ^) and we can write 



^ii^^r) = -In 



1 / TKPdf 



(2.13) 



and obtain the analogous expression of (12.81 ) in the matter: 
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^ _ rg{r") 



dr"'^ 

1 _ fgi'^ ) 

-L // 



(2.14) 



where we introduced the Schwarzschild radius linked to the quantity of matter included in the 
sphere with radius r": 



rg{r") = X T\r^df-, 
Jo 



(2.15) 



obviously when the distance is bigger than the radius of the body, the metric (12.141) is equal to 
If we consider the transformation (11.901 ). which in the case of Schwarzschild solution is 



r"^-rgr" 



(2.16) 



it is possible to obtain the Schwarzschild solution (|2.8I) in isotropic coordinates: 



ds' 



'1 - 


^9 1 

At' 


2 

dt^- 




.1 + 


4r' J 




4r'_ 



{dr" + r"dn) 



(2.17) 



Li both cases, the solutions (12.81) and (12.141) agree with the trace equation of Einstein equation: 
R = —XT. Since in the vacuum the trace of matter tensor is vanishing (except the origin, in 
which the trace is proportional to 5(x)) we can state that the Schwarzschild solution is "Ricci flat": 
R = 0. 



If we add in the Hilbert - Einstein lagrangian (11.131) a term like {—2^/^ A) with A a generic 
constant the field equations (11.101) are modified as follows 



and if we consider a point-like source, we find the Schwarzschild - de Sitter solution 



(2.18) 



36 



Chapter 2 Exact and perturbative solutions in General Relativity 



— H r 



//2 



dr 



112 



- r"^dn 



1 _ Is _|_ ^^-"2 

J- -r 3' 



(2.19) 



In this case the trace of (12.181) is 



R = AA-XT (2.20) 

from which we note that this solution does not admit solution in the vacuum, since also in absence 
of ordinary matter {T^„ = 0) we have a nonvanishing scalar curvature. The contribution is given 
by cosmological constant A. It is also possible in this case to find the analogous of (|2.14l) . 

Finally let us consider as source a radial and static electric field E = Qx /\xf. We know that 
the Lagrangian of electromagnetic field is —-^FapF'^l^ where F^p is the electromagnetic tensor. 
Then, the Hilbert - Einstein lagrangian is 



Che = V^iR - — , (2.21) 

and the Einstein equation (|1.10l) becomes 



= {g,uFo.0F^^ - 4F^„F"J . (2.22) 
The solution for a spherically symmetric system is the Reissner - Nordstrom solution: 



ds' 



1-^ + 



de- 



dr 



112 



J12 



19. I 



01 

r."2 



dQ. 



(2.23) 



In all the above cases shown the Birkhoff theorem holds: The metric tensor generated in vac- 
uum by a matter density distribution with a spherical symmetry is time-independent. Also a time- 
dependent source with a spherical symmetry produces a static metric. The curvature of spacetime 
in the matter, a distance r from the origin, is proportional only to the matter inside the sphere of 
radius r. This conclusion is compatible with the Gauss theorem of classical mechanics. 

One of the goal of the present thesis is to develop similar considerations in the case of /-gravity. 



2.2 Perturbations of the Schwarzschild solution: The Eddington parameters /3 and 7 
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2.2 Perturbations of the Schwarzschild solution: The Edding- 
ton parameters p and 7 

The Schwarzschild solution (12.171) is a mathematically exact solution and is true everywhere. But 
in same cases the physical conditions could permit a "reduction" of them. In fact the Schwarzschild 
radius r-g is a scale-length induced by theory. Then we could be at radial distance r' for the which 
we have rg/r' <^ 1 and the (|2.17l) becomes 



+ 



+ ... 



dt' 



1 + ^ + 



dr'^ + r'''dQ 



12, 



(2.24) 



Since we are interesting to investigate the deviations, induced by /-gravity, from behavior 
(12.241) it is useful to introduce the method taking into account such deviations with respect to GR. 
A standard approach is the Parameterized-Post-Newtonian (PPN) expansion of the Schwarzschild 
metric (12.171) . Eddington parameterized deviations with respect to GR, considering a Taylor series 
in term of Vg/r' assuming that in Solar System, the limit Vg/r' <^ 1 holds [I129II . The resulting 
metric is 



ds^ 



df 



dr" + r'^dVt 



(2.25) 



where a, (3 and 7 are unknown dimensionless parameters (Eddington parameters) which parame- 
terize deviations with respect to GR. The reason to carry out this expansion up to the order (r^/r')^ 
in gtt and only to the order {rg/r') in gij is that, in applications to celestial mechanics, gij always 
appears multiplied by an extra factor v'^ ^ M/r'). It is evident that the standard GR solution for 
a spherically symmetric gravitational system in vacuum, is obtained for a = /3 = 7 = 1 giving 
again the "perturbed" Schwarzschild solution (12.241) . Actually, the parameter a can be settled to 
the unity due to the mass definition of the system itself [I129II . As a consequence, the expanded 
metric (12.251) can be recast in the form : 



ds' 



dt' 



1 + 7^ 



dr'" - r"'dn . 



(2.26) 



where we have restored the standard spherical coordinates by means of the transformation r" 



38 



Chapter 2 Exact and perturbative solutions in General Relativity 



2 

. The two parameters (3, 7 have a physical interpretation. The parameter 7 measures 

the amount of curvature of space generated by a body of mass M at radius r'. In fact, the spatial 
components of the Riemann curvature tensor are, at post-Newtonian order. 



^ij^i = ipl^^^Jki (2.27) 

independently of the gauge choice, where Nijki represents the geometric tensor properties (e.g. 
symmetries of the Riemann tensor and so on). On the other side, the parameter /3 measures the 
amount of non-linearity (~ {fg/r'Y) in the gu component of the metric. However, this statement 
is valid only in the standard post-Newtonian gauge. 



2.3 General remarks on the Newtonian and the post — Newto- 
nian approximation of Einstein equation 

At this point, it is worth discussing some general issues on the Newtonian and post-Newtonian 
limits. Basically there are some general features one has to take into account when approaching 
these limits, whatever the underlying theory of gravitation is. In fact here we are not interested 
in entering the theoretical discussion on how to formulate a mathematically well founded New- 
tonian limit (and post-Newtonian) of general relativistic field theories, nevertheless we point the 
interested reader to ^ [T79l [TM [TSTl [T82l [T83l [T84ll . In this section, we provide the explicit 
form of the various quantities needed to compute the approximations in the field equations in GR 
theory and any metric theory of gravity. We only mention that there is also been a discussion on 
alternative ways to define the Newtonian and Post - Newtonian limit in higher-order theories in the 
recent literature, see for example [I185II . In this work, the Newtonian and Post - Newtonian limit is 
identified with the maximally symmetric solution, which is not necessarily Minkowski spacetime 
in / - theories which could be singular. 

If one consider a system of gravitationally interacting particles of mass M, the kinetic en- 
ergy iiVft;^ will be, roughly, of the same order of magnitude as the typical potential energy 
U = GM"^ /f, with M, f, and v the typical average values of masses, separations, and veloci- 
ties of these particles. As a consequence: 



y2 ^ 



GM 



(2.28) 



2.3 General remarks on the Newtonian and the post — Newtonian approximation of Einstein 

equation 39 



(for instance, a test particle in a circular orbit of radius r about a central mass M will have velocity 
V given in Newtonian mechanics by the exact formula f ^ = GM/r.) 

The post-Newtonian approximation can be described as a method for obtaining the motion of 
the system to higher approximations than the first order (approximation which coincides with the 
Newtonian mechanics) with respect to the quantities GM/f and v'^ assumed small with respect to 
the squared light speed. This approximation is sometimes referred to as an expansion in inverse 
powers of the light speed. 

The typical values of the Newtonian gravitational potential $ are nowhere larger (in modulus) 
than 10^^ in the Solar System (in geometrized units, $ is dimensionless). On the other hand, 
planetary velocities satisfy the condition v"^ < — $, while the matter pressure p experienced inside 
the Sun and the planets is generally smaller than the matter gravitational energy density — p$, in 
other words p/ p < — $. Furthermore one must consider that even other forms of energy in the 
Solar System (compressional energy, radiation, thermal energy, etc.) have small intensities and the 
specific energy density 11 (the ratio of the energy density to the rest-mass density) is related to U 
by n < f/ (IT is ~ 10"^ in the Sun and ~ 10"^ in the Earth [I129II ). As matter of fact, one can 
consider that these quantities, as function of the velocity, give second order contributions : 



Therefore, the velocity v gives 0(1) terms in the velocity expansions, f/^ is of order 0(4), Uv of 
0(3), UU is of 0(4), and so on. Considering these approximations, one has 



~ ~ p/p ~ n ~ 0(2) . 



(2.29) 




(2.30) 



and 



\d/dt\ 



0(1). 



(2.31) 



V 



Now, particles move along geodesies : 



Typical values of p/ p are ^ 10 in the Sun and ^ 10 in the Earth II129I . 
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ds^ '^^ ds ds 

which can be written in details as 



d'^x^' ^„ dx" dx'' ^ 

+ T't- - = 0, (2.32) 



d^x' 



■ dx^ 
— _r' _ 9r* 



dx"" dx" 
dt dt 



+ 



Tl + 2T 



tm 



dx"^ , dx"^ dx^ 

+ 2r 



dt 



dt dt 



dx^ 
'dt 



(2.33) 



In the Newtonian approximation, that is vanishingly small velocities and only first-order terms in 
the difference between g^i, and the Minkowski metric r^^j,, one obtains that the particle motion 
equations reduce to the standard result : 



(2 34) 

dt' - « - 2 dx^ ■ ^ ^ 

The quantity 1 — gu is of order GM /f, so that the Newtonian approximation gives —-ir to the 

dt''' 

order GM jr^, that is, to the order v^/r. As a consequence if we would like to search for the post- 
Newtonian approximation, we need to compute -7-^ to the order v^j r. Due to the Equivalence 

iXv 

Principle and the differentiability of spacetime manifold, we expect that it should be possible to 
find out a coordinate system in which the metric tensor is nearly equal to the Minkowski one 77^,^, 
the correction being expandable in powers of GM /f ~ i;^. In other words one has to consider the 
metric developed as follows : 



' gu{t, x) ~ 1 + g^^\t, x) + g'i\t, x) + 0(6) 

< gu{t,^)^g'i\t,^)+0{5) (2.35) 

^ ^,,(i,x)~-5,, +4'^(i,x)+0(4) 
where 5ij is the Kronecker delta, and for the controvariant form of g^,,, one has 



2.3 General remarks on the Newtonian and the post — Newtonian approximation of Einstein 

equation 41 



' g'\t, x) ~ 1 + x) + x) + 0(6) 

< c/** {t, x) ~ c/(3)ti (^^^ x) + 0(5) (2.36) 

^ ^^^■(t,x)~-5,,+^(2)i,(^,x) + 0(4) 

The inverse of the metric tensor (12.351) is defined by The relations among the higher than 

first order terms turn out to be 



' g^'^'\t,x) = -gS\t,x) 
g('^'\t,x)=gS\t,xf-g',f{t,x) 

(2.37) 

^ g('^^^(t,x) = -gff{t,x) 

In evaluating F^^ we must take into account that the scale of distance and time, in our systems, are 
respectively set by f and f/v, thus the space and time derivatives should be regarded as being of 
order 



d 1 d V 
dx^ f ' dt f 



— (2.38) 



Using the above approximations (12.351) . (12.361) and (12.371 ) we have, from the definition (11.31) . 
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r(2)' = 1 o(2) - o(2) - o(2) 

jk 2 I ^jk,i ^ij,k ^ik,j 
r(3)^ _ 1 _ (3) _ (2) 

^ tj ~ 2 \ ytj,i yit,j yij,t 



ti ~ 2 1 ^ yim.ytt,ni '^yu^t 



The Ricci tensor components (11.41) are 



r(3)* 

^3 



tt 2^**,* 



— \ \9tt!i 9tt^ 9tt\ 



r(2)* - lo(2) 



(2.39) 



d(2) _ 1 (2) 



(4) _ l.W + 0^2) 1 (2) (2) 1 (2) _ 1 (2) (2) 

- 29tt,rnm ^ 2 9 "^"^^"^9^,11 ' 29^^9tt,mn ' 29mm.,tt 49tt,m9tt,m' 



_1 (2) (2) _ (3) 
^9mm,n9tt,n 9tm,tm 



T^{3) ^ 1 (3) _ 1 (2) _ 1 (3) 1 (2) 

-^ti 29ti,mm 29im,mt 29mt,mi ' 29mm,ti 



d(2) ^ 1 (2) _ 1 (2) _ 1 (2) _ 1 (2) 1 (2) 
V. ij 29ij,mm 29im,mj 29jm,mi 29tt,ij ~^ 29mm,ij 



(2.40) 



and the Ricci scalar (11.51) is 



2.3 General remarks on the Newtonian and the post — Newtonian approximation of Einstein 

equation 43 



r d(2) _ d{2) _ p(2) _ ^(2) _ ^(2) 
-iL — -Tt^j -TLmm — Vft mm Unn. 



(2) 



9tt,mm 9nn,m.m + gmn, 



d(4) _ p{4) ^(2)d(2) ^(2) p(2) 



1 (4) 1 (2) (2) 1 (2) (2) , 1 (2) _ 1 (2) (2) 

29tt,mm ' 25'"in,m5'tt,n + 23^'^9tt,mn + 23mm,tt i9tt,m9tt,m' 



0^2) _ (3) _ 1 (2) (2) _ 1 (2) (2) _ (2) 
4^ymm,nytt,n 9tm,tm 29tt ytt,rnm 29'"^'^\ 9mn,ll 9ml,ln^ 



(2) _ (2) (2) 
9nl,lm 9tt,mn ' 9u,m 



(2.41) 



The Einstein tensor components (11.1 II) are 



r r^(2) _ p(2) _ 1 p(2) _ 1 „(2) , 1 



(2) 

2 ymn,mn 



v{3) _ p(3) _ 1 (3) 



1.(2) 



1.(3) , 1.(2) 



/ ^ti ti 29ti,mm 29im,mt 29mt,mi ' 29mm,ti 



^(2) ^ d(2) , D(2) ^ 1 (2) _ 1„(2) _ 1.(2) _ 1.(2) , 1 „(2) , 
ij ij '2 2'^ij,mm 2^i'rn,mj 2-'j'm,mi 2-'tt,ij 2^'rnva,ij' 



1^ 
2 



(2) _ (2) , (2) 

9tt,mm 9nn,mm + 9mn,mn 



By assuming the harmonic gauge 



(2.42) 



9'''^';^ = 



(2.43) 



it is possible to simplify the components of Ricci tensor (12.401) . In fact for /x = one has 



^The gauge transformation is /i^^ = ft,^^ — (^^ ,y — d^^f^ when we perform a coordinate transformation as x'^ = 
+ with 0(C^)<C 1. To obtain our gauge and the validity of field equation for both perturbation /i^j, and /ip^ the 
have satisfy the harmonic condition D^^ = 0. 
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23^^^^^^ ^ g^l - 2g^l^ + g^^t = , (2.44) 



and for 11 = i 



2g^^r^, ^ g'i^l + 2g^Xm ' 9^±, = • (2.45) 
Differentiating Eq. (l2.44|) with respect to t, and (12.451) and with respect to t, one obtains 



(2.46) 



9tt,tj '^9mt,jm ~^ 9mm,tj ~ ^ ) (2.47) 



9tt'!ti + '^9mi.tm 9mm,ti ~ ^ • (2.48) 



On the other side, combining Eg. (12. 471) and Eg. (12. 481) . we get 



o(2) = (2 49) 

Jmm,tt Jmt,tm Jmt,mi ' y-^-' J 



Finally, differentiating Eg. (12.451) with respect to , one has : 



9tt!ij + '^9mi^jm, 9^mm,ij ~ ^ (2.50) 

and redefining indexes as j — > i, i — )■ j since these are mute indexes, we get 



(2) (2) _ (2) ^ 2 51) 



2.3 General remarks on the Newtonian and the post — Newtonian approximation of Einstein 

equation 45 



Combining Eq. (|2.50l) and Eq. (l2.51l) . we obtain 



0(2) , (2) (2) _ (2) ^ 2 52) 



Relations (|2.46l) . (12.491) . (I2.52|) guarantee us to rewrite Eqs. (12.401) as 



-Kit 



(3): 

ti \hg 



5(2)1 



d(4)| _ 1 An(^) -U in(2)^(2) _ i ^(2) _ i I ^ ^(2) 1 2 



(2.53) 



1 A.(2) 



and Eqs. (12.411) becomes 



p(2)| _ 1 An(2) IaJ"^) 

(2.54) 

d(4)| _ 1 A„(4) , 1 (2) (2) 1 (2) l| (2)|2 1 (2) a (2) i (2) . (2) 

^ '\HG — 2^9tt + 29^'^9tt,mn ~ 29tt,tt ~ 2 I V fl'ti I ~ 29tt ^9tt ~ 29mnl^gmn 

where V and A are, respectively, the gradient and the Laplacian in flat space. The Einstein tensor 
components (11.111) in the harmonic gauge are 



r^(2)| _ 1 Ar,(2) I 1 A«(2) 



^(4) 



tt \HG 



ti \HG 



^(2) I 

'^ij \hg 



(3) 



l^9l 



2 ^ij ~ 4 



(2.55) 



A^?;? - Ag., 



(2) 
mm 



On the matter side, i.e. right-hand side of the field equations (11.101) . we start with the general 
definition of the energy-momentum tensor of a perfect fluid 
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Tain = (p + Up + p) UaUp - PQap ■ 



(2.56) 



Following the procedure outlined in [I155II . we derive the explicit form of the energy-momentum 
as follows 



Ttt = P + p{v^ -2U + U)+p 



t;2 ( £ + t;2 + 2V + n I + a - 2UU 



I Tu = -pv' + p 



-v'i ^^ + 2V + + U] +ht 



(2.57) 



Tij = pv'v^ + p5ij + p 



We are now ready to make use of Einstein field equations (11.101) . which we assume in the form 



T, 



T 



(2.58) 



From their interpretation as the energy density, momentum density and momentum flux, then, we 
have Ttu Tu and Tij at various order 



r Tu = rf + ) + 0(4) 



{ Tu = T^p + 0(3) 



t T,, = Tjf + 0(4) 



(2.59) 



where T^i^^ denotes the term in T^,^ of order M/r^ . In particular T^^' is the density of rest 
mass, while t{P is the non-relativistic part of the energy density. What we need is 



.(0) 



5, 



fit/ 



T 



(2.60) 



2.3 General remarks on the Newtonian and the post — Newtonian approximation of Einstein 

equation 47 



But GM/f is of order v"^, so (l235l) and (|239|) give 

' Su = + SS^ + 0(6) 

< Su = SlP + 0(3) (2.61) 
where 5'^^'' denotes the term in S^^ of order M/f^ . In particular 



( 



(2) 



(1) 



(0) 



(0) 



It 



1^(2) , lrpi2) 



(2.62) 



(1) 



Using the (12.531) and (|2.61l) in the field equation (12.581) we find that the field equations in harmonic 
coordinates are indeed consistent with the expansions we are using, and give 



r?{2) y Q 



d(4) 



(0) 



(2) 



d(3) 



(0) 



r?{2) y c 



(0) 



(2.63) 



and in particular 



48 



Chapter 2 Exact and perturbative solutions in General Relativity 



AS 



rp(2) rp{2) 

J- II ~r -t mm 



(2) (2) (2) I I (2) 1 2 

9mngtt,mn + 9tt,tt + I V I 



X 5ij T^^ 



(0) 



(2.64) 



From the first one of (12.641) . we find, as expected, the Newtonian mechanics: 



Air X — X 



2G I dV ,*^ ^ = 2$(x) 
X — x' 



(2.65) 



where $(x) is the gravitational potential which, in the case of point-like source with mass M, is 



$fx) 



GM 



(2.66) 



From the third and fourth equations of (12.641) we find that 



ft?' = -4/<''x'^fc§l = Z.(x) 

(2.67) 



^ij 4tt J ^ |x-x'| 

The second equation of (12.641 ) can be rewritten as follows 



A 



■^tt ' -'- mm 



8$A$ + 2$4t 



(2.68) 



and the solution for g^-f^ is 



2.3 General remarks on the Newtonian and the post — Newtonian approximation of Einstein 

equation 49 



47r _/ |X — X'l TT 

/'rfV^^ = 2 0(x). (2.69) 
27r J |x — x'l 

By using the equations at second order we obtain the final expression for the correction at fourth 
order in the time-time component of the metric: 



-L^lls^^ (2.70) 

47r J |x — x'l 

We can rewrite the metric expression (|2.35|) as follows 



/ 1 + 2$ + 2 e 

^/^^ ~ ^ X n o^^ ^2.71) 

Finally the Lagrangian of a particle in presence of a gravitational field can be expressed as 
proportional to the invariant distance ds^l'^, thus we have : 



1/2 

1 + gi? + + 2^?;S^™ - + gl^lv^v- ] , (2.72) 



which, to the 0(2) order, reduces to the classic Newtonian Lagrangian of a test particle Lnbw = 

\ 1/2 

1 + 2$ — J , where = and |vp = v"^v^. As matter of fact, post-Newtonian physics 
has to involve higher than 0(2) order terms in the Lagrangian. In fact we obtain 



L 



2 



3 

+ 4 



(2.73) 



An important remark concerns the odd-order perturbation terms 0(1) or 0(3). Since, these 
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terms contain odd powers of velocity v or of time derivatives, they are related to the energy dis- 
sipation or absorption by the system. Nevertheless, the mass-energy conservation prevents the 
energy and mass losses and, as a consequence, prevents, in the Newtonian limit, terms of 0(1) 
and 0(3) orders in the Lagrangian. If one takes into account contributions higher than 0(4) order, 
different theories give different predictions. GR, for example, due to the conservation of post- 
Newtonian energy, forbids terms of 0(5) order; on the other hand, terms of 0(7) order can appear 
and are related to the energy lost by means of the gravitational radiation. 

2.4 General remarks on the post — Minkowskian approxima- 
tion of Einstein equation 

We suppose the metric to be close to the Minkowski metric r/^i,: 



g^iu = + h^,^ {2.1 A) 

with h^jj small quantities (0{h)'^ ^ 1). To first order in h, the Christoffel symbols (|1.3I) are 



r^. = 7;V'"'{h^.a,u + Ka,^. - Ku,.) . (2.75) 



As long as we restrict ourselves to first order in h, we must raise and lower all indices using rj^y, 
not g^y-, that is 



n'^^Kr = h\ = h, r/'^--_ = _, etc. (2.76) 

With this assumptions, the Ricci tensor and scalar (|1.4I) - (11.51) are then 



d(i) — h<r _ In /) — 1/) 



(2.77) 



where V°Va ~ y'' = is the d'Alembertian operator in the flat space. The field equation 
(11.101) becomes 



2.4 General remarks on the post — Minkowskian approximation of Einstein equation 
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Gi:}=I^-lR''\u = XT^ (2.78) 

where T^^ is fixed at zero-order in (|2.78l) since in this perturbation scheme the first order on 
Minkowski space has to be connected with the zero order of the standard matter energy momentum 
tensoiB Eqs. (|2.78l) in terms of /i^^, are 



hl,u). - ^0,h,, - -h,, - -[K^'^^ - n,h]r^,^ = XT^J . (2.79) 

Since T^^, is taken to the lowest order in h^^,, so it is independent of h^,,, it has to satisfies the 
ordinary conservation conditions: 

T"!"^ = 0. (2.80) 

Note that it is this form of the conservation law that is needed for the consistency of (12.791 ). because 
(11:801) implies 

G^'^'\^ = (2.81) 
whereas the linearized Ricci tensor satisfies Bianchi identities (|1.6I) of the form 



2 



(2.82) 



By choosing the transformation h^,y = h^^ — ^r]^^ and the gauge condition /i^*"^ = (harmonic 
gauge (|2.43|) ) one obtains that field equations read 



□L, = -2A'T(°). (2.83) 



■'in this perturbation scheme the first order on Minkowski space has to be connected with the zero order of the stan- 
dard matter energy momentum tensor. This formahsm descends from the theoretical setting of Newtonian mechanics 
which requires the appropriate scheme of approximation and coincides with a gravity theory analyzed at the first order 
of perturbations in the curved spacetime metric. 
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One solution is the retarded potential 



h,^it, x)=AG I d'x' ^^'^ (2.84) 



or in terms of perturbation h 





.-x'l) 




X — x' 


r j 



x)=AG j ^ . (2.85) 

The propagation of h^^ is possible with a particle massless. 



Chapter 3 

Spherical symmetry in / — gravity 



Spherical symmetry in /-gravity is discussed in details considering also the relations with the weak 
field limit [D]. Exact solutions are obtained for constant Ricci curvature scalar and for Ricci scalar 
depending on the radial coordinate. In particular, we discuss how to obtain results which can be 
consistently compared with General Relativity giving the well known post-Newtonian and post- 
Minkowskian limits. Furthermore, we implement a perturbation approach to obtain solutions up 
to the first order starting from spherically symmetric backgrounds. Exact solutions are given for 
several classes of /-theories in both = constant and i? = R{r). 



3.1 The Ricci curvature scalar in spherical symmetry 

Starting by the definition of Ricci scalar (|1.5I) and imposing the spherical symmetry (11.891) . the 
Ricci scalar in terms of the gravitational potentials {gu and grr) reads : 



R{t, r) 



2n2 

rr 



gttgr 



i2 

9tt 



r + gtt 



/ / \ 

9tt9rr) 



II 

r9tt 



rgr 



-^9l 



2 _ _|_ / 

9rr 9rr i fg,fj 



(3.1) 



where, for the sake of brevity, we have discarded the explicit dependence in gtt{ti r) and grr{t, r) 
and the prime indicates the derivative with respect to r while the dot with respect to t. If the metric 
(11.891) is time-independent, i.e. gtt{t,r) = a{r), grrit,r) = hir). the (13.11) assumes the simpler 
form : 
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R{r) 



a{r) 



2r^a(r)^6(r)^ 
-4a(r)^ ( b{r)^ - b{r) + rb'{r) 



2b{r) ( 2a'(r) + ra"(r) - ra'{r)b'{r) 



— b{r)a' {r^r'^ 



(3.2) 



where the radial dependence of the gravitational potentials is now explicitly shown. This expres- 
sion can be seen as a constraint for the functions a(r) and 6(r) once a specific form of Ricci scalar 
is given. In particular, it reduces to a Bernoulli equation of index two with respect to the metric 
potential 6(r) : 



b'{r) + 
2a(r 



r^a\rY - 4a(r)^ - 2ra(r)[2a(r)' + ra{r)"] 



ra{r)[4:a{r) + ra'(r)] 



b{r) 



+ 



2 + r'^R{r) 



r 4a(r) + ra'{r) 
A general solution of (13.31) is: 



6(r)2 = b'{r) + h{r)b{r) + l{r)b{rf 



(3.3) 



6(r) 



exp[— J dr h{r)] 



K+ J drl{r) exp[— J dr h{r)] 



(3.4) 



where K is an integration constant while h{r) and /(r) are the two functions which, respectively, 
define the coefficients of the quadratic and the linear term with respect to 6(r), as in the standard 

definition of the Bernoulli equation [I156II . Looking at the equation, we can notice that it is possible 

2 

to have /(r) = which implies to find out solutions with a Ricci scalar scaling as — ^ in term of 
the radial coordinate. On the other side, it is not possible to have h{r) = since otherwise we will 
get imaginary solutions. A particular consideration deserves the limit r — )• oo. In order to achieve 
a gravitational potential b(r) with the correct Minkowski limit, both h{r) and Z(r) have to go to 
zero provided that the quantity r'^R{r) turns out to be constant: this result implies b'{r) = 0, and, 
finally, also the metric potential b{r) has a correct Minkowski limit. In general, if we ask for the 
asymptotic flatness of the metric as a feature of the theory, the Ricci scalar has to evolve at infinity 
as r~" with n ^ 2. Formally it has to be: 



lim r^i?(r) 



(3.5) 



3.2 Solutions with constant curvature scalar 
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with n E N. Any other behavior of the Ricci scalar could affect the requirement of the correct 
asymptotic flatness. This result can be easily deduced from (13.31) . In fact, let us consider the 
simplest spherically symmetric case: 



a{r)dt^ - - r'^dVL . 
a[r) 



(3.6) 



The Bernoulli equation (13.31 ) is easy to integrate and the most general metric potential a(r), com- 
patible with the Ricci scalar constraint (13.21) . is : 



a(r) 



1 



ki ko 
— + — H 



r'^R{r)dr 



dr 



(3.7) 



where ki and k2 are integration constants. Actually one gets the standard result a(r) = 1 (Minkowski) 
for r — )■ oo only if the condition (13. 5|) is satisfied, otherwise we get a diverging gravitational po- 
tential. 



3.2 Solutions with constant curvature scalar 



Let us assume a scalar curvature constant (R = Rq). The field equations (11.251) and (|1.26l) reduce 
to: 



/'(i?o)i?o - 2/(i?o) = XT 



(3.8) 



Such equations can be arranged as: 



(3.9) 



Ro = qXT-A\ 



where A 



f{Ro) 
V'{Ro) 



and q ^ = f'{Ro). Since we are analyzing the weak field limit of HOG, it 
is reasonable to consider Lagrangians which work as the Hilbert - Einstein one when i? — )■ (this 
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even means that we can suitably put the cosmological constant to zero), that is: 



lim/~i?. (3.10) 

In such a case, the trace equation of (13.91 ) indicates that in the vacuum case (T^j, = 0), one obtains 
a solutions with constant curvature R = Rq. 

Let us now suppose that the above theory, for small curvature values, evolves to a constant as 
limj^_>o f = fo- Even in this case, considering only the trace equation of (13.81) . some interesting 
features emerge. In fact if we consider the expression for / : 

/ = fo + fiR + HR). (3.11) 

where fo and fi are a coupling constants, while -FiR) is a generic analytic function of R satisfying 
the condition 

\imR~'^J'(R) = constant, (3.12) 
it is evident that no zero - curvature solutions are obtainable since : 



r(i?o)i?o-2^(i?o)-/ii?o-2/o = XT. (3.13) 

Furthermore, in this case, even in absence of matter, there are no Ricci flat solution of the field 
equations since the higher order derivative terms give curvature constant solutions corresponding 
to a sort of effective cosmological constant. Of course, this is not the case for GR. In fact in 
the standard Einstein theory, constant curvature solutions different from zero are in order only 
in presence of matter because of the proportionality of the Ricci scalar to the trace of energy - 
momentum tensor of matter. Besides, one can get a similar situation in presence of a cosmological 
constant put by hands into the dynamics. 

In other words, the difference between GR and HOG is that the Schwarzschild - de Sitter solu- 
tion is not necessarily given by a A - term while the effect of an "effective" cosmological constant 
can be achieved by the higher order derivative contributions. This result has been extensively 
investigated in several recent papers as, for example, H 14011 . For a discussion, see also p4l. 



3.2 Solutions with constant curvature scalar 
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Let us consider now the search for a general solution of (11.251) and (11.261) considering a spher- 
ically symmetric metric as (11.891) . By substituting the metric into the field equations. One obtains 
that the {t, r} - component of (11.251) gives ^g^'^^^^f^-^ = 0. This means that grr{t, r) has to be time 
independent and we can write grri^, r) = h{r). On the other hand, by the {6*, 6} - component of 
(|1.25l) . one gets the relation = ^(r) with ^(r) a given time independent function : 



gtt{t,r) = a(t)exp 



2 / dr 



[1 - r\X + qA:p)]b{7 



(3.14) 



where A and q are defined as above and p is the pressure of a perfect fluid being the stress-energy 
tensor of matter 



Tf,u = {p + p)ui,Uu-pgf,u] (3.15) 

p is the energy density and = dx^/ds is the 4- velocity. Therefore, the function guit, r) has to 
be a separable function of time and radial coordinates, i.e. gtt{t, r) = d{t)a{r). As a matter of 
facts, the metric (11.891) becomes 



ds"^ = d{t)a{r)dt^ - b{r)dr^ - r^dn (3.16) 

which can be recast as 



ds'^ = a{r)dP - b{r)dr^ - r^dn , (3.17) 

by a suitable time - redefinition di = ^/d{t)dt. Nevertheless we redefine the time i as t. 

This exact result states that whenever one has a constant scalar curvature spacetime, any spher- 
ically symmetric background has to be necessarily static. In other words, this means that the 
Birkhoff theorem holds for the / - theories with constant curvature as it has to be (see II157II ). It 
has to be noticed that such a result is in striking contrast with what has been argued elsewhere 
about the fact that Birkhoff theorem does not hold, in general, for HOG theories. 

A remark is in order at this point. We have obtained this result considering a constant Ricci 
scalar spacetime and deducing some conditions on the form of gravitational potentials. Never- 
theless one can even reverse the problem. It is possible to argue that whenever the gravitational 
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potential gtt{t, r) is described by a separable functions and grr{t, r) is time - independent, by the 
definition of the Ricci scalar, one gets that R = Rq and at the same time the final solutions of the 
field equations will be static if the spherical symmetry is invoked. 

Actually, for a complete analysis of the problem, one should take into account the remaining 
field equations descending from (11.251) and (|1.26l) which have to be solved by taking even into 
account the definition of the Ricci scalar (|3.2I) . In summary, we have to solve the system : 



f Rtt + Xa{r)-qX[p + p{l-a{r))] = 
Rrr - Xb{r) - qXpb{r) = 
Ro - qX{p - 3j9) + 4A = 

^ R{a{r),h{r)) = Ro 
A general solution of the above set of equations is achieved for p = —p and reads 



(3.18) 



ds^= l+'^ + '-^P^Adt^ ^^V^^ r^dn. (3.19) 

\ / r 3 

In other words, any / - theory in the case of constant curvature scalar (R = Rq) exhibits solutions 
with cosmological constant as the solution Schwarzschild - de Sitter (12.191) if we consider the 
relation A = qXp — A = ^^jt^"^^- This is one of the reason why the dark energy issue can be 
addressed using these theory [|2S[lk|57l[51[59l[60l|6B[6l[63l. This fact is well known using the 
FRW metric [gl. 

If we neglect the cosmological constant /q and fi is set to zero in the (|3.1 II) . we obtain a 
new class of theories which, in the limit i? — )■ 0, does not reproduce GR (from (13.121) . we have 
lim^^o / ~ R^)- In such a case analyzing the whole set of (11.251 ) and (11.261 ). one can observe that 
both zero and constant 7^ curvature solutions are possible. In particular, if i? = _Ro = the field 
equations are solved for every form of gravitational potentials entering the spherically symmetric 
background (|3.17l) . provided that the Bernoulli equation (|3.3I) . relating these functions, is fulfilled 
for the particular case R{r) = 0. The solutions are thus defined by the relation 
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/ - theory: 
R 



R 



ii+R 



Field equations: 



R^^ = 0, with R = 



Rfiu — 



= 



with R = 0, ^1^0 



with R 



1 

n-l 



("-2)6 

with R = 0, ^1 = 



R^lv + ^g^iu = with = i?o, ^1 = 0, n = 2 
RfMu + ^Qixv = with R = 



(m+2)g2 



1 

m + l 



RtMu + ^9fj.u = 0, with R = Rq so that 
eii?r' + (2 - n)6i?o"'"" + ("^ + 2)^3 = 



Rfiu = 



with R = 



RtJiy + Asf^!. = with i? = 
^/.i. + A^^^ = 0, with = 



Table 3.1: Examples of / - models admitting constant and zero scalar curvature solutions. In the right hand 
side, the field equations are given for each model. The power n, m are natural numbers while are generic 
real constants. 



exp[— J dr h{r)\ 



J r[a{r)+ra' (r)] 



In table we give some examples of / - theories admitting solutions with constant 7^ or null 
scalar curvature. Each model admits Schwarzschild, Schwarzschild - de Sitter, and the class of 
solutions given by (13.201 ). 
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3.3 Solutions with curvature scalar function of r 

Up to now we have discussed the behavior of /-gravity seeking for spherically symmetric solutions 
(11.891 ) with constant scalar curvature. This situation is well known in GR and give rise to the 
Schwarzschild solution (i? = 0) and the Schwarzschild - de Sitter solution (i? = Rq ^ 0). The 
problem can be generalized in / - gravity investigating considering the Ricci scalar as an arbitrary 
function of the radial coordinate r. 

This approach is interesting since, in general, HOG theories are supposed to admit such kind 
of solutions and several examples have been found in literature [[TOl |32l |119[ |139[|140L Here we 
want to face the problem from general point of view. 

If we choose the Ricci scalar i? as a generic function of the radial coordinate {R = R{r)), it is 
possible to show that also in this case the solution of the field equations (|1.25l) and (11.261) is time 
independent (if T^,y = 0). In other words, the Birkhoff theorem has to hold. The crucial point of 
the approach is to study the off - diagonal {t, r} - component of (11.251) as well as in the case of GR. 
This equation, for a generic / reads : 



^[^'f')9rr{t,r) = 0, (3.21) 

and two possibilities are in order. Firstly, we can choose gtt{t, r) ^ 0. This choice implies that 
/' ~ r^^. If this is the case, the remaining field equation turn out to be not fulfilled and it can be 
easily recognized that the dynamical system encounters a mathematical incompatibility. 

The only possible solution is given by grr{t, r) = and then the gravitational potential has 
to be grrit, r) = h{r). Considering also the {6, 6} - equation of (11.251) one can determine that the 
gravitational potential gtt{t, r) can be factorized with respect to the time, so that we get solutions 
of the type (13.161) which can be recast in the stationary spherically symmetric form (13.171) after a 
suitable coordinate transformation. 

As a matter of fact, even the more general radial dependent case admit time - independent 
solutions. From the trace equation and the {0,9}- component, we deduce a relation which links 
a(r) and b{r) : 



2 r (fl/-2/)h(r) ^ 

«(0 = ^^-^1^2 ' (3-22) 

(with /" 7^ 0) and one which relates b{r) and / (see also PI 3911 1401 for a similar result) : 
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b{r) 



6[f'{rR'f"y - rR'^f"^] 



(3.23) 



rfirR'f" - 4/') + 2f'{rR{f' - rR'f") - 3R'f") ' 



As above, three further equations has to be satisfied to completely solve the system (respec- 
tively the {t, t} and {r, r} components of the field equations and the Ricci scalar constraint) while 
the only unknown functions are / and the Ricci scalar R{r). 

If we now consider a fourth order model of the form f = R + F{R), with F{R) ^ i? we are 
capable of satisfying the whole set of equations up to third order in F. In particular, we can solve 
the whole set of equations: the relations (13.221) and (|3.23l) will give the general solution depending 
only on the forms of F and R = R{r), that is: 



Once the radial dependence of the scalar curvature is obtained, (13.241) allow to write down the 
solution of the field equations and the gravitational potential, related to the function a(r), can 
be deduced. Furthermore one can check the physical relevance of such a potential by means of 
astrophysical data, see for example the analysis in [I152II . 

3.4 Perturbing the spherically symmetric solutions 

The search for solutions in / - gravity, in the case of Ricci scalar dependent on the radial coordinate, 
can be faced by means of a perturbation approach. There are several perturbation techniques 
by which HOG can be investigated in the weak field limit. A general approach is starting from 
analytical / - theories assuming that the background model slightly deviates from the Einstein GR 
(this means to consider f = R + J^{R) where J^{R) <^ R as above). Another approach can be 
developed starting from the background metric considered as the 0th - order solution. Both these 
approaches assume the weak field limit of a given HOG theory as a correction to GR, supposing 
that zero order approximation should yield the standard lore. 

Both these methods can provide interesting results on the astrophysical scales where spheri- 
cally symmetric solutions characterized by small values of the scalar curvature, can be taken into 
account. 




(3.24) 
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In the following, we will consider the first approach assuming that the background metric 
matches, at zero order, the GR solutions. 

In general, searching for solutions by a perturbation technique means to perturb the metric 



(3.25) 



This implies that the field equations (11.251) and (11.261) split, up to first order, in two levels. Besides, 
a perturbation on the metric acts also on the Ricci scalar R (11.51) : 



(3.26) 



and then we can Taylor expand the analytic / about the background value of R, i.e.: 



R - 



R - i?W 



df 



= En 



/"+3(fl(0)) 



R - 



R - 



+fm + f"mRW + f!X^Rii)^ 



f///(o) I f/y(o) d(i) _ df" 



df 
dRW 



(3.27) 



However the above condition J-'{R) <^ R has to imply the validity of the linear approximation 
^ ' ^ ^' "^^^^ demonstrated by assuming f = 1 + and /" = J-'". Immediately we 

obtain that the condition is fulfilled for 
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For example, given a Lagrangian of the form / = R+ the (|3.28l) means 



while, for / = i? + ai?^ the dim]) means 



^^'T « 1 , (3-29) 



< 1 . (3.30) 



1 + 2ai?(o) 



This means that the validity of the approximation strictly depends on the form of the models and 
the value of the parameters, in the previous case Rq and a. For the considerations below, we will 
assume that it holds. A detailed discussion for the Palatini formalism is in [|142[ 11861 . 

The zero order field equations read : 

- IgfJf^'' + = ^ T^'J (3.31) 



where 



In V^,,R^A ] -r^'^U^jR^J - g^^Jg^'^'-RfR^'J \ • (3.32) 



At first order one has : 



/'^°^{4l^-^^7i°Ji?^^^}+r^°^^^^^4^^^ = (3.33) 



with 
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-17(1) _ _/-//(0) I _ -pW^ /?(!)_ f(I)'' C'(0)_„(0) (0)pa d(1) , „(l)p^T d(0) , 



P ,0- 



T^Z^Rf - 9t} [9^'^'\,R^? + g^'^'^'R'il + 9^'^''' In v^!>SM 1+ 



_^/W0)^(i) 1^(0)^(0) _ ,(0),(0W^(0)^(0) 



(3.34) 



A part the analyticity, no hypothesis has been invoked on the form of /. As a matter of fact, / can 
be completely general. At this level, to solve the problem, it is required the zero order solution of 
(13.311) which, in general, could not be a GR solution. This problem can be overcome assuming the 
same order of perturbation on the /, that is: 



/ = r + :f{r) 



(3.35) 



where is a generical function of the Ricci scalar as above. Then we have 



(3.36) 



and the (13.311) reduce to the form 



= = A'TW. (3.37) 
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On the other hand, the (|3.33l) reduce to 

R^} - la'^jR''' - la^jR''' - ^a^^F^' + ^"''R^l + = ^ ^'J (3-38) 

where 

-9^^^ [g^'^'^^.R^' + g^'^^^R^l + 9^''"^ in v^!>!°)) } • (3.39) 

The new system of field equations is evidently simpler than the starting one and once the zero order 
solution is obtained, the solutions at the first order correction can be easily achieved. In Tables [X2l 
and l3.3[ a list of solutions, obtained with this perturbation method, is given considering different 
classes of / - models. 

Some remarks on these solutions are in order at this point. In the case of / models which 
are evidently corrections to the Hilbert - Einstein Lagrangian as A + _R + e-Rlni? and R + eR^, 
with e ^ 1, one obtains exact solutions for the gravitational potentials a(r) and h{r) related by 
a(r) = 6(r)~^. The first order expansion is straightforward as in GR. If the functions a(r) and 
6(r) are not related, for / = A + i? + ei? In i?, the first order system is directly solved without any 
prescription on the perturbation functions x(r) and y{r). This is not the case for f = R + ei?" 
since, for this model, one can obtain an explicit constraint on the perturbation function implying 
the possibility to deduce the form of the gravitational potential $(x) from a(r) = 1 + 2 $(x). 
In such a case, no corrections are found with respect to the standard Newtonian potential. The 
theories / = i?" and / = ^^^^^^^-^ show similar behaviors. The case f = R^ is peculiar and it 
has to be dealt independently. 
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/ - theory: 
spherical potentials: 
solutions: 
first order metric: 



solutions: 



A + R + eRlnR 

a{r) = b{r)-^ = 1 + ^ - ik! + Sx{r) 



x(r) = ^ + ^MM-2A)-iy 



a{r) = 1 - ^ + 5x{r), b{r) = + Sy{r) 



x{r) = {Ar^-6)\h + J 



dr 



36r(5(Ar2-6) 

+18)y{r) + r{36reA[log(-2A) - 1] 
+5(Ar2-6)V(^)} 



fc2^-6r^eA[ln(-2A)^l] 
r<5(r-2A-6)2 



45(2AV^ - 15Ar2 



/ - theory: 
spherical potentials: 
solutions: 
first order metric: 
solutions: 



air I 



6(r)-^ = 1 + ^ + Sx{r) 
x{r) = ^ 
a{r) = l + 6^, b{r) = l + 6^ 
x{r) = ki + k2r, y{r) = k^ 



f - theory: 



first order metric: 



solutions: 



R/{Ro + R) 
a{r) = l + S^, b{r) = l + 6^ 



'-^^k, - ^%7^fc2 + k,r 



x{r) = 

VV) — 5h372 '^1 ^372 1^2 



r 



,1/2, 



Table 3.2: A list of exact solutions obtained via the perturbation approach for several classes of / - theories; 
ki are integration constants; the potentials a(r) and h{r) are defined by the metric (I3.17I ). 
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/ - theory: 
spherical potentials: 



a(r) = 6(r) 



-1 _ 1 _|_ fei _|_ Rpr'^ 
r 12 



+ 5x{r) 



solutions: 



first order metric: 



solutions: 



first order metric: 



solutions: 



'n = 2, Ro^Q and x{r) = + ^ 



exp 



flQrp In(r-ro) 



with ro satisfying the condition 

6ki + 8ro + Ro4 = 

n > 2, System solved only whit Rq = 
and no prescriptions on a;(r) 



n = 2 y{r) = -^-^ + lrx'{r) + h, 
R{r) = 6Ro 

n^2 y{r) = -\ j dr r'^R{r) - ^ + \rx\r) + ki 
with R{r) whatever 

a(r) = 1 — y + 5x{r), h{r) = j^^- + 6y{r) 



n = 2 y{r) 



rk-i 



+ 



(3rg-4r)(r9-r)2 



77. 7^ 2 whatever functions x(r) , ?/(r) and R{r) 



Table 3.3: A list of exact solutions obtained via the perturbation approach for several classes of / - theories; 
ki are integration constants; the potentials a(r) and b{r) are defined by the metric (I3.17I ). 



Chapter 4 

The Noether Symmetries of / — gravity 



We search for spherically symmetric solutions of /-theories of gravity via the Noether Symmetry 
Approach [B]. A general formalism in the metric framework is developed considering a point-like 
/ - Lagrangian where spherical symmetry is required. New exact solutions are given. 

4.1 The point-like / Lagrangian in spherical symmetry 

As hinted in the plan of thesis, the aim of this chapter is to work out an approach to obtain spher- 
ically symmetric solutions in HOG by means of Noether Symmetries. In order to develop this 
approach, we need to deduce a point-like Lagrangian from the action (11.241) . Such a Lagrangian 
can be obtained by imposing the spherical symmetry in the field action (11.241) . As a consequence, 
the infinite number of degrees of freedom of the original field theory will be reduced to a finite 
number. The technique is based on the choice of a suitable Lagrange multiplier defined by assum- 
ing the Ricci scalar, argument of the function / in spherical symmetry. Elsewhere, this approach 
has been successfully used for the FRW metric with the purpose to find out cosmological solutions 
ll26l [T74l[T87l[T88]| . 

In general, a spherically symmetric spacetime can be described assuming that the metric (11.861 ) 
is time independent : 

ds"^ = A{r)df - B{r)dr^ - M{r)dn , (4.1) 

where gtt{t' , r) = A{r), grr{t' ,r) = B{r) and gQQ{t' ,r) = M(r). Obviously the conditions 
M(r) = and B{r) = A^^{r) are requested to obtain the standard Schwarzschild case of GR. 
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Our goal is to reduce the field action (11.241) to a form with a finite degrees of freedom, that is the 
canonical action of the form 



A= j drC{A,A',B,B',M, M',R,R') 



(4.2) 



where the Ricci scalar R and the potentials A, B, M are the set of independent variables defining 
the configuration space. Prime indicates now the derivative with respect to the radial coordinate r. 
In order to achieve the point-like Lagrangian in this set of coordinates, we write, in the vacuum, 
the (fL24l) as 



A = / d^x^/^ 



f{R)-X{R-R) 



(4.3) 



where A is the Lagrangian multiplier and R is the Ricci scalar (11.51) expressed in terms of the metric 



R 



A" 



M" AM' 



A 



12 



M 



12 



AB' B'M' 



AB ' BM ' ABM 2A^B 
which can be recast in the more compact form 



25M2 2A52 B^M M 



(4.4) 



, A" M" 

R = R* + — + 2- 



AB BM 



(4.5) 



where R* collects first order derivative terms. The Lagrange multiplier A is obtained by varying 
the action (|4.3I) with respect to R. One gets A = /J3. By expressing the determinant g and R in 
terms of A, B and M, we have, from (14.31 ), 



A = / drA'/^B'^Hl 



A" M" \ 
f-f„{R-R-- — -2 — 



drlA^I^B^'^M 



f-fniR-R*) 



The two lines differs for a divergence term which we discard integrating by parts. Therefore, the 



'in this chapter the derivatives of / will be indicated like = //j. 
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point-like Lagrangian becomes : 



2M51/2 

-i?'M' - A^''^B^'\2 + Mi?)/ij - Mf] 



51/2 



(4.7) 



which is canonical since only the configuration variables and their first order derivatives with re- 
spect to r are present. Eq. (|4.7I) can be recast in a more compact form introducing the matrix 
formalism : 



C = q''fq' + V 



(4.8) 



where q = {A, B, M, R) and q' = (A' , B' , M' , R') are the generalized positions and velocities 
associated to C The index indicates the transposed column vector. The kinetic tensor is given 

by Tij = . V = V(q) is the potential depending only on the configuration variables. The 

dq'idq'j 

Euler - Lagrange equations read 



dr 



dr 



fq']-V,V-q"iV,f]q' 



2fq" + 2 g' ■ V,r g' - V,V - q"i V,f g' = 



(4.9) 



which furnish the equations of motion in term of A, B, M and R, respectively. The field equation 
for R corresponds to the constraint among the configuration coordinates. It is worth noting that 

is zero. This result clearly depends on the absence 



the Hessian determinant of (14.71) . „ , „ , 

dq-dq'j 

of the generalized velocity B' into the point - like Lagrangian. As matter of fact, using a point-like 
Lagrangian approach implies that the metric variable B does not contributes to dynamics, but the 
equation of motion for B has to be considered as a further constraint equation. 

Beside the Euler - Lagrange equations (14.91 ). one has to take into account the energy Ec '■ 



Er 



q' ■ Vg'C - C 



(4.10) 
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which can be easily recognized to be coincident with the Euler-Lagrangian equation for the com- 
ponent B of the generalized position q. Then the Lagrangian (14.71) has three degrees of freedom 
and not four, as we would expect "a priori". 

Now, since the motion equation describing the evolution of the metric potential B does not 
depends on its derivative, it can be explicitly solved in term of i? as a function of other coordinates : 



2AM[{2 + MR)fn-Mf] ' ^ ' 

By inserting the (14.1 II ) into the Lagrangian (14.71) . we obtain a non-vanishing Hessian matrix re- 
moving the singular dynamics. The new Lagrangian readj^ 



(4.12) 



with 



x[2M^fRnA'R' + 2MM\fRA' + 2AfnRR') + AJrM'^] . (4.13) 

Since — — = 0, L is canonical (L is the quadratic form of generalized velocities. A' , M' and R' 

or 

and then coincides with the Hamiltonian), so that we can consider L as the new Lagrangian with 
three degrees of freedom. The crucial point of such a replacement is that the Hessian determinant 
is now non - vanishing, being : 



= 3AM[(2 + MR)fR - MfffRfRR^ . (4.14) 

Obviously, we are supposing that (2 + MR)/^ — Mf ^ 0, otherwise the above definitions of 
B, (14.111) . and L, (14.131) . lose of significance, besides we assume /rr 7^ to admit a wide class 
of HOG models. The case f = R requires a different investigation. In fact, considering the 

^Lowering the dimension of configuration space through the substitution ( I4.1 It does not affect the dynamics, since 
i? is a non-evolving quantity. In fact, introducing the ( I4.11I ) directly into the dynamical equations given by ( |4.7l l, they 
coincide with those derived by ( 14.131 ). 



dq[dq'. 
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GR point - like Lagrangian needs a further lowering of degrees of freedom of the system and the 
previous results cannot be straightforwardly considered. From (14.71) . we get : 



4 1/2 



^1/2^1/2 



A'M' - 2A'/^B'/^ , 



(4.15) 



whose Euler-Lagrange equations provide the standard equations of GR for Schwarzschild metric. 
It is easy to see the absence of the generalized velocity B' in (14.151) . Again, the Hessian determinant 
is zero. Nevertheless, considering, as above, the constraint (14.1 II) for B, it is possible to obtain a 
Lagrangian with a non- vanishing Hessian. In particular one has : 



5, 



GR 



M'2 A'M' 

+ 



AM 



2A 



(4.16) 



GR 



.1/2 
■'GR 



M'{2MA' + AM') 
M 



(4.17) 



and then the Hessian determinant is 



GR 



dq'idq'j 



-1, 



(4.18) 



which is nothing else but a non-vanishing sub-matrix of the / Hessian matrix. 

Considering the Euler-Lagrange equations coming from (|4.16l) and (14.171) . one obtains the 
vacuum solutions of GR (12.81) . that is : 



A = kA- 



r + ki 



B 



kok. 



2^4 



A 



M = k2{r + kif 



(4.19) 



In particular, the standard form of Schwarzschild solution is obtained for /ci = 0, A;2 = 1, /cs = Vg 
and ^4 = 1. 

A formal summary of the field equations descending from the point - like Lagrangian and their 
relation with respect to the ones of the standard approach is given in Tab. (|??1) . 

The explicit form of field equations (11.251) - (11.261) in the vacuum with the metric (14.11) are 
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Field equations approach 



6 J d^x 



= 

i 



H = g'^^Ha^ = 

; 

Hj.j. — 
Hge = csc^ 6lif^^ = 



Point-Uke Lagrangian approach 

; 

6jdrC = 
i 

- = 



i 

d dC dC_ n 

dr dB' dB ^ ^ 

d dC dC n 

dr dM' dM ~ ^ 

A combination of the above equations 



Table 4.1: The field-equations approach and the point-Uke Lagrangian approach differ since the symmetry, 
in our case the spherical one, can be imposed whether in the field equations, after standard variation with 
respect to the metric, or directly into the Lagrangian, which becomes point-like. The energy Ec corresponds 
to the tt - component of H^^y. The absence of B' in the Lagrangian implies the proportionality between the 
constraint equation for B and the energy function Ec- As a consequence, the number of independent 
equations is three (as the number of unknown functions). Finally it is obvious the correspondence between 
99 component and field equation for M. 



Htt = 2A^B^Mf + {BMA'^ - A[2BA'M' + M{2BA" - A'B')]}fR + 
+ {-2A^MB'R! + AA^BM'R! + AA^BMR")fRR + 

+4A^BMR'^fRRR^0, (4.20) 



Hrr = 2A^B'^M^f + {BM^A'^ + AM'^A'B' + 2A^MB'M' + 2A^BM'^ + 

-2ABNPA" - AA^BMM")fR + {2ABA'PAR! -t- 

+AA^BMM'E!)fRR = , (4.21) 



4.2 The Noether Symmetry Approach 75 



Hee = 2AB^Mf + {AAB^ ~ BA'M' + AB'M' -2ABM")fR + 

+ {2BMA'R' - 2AM B'R' + 2ABM'R' + AABMR")fRR + 
+4ABMR''fRnR = 0, (4.22) 



H^^ = sm^eHee = 0. (4.23) 



H = AAB^Mf - 2AB'^MRfR + 3{BMA'R' - AMB'R' + 

+2ABM'R' + 2ABMR!')fRR + 6ABMR''^fRRR = (4.24) 

4.2 The Noether Symmetry Approach 

In order to find out solutions for the Lagrangian (14.131) . we can search for symmetries related to 
cyclic variables and then reduce dynamics. This approach allows, in principle, to select /-gravity 
models compatible with spherical symmetry. As a general remark, the Noether Theorem states that 
conserved quantities are related to the existence of cyclic variables into dynamics [|167lll89lll90ll . 
Let us give a summary of the approach for finite dimensional dynamical systems. 
Let C{q\ cf) be a canonical, non-degenerate point- like Lagrangian where 



— = deti^i,- = det 

OA 



7^0, (4.25) 



with Hij as above, the Hessian matrix related to C. The dot indicates derivatives with respect to the 
affine parameter A which, ordinarily, corresponds to time t. In our case, it is the radial coordinate 
r. In standard problems of analytical mechanics, C is in the form 



£ = r(q,q)-V(q), (4.26) 
where T and V are the "kinetic" and "potential energy" respectively. T is a positive definite 
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quadratic form in q. The energy function associated with L is 



E,^—ci-L^ (4.27) 

which is the total energy T . It has to be noted that is, in any case, a constant of motion. 
In this formalism, we are going to consider only transformations which are point-transformations. 
Any invertible and smooth transformation of the "positions" = Q*(q) induces a transformation 
of the "velocities" such that 



QXq) = ^g^ (4.28) 

the matrix i7 = WdQ" jdq'W is the Jacobian of the transformation on the positions, and it is assumed 
to be nonzero. The Jacobian J of the "induced" transformation is easily derived and ^7 7^ — )■ 
^7 7^ 0. Usually, this condition is not satisfied in the whole space but only in the neighbor of a 
point. It is local transformation. If one extends the transformation to the maximal submanifold 
such that 7^ 0, it is possible to get troubles for the whole manifold due to possible different 
topologies II190L 

A point transformation = Q^{q) can depend on one (or more than one) parameter. Let 
us assume that a point transformation depends on a parameter e, i.e. = Q\q,€), and that it 
gives rise to a one-parameter Lie group. For infinitesimal values of e, the transformation is then 
generated by a vector field: for instance, as well known, d/dx represents a translation along the 
X axis, x{d/dy) — y{d/dx) is a rotation around the z axis and so on. In general, an infinitesimal 
point transformation is represented by a generic vector field on Q 



X = a^(q)|-. (4.29) 



The induced transformation (14.281) is then represented by 



d_ 



(4.30) 



is called the "complete lift" of X [I190II . A function /(q, q) is invariant under the transformation 
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if 



dX 







(4.31) 



where Lx<=f is the Lie derivative of /. In particular, if Lx'^C = 0, X'^ is said to be a symmetry for 
the dynamics derived by C 

In order to see how Noether's theorem and cyclic variables are related, let us consider a La- 
grangian C and its Euler-Lagrange equations 



d dC dC 



0. 



dX dq^ dq^ 

Let us consider also the vector field (14.301) . Contracting (14.321) with the a*'s gives 



(4.32) 



d dC dC^ Q 

dX dq'' dq^ 



Being 



a 



d_dC 

dX dcf 



dX \ dcf 



da^ \ dC 

IxJW 



from (|4.33l) . we obtain 



The immediate consequence is the Noether Theorei^ 
If LxC = 0, then the function 



is a constant of motion. 



(4.33) 



(4.34) 



(4.35) 



(4.36) 



In the following, with abuse of notation, we shall write X instead of X*^, whenever no confusion is possible. 
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Remark. Eq. (l4.36l) can be expressed independently of coordinates as a contraction of X by a 
Cartan one-form 



(4.37) 



For a generic vector field Y = y^d/dx^, and one-form /3 = I3idx\ we have, by definition, iy/S 
y^(3i. Thus Eq. (l4.36l) can be written as 



ix^c — So 



(4.38) 



By a point-transformation, the vector field X becomes 



d 



(4.39) 



We see that X is still the lift of a vector field defined on the "space of positions". If X is a symmetry 
and we choose a point transformation such that 



i^dQ^ = 1 ; ixdQ' = i ^ 1 



(4.40) 



we get 



(4.41) 



Thus is a cyclic coordinate and the dynamics can be reduced II1671I1891 . 
Remarks: 



1 . The change of coordinates defined by (|4.40l) is not unique. Usually a clever choice is very 
important. 



2. In general, the solution of (14.401) is not well defined on the whole space. It is local in the 
sense explained above. 
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3. It is possible that more than one X is found, say for instance Xi, X2. If they commute, i.e. 
[Xi, X2] = 0, then it is possible to obtain two cyclic coordinates by solving the system 

tx.dQ^ = 1; tx.dQ^ = 1; tx.dQ' = 0; 2^1; tx.dQ' = 0; « ^ 2 . (4.42) 

The transformed fields will be d/dQ^, d/dQ"^. If they do not commute, this procedure is 
clearly not applicable, since commutation relations are preserved by diffeomorphisms. Let 
us note that X3 = [Xi, X2] is also a symmetry, indeed, being Lx^^C = Lx^Lx^C—Lx^Lx^C, = 
0. If X3 is independent of Xi, X2, we can go on until the vector fields close the Lie algebra. 
The usual way to treat this situation is to make a Legendre transformation, going to the 
Hamiltonian formalism and to a Lie algebra of Poisson brackets. If we look for a reduction 
with cyclic coordinates, this procedure is possible in the following way: 

• we arbitrarily choose one of the symmetries, or a linear combination of them, and get 
new coordinates as above. After the reduction, we get a new Lagrangian £(Q); 

• we search again for symmetries in this new space, make a new reduction and so on 
until possible; 

• if the search fails, we try again with another of the existing symmetries. 

Let us now assume that C is of the form (I4.26h . As X is of the form (14.301) . LxC will be a 
homogeneous polynomial of second degree in the velocities plus a inhomogeneous term in the g*. 
Since such a polynomial has to be identically zero, each coefficient must be independently zero. If 
n is the dimension of the configuration space, we get l + n(n + l)/2 partial differential equations 
(PDE). The system is overdetermined, therefore, if any solution exists, it will be expressed in terms 
of integration constants instead of boundary conditions. It is also obvious that an overall constant 
factor in the Lie vector X is irrelevant. In other words, the Noether Symmetry Approach can be 
used to select functions which assign the models and, as we shall see below, such functions (and 
then the models) can be physically relevant. This fact justifies the method at least a posteriori. 

4.3 The Noether Approach for / — gravity in spherical symme- 
try 

Since the above considerations, if one assumes the spherical symmetry, the role of the ajfine pa- 
rameter is played by the coordinate radius r. In this case, the configuration space is given by 
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Q = {A, M, R} and the tangent space by T Q = {A, A', M, M', R, R'}. On the other hand, 
according to the Noether theorem, the existence of a symmetry for dynamics described by the 
Lagrangian (14.131) implies a constant of motion. Let us apply the Lie derivative to the (14.131) . we 
haveE : 



a ■ VgL + 2( Vaa ] L 



(4.43) 



which vanish if the functions a satisfy the following system 



a ■ V,L + 2(V,a)*L = a,^^ + 2^U^ = . (4.44) 

oil c/Qk 



The system (|4.44l) assumes the following explicit form 



From now on, q indicates the vector {A, M, R). 
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' rl^^h + M^hn^ =0 



A_ 

M 



(2 + MR)asfRR - '-f^fR 



fR + 



+T 



ai I cy dai I 2A da-j If \ A 



M 



M dM 



«3 _|_ A dctj \ f 



M 



a2{f-RfR)fR-T 



'^^ dM dA 



as 

+ 



9a 1 j_ A dai \ e 
M dA l-lR 



M(2 + MR)asfRR - 2«2/ij 



da2 
R dR 



+ 



+ (^2a2 + + 2A^ + Mff ^ a^fRR + 



Ml 



RRFl 



2A[{2 + MR)asfRR - (/ - RfR)a2]fRR + T 



+ (2ai + 2Aff + M|f + 2Af^j /rr + 2Aa3 W 



dai I A da2 \ f i 



dR 



|9aa I /M-Soi , 94^1 







(4.45) 



where T = (2 + MR)fR — Mf. Solving the system (14.451) means to find out the functions 
which assign the Noether vector. However the system (14.441) implicitly depends on the form of / 
and then, by solving it, we get also / theories compatible with spherical symmetry. On the other 
hand, by choosing the / form, we can explicitly solve (14.441) . As an example, one finds that the 
system (14.441) is satisfied if we chose 



a 



ai,a2,as) = i {3 — 2s)kA, —kM, kR 



(4.46) 
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with s a real number, k an integration constant and /o a dimensional coupling constano This 
means that for any f = exists, at least, a Noether symmetry and a related constant of motion 

So: 



Eq = a - Vg'L = 

= 2skMR^'''^ [2s + {s- l)MR] [{s - 2)RA' - {2s'^ - 3s + 1)AR'] . (4.47) 

A physical interpretation of Sq is possible if one gives an interpretation of this quantity in GR. In 
such a case, with s = 1, the above procedure has to be applied to the Lagrangian (14.171) . We obtain 
the solution 



a^^ = (-kA, kM) . (4.48) 

The functions A and M give the Schwarzschild solution (14.191) , and then the constant of motion 
acquires the form 



So = r3. (4.49) 

In other words, in the case of Einstein gravity, the Noether symmetry gives as a conserved quantity 
the Schwarzschild radius or the mass of the gravitating system. Another solution can be find out 
for R = Rq where Rq is a constant. In this case, the field equations (11.251) reduce to 

R,.u + k^g^^ = , (4.50) 
where ko = — |/(_Ro)///?(-Ro)- The general solutionis 



A(r) = -i^ = l + ^ + ^r^ M = r' (4.51) 
B{rj r 12 

with the special case 

^The dimensions are given by B.^^'^ in term of the Ricci scalar For the sake of simplicity we will put /o = 1 in 
the forthcoming discussion. 



4.3 The Noether Approach for / — gravity in spherical symmetry 83 



^W = 7^ = l + -' M = r\ R = 0. (4.52) 

B{r) r 

The solution (14.511) is the well known Schwarzschild-de Sitter one which is a solution in most 
of modified gravity theories. It evades the Solar System constraints due to the smallness of the 
effective cosmological constant. However, other spherically symmetric solutions, different from 
this, are more significant for Solar System tests. In the general case / = R^, the Lagrangian (14.131) 
becomes 

sR^'-^[2s + {s~l)MR] 
^ ~ M 

X [2(s - 1)M^A'R' + 2MRM'A' + 4(s - 1)AMM'R' + ARM'^] , (4.53) 

and the expression (|4.1 II) for B is 



s[2{s - l)M^A'R! + 2MRM'A' + 4(s - l)AMM'R! + ARM'"^] 
~ 2AMR[2s + (s - l)MR] ^ ' ' 

As it can be easily checked, GR is recovered when s = 1. Using the constant of motion (14.471) . we 
solve in term of A and obtain 



4s2-9s+5 



/Jl 2-s df 
2ksis-2)M[2s + (s^l)MR] ^ 

for s ^ 2, with ki an integration constant. For s = 2, one finds 



A = ^ (A 

12kr\4 + rm)RR' ' ^ ' ^ 

These relations allow to find out general solutions for the field equations giving the dependence of 
the Ricci scalar on the radial coordinate r. For example, a solution is found for 



3 = 5/4, M = r'^, R = 5r-^, 



(4.57) 
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obtaining the spherically symmetric metric 

ds^ = ^(k2 + Mdt^ - - ( — ^ I dr^ - r^dn , (4.58) 

with ^2 = ^ffc- It is worth noting that such exact solution is in the range of s values ruled out by 
Solar System observations, as pointed out in ni45[|146lll47ll . 



4.4 Perspectives of Noether symmetries approach 

In this chapter, we have discussed a general method to find out exact solutions in Extended Theories 
of Gravity when a spherically symmetric background is taken into account. In particular, we 
have searched for exact spherically symmetric solutions in /-gravity by asking for the existence 
of Noether symmetries. We have developed a general formalism and given some examples of 
exact solutions. The procedure consists in: i) considering the point-like / Lagrangian where 
spherical symmetry has been imposed; ii) deriving the Euler-Lagrange equations; Hi) searching 
for a Noether vector field; iv) reducing dynamics and then integrating the equations of motion 
using conserved quantities. Viceversa, the approach allows also to select families of / models 
where a particular symmetry (in this case the spherical one) is present. As examples, we discussed 
power law models and models with constant Ricci curvature scalar. However, the above method 
can be further generalized. If a symmetry exists, the Noether Approach allows, as discussed in 
§ 14. 2[ transformations of variables where the cyclic ones are evident. This fact allows to reduce 
dynamics and then to get more easily exact solutions. For example, since we know that f = W - 
gravity admit a conserved quantity, a coordinate transformation can be induced by the Noether 
symmetry. We ask for the coordinate transformation : 



L = L(g, g') = L(A, M, R, A', M\ R')^L = L(M, R, A', M\ R') , (4.59) 

for the Lagrangian (14.131) , where the Noether symmetry, and then the conserved quantity, corre- 
sponds to the cyclic variable A. If more than one symmetry exists, one can find more than one 
cyclic variables. In our case, if three Noether symmetries exist, we can transform the Lagrangian 
L in a Lagrangian with three cyclic coordinates, that is A = A{q), M = M{q) and R = R{q) 
which are function of the old ones. These new functions have to satisfy the following system 
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r (3-2.)Ai-MM + i?|| = i, 

< (4.60) 

[ (3-2.)Af -M|| + i?f = 0, 

with i = 2,3 (we have put k = 1). A solution of (14.601) is given by the set (for s ^ 3/2) 

A = + Fa{A^ M''-^ , A^i M^^) 

(4.61) 

q- = Fi{A^- M''^ A^ M^^) 
and if s = 3/2 

I = - In M + FA{A)GAiMR) 

(4.62) 

Qi = Fi{A)G,{MR) 

where Fa, Fi, Ga and Gi are arbitrary functions and 7]a, fli, and integration constants. 

These considerations show that the Noether Symmetries Approach can be applied to large 
classes of gravity theories. Up to now the Noether symmetries Approach has been worked out in 
the case of FRW - metric. In this chapter, we have concentrated our attention to the development 
of the general formalism in the case of spherically symmetric spacetimes. Therefore the fact that, 
even in the case of a spherical symmetry, it is possible to achieve exact solutions seems to suggest 
that this technique can represent a paradigmatic approach to work out exact solutions in any theory 
of gravity. At this stage, the systematic search for exact solution is well beyond the aim of this 
thesis. A more comprehensive analysis in this sense will be the argument of forthcoming studies. 
The results presented in this chapter point out that it does not hold in general for the specific / 
theories considered. However, the above technique could be a good approach to select suitable 
classes of theories where such a theorem holds. 



Chapter 5 



/ — gravity and scalar- tensor gravity: 
affinities and differences 

In the last years a very strong debate has been pursued about the Newtonian limit of HOG models. 
According to some authors the Newtonian limit of / - gravity is equivalent to the one of Brans- 
Dicke gravity with ubd = 0, so that the PPN parameters of these models turn out to be ill defined. 
In this chapter we show that this is indeed not true. We discuss that HOG models are dynamically 
equivalent to a O'Hanlon Lagrangian which is a special case of Scalar- tensor theory characterized 
by a self-interaction potential and that, in the low energy and small velocity limit, this will imply 
a non-standard behaviour [H]. This result turns out to be completely different from the one of a 
pure Brans-Dicke model and in particular suggests that it is completely misleading to consider the 
PPN parameters of this theory with lobd = in order to characterize the homologous quantities 
of /-gravity. 

By using the definition of the PPN-parameters 7 and (3 (12.251) in term of /-theories, we show 
that a family of third-order polynomial theories, in the Ricci scalar R, turns out to be compatible 
with the PPN - limit and the deviation from GR, theoretically predicted, can agree with experimen- 
tal data [A]. 

5.1 PPN — parameters in Scalar — Tensor and Fourth Order 
Gravity 

If one takes into account a more general theory of gravity, the calculation of the PPN - limit can be 
performed following a well defined pipeline shown in the §|23] which straightforwardly generalizes 
the standard GR. A significant development in this sense has been pursued by Damour and Esposito 
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- Farese [|23l I19U I192[ I193II which have approached to the calculation of the PPN-limit of scalar - 
tensor gravity by means of a conformal transformation (see § 11.31) to the standard Einstein frame. 
This scheme provides several interesting results up to obtain an intrinsic definition of 7, (3 in term 
of the non - minimal coupling function F{(f)). The analogy between scalar - tensor gravity and 
higher order theories of gravity has been widely investigated Il33ll43ll48l . 

Starting from this analogy, the PPN results for scalar - tensor gravity can be extended to HOG 
II137L In fact, identifying — )■ i? [|48l , it is possible to extend the definition of the scalar-tensor 
PPN - parameters (231 [HI to the case of HOG : 



. _ f" , .If f-f" .... 

7 f + 2/"2' ^ 4 V2/' + 3/"V (ii?- ^ 

In II137II , these definitions have been confronted with the observational upper limits on 7 and (3 
coming from Mercury Perihelion Shift [I195II and Very Long Baseline Interferometry [I196II . Ac- 
tually, it is possible to show that data and theoretical predictions from (|5.1I) agree in the limits 
of experimental measures for several classes of fourth order theories. Such a result tells us that 
extended theories of gravity are not ruled out from Solar System experiments but a more careful 
analysis of theories against experimental limits has to be performed. A possible procedure could be 
to link the analytic form of a generic fourth order theory with experimental data. In fact, the match- 
ing between data and theoretical predictions, found in [I137II . holds provided some restrictions for 
the model parameters but gives no general constraints on the theory. In general, the function / 
could contain an infinite number of parameters {i.e. it can be conceived as an infinite power series 
[|43l ) while, on the contrary, the number of useful relations is finite (in our case we have only two 
relations). An attempt to deduce the form of the gravity Lagrangian can be to consider the relations 
(15.11) as differential equations for /, so that, taking into account the experimental results, one could 
constrain, in principle, the model parameters by the measured values of 7 and 13. This hypothesis 
is reasonable if the derivatives of / function are smoothly evolving with the Ricci scalar. Formally, 
one can consider the r.h.s. of the definitions (15.11) as differential relations which have to be matched 
with values of PPN - parameters. In other words, one has to solve the equations (15.11) where 7 and 
(3 are two parameters. Based on such an assumption, on can try to derive the largest class of / 
- theories compatible with experimental data. In fact, by the integration of (15.11) . one obtains a 
solution parameterized by (3 and 7 which have to be confronted with the experimental quantities 
I3exp and '^exp- 



Assuming /' + 2/"^ 7^ and defining A 
for/: 



1-7 



27-1 



we obtain from (15.11 ) a differential equation 
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f"' = Af'. (5.2) 

The general solution of such an equation is a third order polynomial / = aR^ + hR^ + cR + d 
whose coefficients have to satisfy the conditions : a = 6 = c = and ci 7^ (trivial solution) 
or a = ^ and 6 = ± with c, d 7^ 0. Thus, the general solution for the non-trivial case, in 
natural units, reads 



/ 



1 


1-7 




1-7 


12 


27 - 1 




27-1 



R^ + cR + d. 



(5.3) 



It is evident that the integration constants c and d have to be compatible with GR prescriptions 
and, eventually, with the presence of a cosmological constant. Indeed, when 7 — )• 1, which implies 
f ^ c R + d, the GR - limit is recovered. As a consequence the values of these constants remain 
fixed (c= 1 and d = A, where A is the cosmological constant). Therefore, the fourth order theory 
provided by (15.31 ) becomes 



/: 



± 



1 


1-7 




1-7 


12 


27-1 




27 - 1 



R^ + R + A, 



(5.4) 



where we have formally displayed the two branch form of the solution depending on the sign of 
the coefficient entering the second order term. Since the constants a, b, c, d of the general solution 
satisfy the relation 3ac — b"^ = 0, one can easily verify that it gives : 



d'j 
dR 



d f 



112 



/± dRf' + 2f 



112 



0, 



(5.5) 



where the subscript refers the calculation to the solution (|5.4I) . This result, compared with the 
second differential equation (|5.1I) . implies 4(/3 — 1) = , which means the compatibility of the 
solution even with this second relation. 
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5.2 Comparing with experimental measurements 

Up to now we have discussed a family of fourth order theories (15.41) parameterized by the PPN - 
quantity 7; on the other hand, for this class of Lagrangians, the parameter /3 is compatible with GR 
value being unity. 



Mercury Perihelion Shift 


I27-/3- 1| < 3 X 10-=^ 


Lunar Laser Ranging 


4/3-7-3 = -(0.7 ±1) X 10-=^ 


Very Long Baseline Interf. 


I7 - 1| = 4 X 10"^ 


Cassini Spacecraft 


7-1 = (2.1 ±2.3) X 10-^ 



Table 5,1: A schematic resume of recent experimental constraints on the PPN - parameters. They are the 
perihelion shift of Mercury II195II . the Lunar Laser Ranging 111971 . the upper limit coming from the Very 
Long Baseline Interferometry II196I and the results obtained by the estimate of the Cassini spacecraft delay 
into the radio waves transmission near the Solar conjunction II198II . 

Now, the further step directly characterizes such a class of theories by means of the experi- 
mental estimates of 7. In particular, by fixing 7 to its observational estimate 'jexp, we will obtain 
the weight of the coefficients relative to each of the non-linear terms in the Ricci scalar of the La- 
grangian (15.41) . In such a way, since GR predictions require exactly 7exp = (3exp = 1, in the case 
of HOG, one could to take into account small deviations from this values as inferred from experi- 
ments. Some plots can contribute to the discussion of this argument. In figure [STl the Lagrangian 
(15.41) is plotted. It is parameterized for several values of 7 compatible with the experimental bounds 
coming from the Mercury perihelion shift (see Table 1 and [I195II ). The function is plotted in the 
range R> 0. Since the property f+{R) = —/_(—/?) holds for the function (15.41) . one can easily 
recover the shape of the plot in the negative region. As it is reasonable, the deviation from GR 
becomes remarkable when scalar curvature is large. 

In order to display the differences between the theory (15.41) and Hilbert-Einstein one, the ratio 
f / R is plotted in figure 15.21 Again it is evident that the two Lagrangians differ significantly for 
great values of the curvature scalar. It is worth noting that the formal difference between the PPN - 
inspired Lagrangian and the GR expression can be related to the physical meaning of the parameter 
7 which is the deviation from the Schwarzschild - like solution. It measures the spatial curvature 
of the region which one is investigating, then the deviation from the local flatness can be due to 
the influence of higher order contributions in Ricci scalar. On the other hand, one can reverse 
the argument and notice that if such a deviation is measured, it can be recast in the framework 
of HOG, and in particular its "amount" indicates the deviation from GR. Furthermore, it is worth 
considering that, in the expression (|5.4I) . the modulus of the coefficients in 7 (i.e. the strength of 
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f (R) 




20 40 60 80 100 
R 



Figure 5.1: Plot of the two branch solution provided in (15.41) . The /+ (dotted line) branch family is up to 
GR solution (straight line), while the one indicated with /_ (dotted-dashed line) remains below this line. 
The different plots for each family refer to different values of 7 fulfilling the condition I7 — 1| <10^^ and 
increased by step of 10^^. 

the term) decreases by increasing the degree of R. In particular, the highest values of cubic and 
squared terms in R are, respectively, of order 10^^ and 10^^ (see figure [53l) then GR remains a 
viable theory at short distances (i.e. Solar System) and low curvature regimes. 

5i I 

4 



F (R) 




20 40 60 80 100 
R 



Figure 5.2: The ratio f /R. It is shown the deviation of the HOG from GR considering the PPN - limit. 
Dotted and dotted-dashed lines refer to the /+ and /_ branches plotted with respect to several values of 7 
(the step in this case is 2.5 x 10^^). 

A remark is in order at this point. The class of theories which we have discussed is a third order 
function of the Ricci scalar R parameterized by the experimental values of the PPN parameter 7. 
In principle, any analytic / can be compared with the Lagrangian (15.41 ) provided suitable values 
of the coefficients. However, more general results can be achieved relaxing the condition (3 = 1 
which is an intrinsic feature for (|5.4I) (see for example 1113711 ). These considerations suggest to take 
into account, as physical theories, functions of the Ricci scalar which slightly deviates from GR, 
i.e. / = foR^^~^^^ with e a small parameter which indicates how much the theory deviates from 
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Figure 5.3: Plot of the modulus of coefficients: a (i?^) (line) and b (R^) (dashed line). The choice of 
plotting the modulus of coefficients is the consequence of two solution for /. 

GR [I145L In fact, supposing e sufficiently small, it is possible to approximate this expression 

\R\^'^'^ ~ \R\ (^l + e\n\R\ + f!]!^ + ...y (5.6) 

This relation can be easily confronted with the solution (15.41) since, also in this case, the corrections 
have very small "strength". 

We can conclude this paragraph having shown how a polynomial Lagrangian in the Ricci scalar 
R, compatible with the PPN - limit, can be recovered in the framework of HOG. The approach is 
based on the formulation of the PPN - limit of such gravity models developed in analogy with 
scalar-tensor gravity [1137L In particular, considering the local relations defining the PPN fourth 
order parameters as differential expressions, one obtains a third-order polynomial in the Ricci 
scalar which is parameterized by the PPN - quantity 7 and compatible with the limit (3 = 1. The 
order of deviation from the linearity in R is induced by the deviations of 7 from the GR expectation 
value 7 = 1. Actually, the PPN parameter 7 may represent the key parameter to discriminate 
among relativistic theories of gravity. In particular, this quantity should be significatively tested 
at Solar System scales by forthcoming experiments like LATOR [11991 . From a physical point 
of view, any analytic function of R, by means of its Taylor expansion, can be compared with 
(15.41) . Therefore, a theory like / = foR^^~^^\ indicating small deviations from standard GR, is 
in agreement with the proposed approach, so, in principle, the experimental 7 could indicate the 
value of the parameter e. In conclusion, one can reasonably state that generic fourth-order gravity 
models could be viable candidate theories even in the PPN - limit. In other words, due to the 
presented results, they cannot be a priori excluded at Solar System scales. 
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5.3 Newtonian limit of / — gravity by O'Hanlon theroy anal- 
ogy 

Recently several authors claimed that HOGs and among these, in particular, HOG models are 
characterized by an ill defined behaviour in the Newtonian regime. In particular, in a series of 
papers [|135[|143[[T44lll49[l200ll201ll it is addressed that Post-Newtonian corrections of the grav- 
itational potential violate experimental constraints since these quantities can be recovered by a 
direct analogy with Brans-Dicke gravity [|25l simply supposing the Brans - Dicke characteris- 
tic parameter cobd vanishing. Actually despite the calculation of the Newtonian and the post- 
Newtonian limit of / - theory, performed in a rigorous manner, have showed that this is not the 
case [|137lll511ll85[|202ll . it remains to clarify why the analogy with Brans - Dicke gravity seems 
to fail its predictions. The issue is easily overcame once the correct analogy between /-gravity and 
the Brans-Dicke model is taken into account. 

In literature, it is suggested that HOG models can be rewritten in term of a scalar-field La- 
grangian non minimally coupled with gravity but without any kinetic term by fact implying cobd = 
(as shown in the chapter 1). Actually, the simplest case of scalar - tensor gravity models has been 
introduced some decades ago by Brans and Dicke in order to give a general mechanism capable 
of explaining the inertial forces by means of a background gravitational interaction. The explicit 
expression of such gravitational action is (|1.30l) . while the general action of /-gravity is (|1.24l) . As 
said above, /-gravity can be recast as a scalar-tensor theory by introducing a suitable scalar field (j) 
which nonminimally couples with the gravity sector. It is important to remark that such an analogy 
holds in a formalism in which the scalar field displays no kinetic term but is characterized by means 
of a self-interaction potential which determines the whole dynamics (O'Hanlon Lagrangian) [l33l . 
We can resume the actions as follow 



<pR + V{(p) + XCr 



(5.7) 



This consideration, therefore, implies that the scalar field Lagrangian equivalent to the purely 
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geometrical /-gravity turns out to be quite different with respect to the ordinary Brans - Dicke 
definition (11.301) . This point represents a crucial aspect of our analysis. In fact, as we afterwards 
will see, such a difference will imply completely different results in the Newtonian limit of the two 
models and, consequently, the impossibility of extend predictions from the PPN approximation of 
Brans-Dicke models to /-gravity. Considering natural units, the O'Hanlon Lagrangian [I203II is the 
third of (15.71 ) or (11.421 ). The field equations are obtained by varying Eq. (11.421 ) with respect to both 
and (p which now represent the dynamical variables (the same field equations are given setting 
uj{(f)) = and = d) in the (11.311) ). Thus, one obtains 



< R+^ = o (5.8) 

where the second line of (|5.8I) the field equation for 0. While the third equation is a combination 
of the trace of the first one and of the second one. The two schemes can be mapped one into the 
other considering the following equivalences 



< V{<P) = f-f'R (5.9) 

^ 0^™ - 2Vi<P) = f'R - 2/ 

where we are the Jacobian matrix of the transformation <^=^ /' is non-vanishing. Henceforth we 
can consider instead of (11.251) - (|1.26l) a new set of field equations determined by the equivalence 
of /-gravity with the O'Hanlon approach: 



<PR,, + I iv{<f)) + 0^ j g,, - 0;^, = XJ:,, 



(5.10) 



□ + 
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where S^^ = T^^ - \Tgf,^,. 

Let us, now, calculate the Newtonian limit of Eqs. (I5.10I ). Considering the perturbations of 
the metric tensor g^^i, and of the scalar field cf) with respect to a background value, we search for 
solution at the 0(2) - order in term of the metric entries and of the scalar field itself (see § 12.31) : 



\ 



(5.11) 



+ 



(5.12) 



where we neglected the vectorial component in the metric. The differential operators turn out to 
be approximated as 



U^dl-^ and V^V. ^ dl, . 



(5.13) 



Actually in order to simplify calculations we can exploit the intrinsic gauge freedom intrinsic in the 
metric definition. In particular, we choose the harmonic gauge (12.431) and the expressions of Ricci 
tensor components given by Eqs. (12.531 ). In relation with the adopted approximation we coherently 
develop the self-interaction potential at second order. In particular, the quantities in (15.101) read : 



0^ - 2^(0) ^ 0(0)^ - 2n0(°))+ 



0(2) 



0(2) 



(5.14) 



The field equations (15.101) . solved at 0-th order of approximation, provide the two solutions 



V(^(o,)=o and '^^^ 



(5.15) 



which fix the 0-th order terms in the development of the self-interaction potential; therefore we 
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have 



(5.16) 

0^ - 2^(0) ^ = 3AV(^) 

where constant factors 0*^°^ have been condensed within the quantity 3A^. Such a constant 

can be easily interpreted as a mass term as will become clearer in the following. Now, taking into 
account the above simplifications, we can rewrite field equations (|5.10l) at the at 0(2) - order in the 
form : 



A.;? = ^4°)-A^^, (5.17) 



A (2) _ 2 A' (0) 2 lo'^Ji 



A0(2) - a20(2) = -!^rW . (5.19) 

3 



The scalar field solution can be easily obtained from the (15.191 ) as 



3 J (27r)3/2 k2 + A2 



while for (jfff ■* and g^f we have 



g)f\^) = — / d'x.' , " X + ^ / c^^x V ^ \ , (5.21) 



'The factor 3 is introduced to simplify an analogous factor present in the field equations ( 15.101 ). 
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27r0(o) J |x-x'| 47r0(o) J |x - x' 



2 

+ 



0(0) 



X2 " ^ ^ ' V X2 



X 







(5.22) 



The above three solutions represent a completely general result. In particular adopting the trans- 
formation rules (15.91) . one can straightforwardly obtain the solutions in the pure / - scheme. 

Let us analyze the above results with a simple example. We can consider a HOG Lagrangian 
of the form / = aiR + a2-R^ so that the "dummy" scalar field reads = ai + 2a2-R. The self- 
interaction potential turns out the be = —^"^-^^ satisfying the conditions V{ai) = and 
V'{ai) = 0. In relation with the definition of the scalar field, we can opportunely identify oi with 
a constant value (f)^^^ = ai. Furthermore, the scalar field "mass" can be expressed in term of the 
Lagrangian parameters as = ^ = Since the Ricci scalar at the second order 
reads 



|!1._L/^I!^, (5.23, 

if we consider a point-like mass M, the energy-momentum tensor components become respectively 
Ttt = p, T p while p = M5{x), therefore we obtain 



The immediate consequence is that the solution for the scalar field </> at second order is 



= 2a2i?(^) = a/I^^ (5.25) 
V 2 3 |x| 

while the complete scalar field solution up to the second order of perturbation is given by 

Once the behavior of the scalar field has been obtained up to the second order of perturbation, in 
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the same way, one can deduce the expressions for g^^^ and glj\ where Ti[^^ = |p and 
hp6ij = ^I][f^6ij. As matter of fact the metric solutions at the second order of perturbation are 



gtt 



9ij 



I 2_£g_ _ /TT 



3a |x| 



7f 1 rge 



-A|x| 



2 3a 



|x| 



1 + - 
^ ^ 3a |x| 



Els 
2 3a 



3a 



1 



|x| ' A|x|: 



|x| 



A2|xP A|x|5 



+ 



+ 



3 

A^l^ 



p-A|x| _ _3 

A^ 



-A|x| 



A2|xP 



(5.27) 



This quantity, which is directly related to the gravitational potential, shows that the gravitational 
potential of the O'Hanlon Lagrangian is non-Newtonian like. Such a behavior prevents from ob- 
taining a natural definition of the PPN parameters as corrections to the Newtonian potential. As 
matter of fact since it is indeed not true that a generic /-gravity model corresponds to a Brans- 
Dicke model with cubd = ^ coherently to its Post-Newtonian approximation. In particular it 
turns out to be wrong considering the PPN parameter 7 = (see, for example, II129II ') of 

Brans - Dicke gravity and evaluating this at to bd = so that one has 7 = 1/2 as suggested in 
[fT35l[T43l[T49ll200l . 

Differently, because of the presence of the self-interaction potential V"(0), in the O'Hanlon 
Lagrangian, a Yukawa like correction appears in the Newtonian limit appears. Such a correction 
in a completely different way even at the post-Newtonian limit. As matter of fact, one obtains 
a completely different gravitational potential with respect to the ordinary Newtonian one and as 
matter of fact the fourth order corrections in term of the v/c ratio (Newtonian level), have to be 
evaluated in a complete new general way. In other words, considering a Brans-Dicke Lagrangian 
and an O'Hanlon one, despite their similar structure, will imply completely different predictions 
in the weak field and small velocity limits. Such a result represents a significant argument against 
the claim that HOG models can be ruled out only on the bases of the analogy with Brans-Dicke 
PPN parameters. 

an important consideration is in order now. The definition of PPN-parameters 7 and (3, in 
the GR realm, is intended as a correction to the Newtonian-like behaviour of the gravitational 
potentials (|2.25l) . In particular, the PPN parameter 7 is related to the second order correction of to 
the gravitational potential while (3 is linked with the fourth order level of perturbation. Actually, if 
we consider the limit f ^ R, from Eqs. (|5.27l) . we have 
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9tt 



9ij 



2_rg_ 
3a |x| 



1 + 



3a |x| ; 



(5.28) 



Since a is an arbitrary constant, in order to match the Newtonian gravitational potential of GR, we 
should fix a = 2/3. This assumption implies 



9tt- i- 1^ 

< (5.29) 

9ij ~ ~^ 2 lx| ^ 

which suggest that the PPN parameter 7, in this limit, results 1/2 which is in striking contrast 
with GR predictions (7 ~ 1). Such a result is in reality not surprising. In fact, the GR limit 
of the O'Hanlon Lagrangian requires (p ~ const and V^(0) — ?■ but such approximations induce 
mathematical inconsistencies in the field equations of / - gravity once these have been obtained by 
a general O'Hanlon Lagrangian. In reality this is a general issue of O'Hanlon Lagrangian. In fact 
it can be demonstrated that the field equations (15.101) do not reduce to the standard GR ones (for 
V{(f)) — > and ~ const) since we have : 



{^^^'^ ai ^M!^ 
(5.30) 
= fT 

But S^j, components read T^u = |p and Sjj = |p5ij = ^^tt^ij in place of Su = \p and Sij = 
\p6ij = SttSij as usual, while the GR field equations are the (12.581) . Such a pathology is in 
order even when the GR limit is performed from a pure Brans - Dicke Lagrangian. In such a 
case, in order to match the Hilbert - Einstein Lagrangian, one needs (p ~ const and cubd = 
0, the immediate consequence is that the PPN parameter 7 turns out to be 1/2, while it is well 
known that Brans - Dicke model fulfils low energy limit prescriptions in the limit a; — )■ 00. Even 
in this case, the problem, with respect to the GR prediction, is that the GR limit of the model 
introduces inconsistencies in the field equations. In other words, it is not possible to impose the 
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same transformation which leads the Brans-Dicke theory into GR at the Lagrangian level on the 
solutions obtained by solving the field equations descending from the general Lagrangian. The 
relevant aspect of this discussion is that considering a / model, in analogy with the O'Hanlon 
Lagrangian and supposing that the self-interaction potential is negligible, introduces a pathological 
behaviour of the dynamical solutions and induces to obtain a PPN parameter 7 = 1/2. This is 
what happens when an effective approximation scheme is introduced in the field equations in order 
to calculate the weak field limit of HOG by means of Brans-Dicke model. Such a result seems, 
from another point of view, to enforce the claim that HOG models have to be carefully investigated 
in this limit and their analogy with scalar-tensor gravity should be opportunely considered. 



5.4 Differences of a generic scalar — tensor theory in the Jor- 
dan and Einstein frames 



up to now Since along the chapter we have discussed the weak field and small velocity limit 
of HOG models in term of Brans-Dicke like Lagrangian remaining in the Jordan frame. There 
we show what are the predictions of the weak field and small velocity limit when a conformal 
transformation (|1.32l) is applied on the O'Hanlon Lagrangian. In other words, we discuss HOG 
models conformally transformed in the Einstein frame. The generic scalar-tensor action Ajp in 
the Jordan frame (|1.29l) is linked to the generic action in the Einstein frame (11.351) via the 
transformations (11.361) between the quantities in the two frames. In the case of the O'Hanlon 
theory in the Jordan frame, (11.421) . i.e. F((f)) = and co{(f)) = 0, the action (|1.35l) in the Einstein 
frame is simplified and the transformation between the two scalar fields reads 

fiMV = -Y^- (5.31) 
If, now, we suppose i^{^p) = —i^o < we have 

= ke^^"", (5.32) 

where Y = ^ ^ and k is an integration constant. We obtain the transformed of (11.421 ) in the 
Einstein frame is 
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^EF 



k 



The field equations on the other side are 



(5.33) 



(5.34) 



where the matter tensor, which now coupled with the scalar field </?, in the Einstein frame reads 



and 



-1 Sj^gCm) _ ^2y^ 



^ ^9 



e"' '^9 pa 



A2y 



k^ 



(5.35) 



(5.36) 



Actually, in order to calculate the weak field and small velocity limit of the model in the Einstein 
frame, we can develop the two scalar fields at the second order ~ 0^°) + 0^^) and </? ~ </7^°^ + </?^^^ 
with respect to a background value. This choice gives the relations : 



(0) 



(2) 



0(2) 



(5.37) 
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Let us consider the conformal transformation gf^,^ = jgf_iL, (|1.32l) . From this relation and consid- 
ering the (15.321) one obtains, if 0^°-' = A, that 



~(2) _ (2) 0(2) 

9tt — 9tt + ^(0) 



(5.38) 



(2) _ (2) _ 0(2) 



As matter of fact, since ^Jf^ = 2$^^, ^Jf = 2^-^^ 5^ and ^Jf^ = 2^^^,glf = 2^^^ from 
(15.381) it descends a relevant relation which links the gravitational potentials of Jordan and Einstein 
frame : 



^ ^ 20(0) ^ =^ 2 

(5.39) 

^ 20(0) ^ +2^ 

If we introduce the variations of two potentials: A$ = ^•^^ ~ and = ^-^^ — '^^^ we 
obtain the most relevant result of this paragraph: 




A$ = - AM/ = - ^ = oc a2 (X /"(O) . (5.40) 

From the above expressions, one can notice that there is an evident difference between the 
behaviour of the two gravitational potentials in the two frames. Such a result suggests that, at the 
Newtonian level, it is possible to discriminate between the two mathematical frame thus one can 
deduce what is the true physical one. In particular, once, the gravitational potential is calculated in 
the Jordan frame and the dynamical evolution of is taken into account at the suitable perturbation 
level, these can be substituted in the first of (15.391) so that to obtain its Einstein frame evolution. 
The final step is that the two potentials have to be matched with experimental data in order to 
investigate what is the true physical solution. A similar result has been provided in a recent paper 



Chapter 6 



The Newtonian limit of Fourth Order 
Gravity theory 

The Newtonian limit of HOG is worked out discussing its viability with respect to the standard re- 
sults of GR. We exclusively investigate the limit in the metric approach, refraining from exploiting 
the formal equivalence of higher - order theories by considering the analogy with specific scalar 
- tensor theories, i.e. we work in the Jordan frame in order to avoid possible misleading inter- 
pretations of the results. Considering the Taylor expansion of a generic /-gravity, it is possible 
to obtain general solutions in term of the metric coefficients up to the third order of approxima- 
tion. Furthermore, we show that the Birkhoff theorem is not a general result for /-gravity since 
time-dependent evolution for spherically symmetric solutions can be achieved depending on the 
order of perturbations [C], [G] . Furthermore we provide explicit solutions for several different 
types of Lagrangians containing powers of the Ricci scalar as well as combinations of the other 
curvature invariants [E]. In particular, we develop the Green function method for fourth - order 
theories in order to find out solutions. Finally, the consistency of the results with respect to GR is 
discussed. In particular, the solution relative to the gu component gives a gravitational potential 
always corrected with respect to the Newtonian one of the linear theory f = R. 

6.1 The Newtonian limit of / — gravity in spherically symmet- 
ric background 

Exploiting the formalism of Newtonian and post-Newtonian approximation described in paragraph 
(12.31) , we can develop a systematic analysis in the limits of weak field and small velocities for the 
/-gravity. We are going to assume, as background, a spherically symmetric spacetime and we are 
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going to investigate the vacuum case. Considering the metric (11.891 ). we have, for a given g^^ 



{ gu{t,r) ~ l + gif\t,r) + g]:t\t,r) 
9rr{t,r) ~ -l+gi^\t,r) 
geeit.r) = -r^ 



while considering Eqs. (12.361) : 



(6.1) 



The determinant reads 



4 - 0^') 

L ijrr 



(6.2) 



g ^ sin^ ^^{-1 + \g^ - g^^] + g^ - g^^]} 



(6.3) 



The Christoffel symbols (12.391) are 



r(3) 



r(3) 



tr 



r(3) 



yrr,t 



yrr,t 



(2) (2) (2) , (4) 

^ tt'^'- tt — 2 2 



, ^^ . ^^- (4) „(2) (2) 

r(2) _|_ F(4)* _ 3tt,T _|_ 9u,r-9tt 9tt,r 



(2) (2) (2) 

yrr,r yrr yrr,r 



r (2) (2)2 

ee ~ ' ' yrr ' grr 



(6.4) 



while the Ricci tensor component, (12.401 ). are 
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(2) 



(4) 



(3) 
tr 



d(2) 
ihrr- 



2r 



Ar 



(2) 

9rr,t 



(6.5) 



2r 



(2) _ 2gff+r[g^,'^,+g^^.j,] 



p(2) _ .,;^2 /ip(2) 
it - - — Sm dKaa 



and, finally, the Ricci scalar expression is 



(6.6) 



In order to derive the Newtonian and post-Newtonian approximation for a generic function /, 
one should specify the / - Lagrangian into the field equations (|1.25l) . This is a crucial point because 
once a certain Lagrangian is chosen, one will obtain a particular approximation referred to such 
a choice. This means to lose any general prescription and to obtain corrections to the Newtonian 
potential, $(x), which refer "univocally" to the considered / - function. Alternatively, one can 
restrict to analytic / - functions expandable with respect to a certain value R = Rq = constant or, 
at least, its non-analytic part, if exists at all, goes to zero faster than i?", with n > 2 at i? — )■ 0. In 
general, such theories are physically interesting and allow to recover the GR results and the correct 
boundary and asymptotic conditions. Then we assume 
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Ro 



(6.7) 



One has to note that the development (|6.7I) . also if similar to (13.271) . is very different from the one 
in the Chapter [S] in fact R^^^ is a general space-time function linked to the background metric g'^l^} 
(13.251) . while here Rq is a constant value of scalar curvature. Besides the coefficients /o, /i, /2, /s 
are not proportional, respectively, to zero-th, first, second, third coefficient of Taylor development 
of /. In fact, we have 



/o = /(i?o) - i?o/'(i?o) + \Rlf"m - lRlf"'{Ro 



fi = f'iRo) - Rof'iRo) + \Rir"m 



/2 = kf'iRo) ~ kRof"'{Ro) 



(6.8) 



If we consider a flat background, then Rq = and the coefficients /o, /i, /2, /s are the terms of 
Taylor series. But if we are finding the solutions at Newtonian and (possibility) post-newtonian 
level we have to consider a vanishing background scalar curvature. It is possible to obtain the 
Newtonian and post-Newtonian approximation of / - gravity considering such an expansion (|6.7I) 
into the field equations (11.251) and expanding the system up to the orders 0(0), 0(2), 0(3) and 
0(4). This approach provides general results and specific (analytic) Lagrangians are selected by 
the coefficients fi in (16.71) . Developing the equations in the case of vanishing matter, i.e. T^,^ = 0, 
we have 



r Hj^'j = 0, = 



Hi^J = 0, ^ 



Hi^J = 0, ^ 



(6.9) 



ifit^ = 0, = 
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and, in particular, from the 0(0) order approximation, one obtains 



/o = 0, (6.10) 

which trivially follows from the above assumption that the space-time is asymptotically Minkowskian 
(asymptotically flat background). This result suggests a first consideration. If the Lagrangian is 
developable around a vanishing value of the Ricci scalar (Rq = 0) the relation Ii6.10\l will imply 
that the cosmological constant contribution has to be zero whatever is the f - gravity theory. 

If we now consider the 0(2) - order approximation, the equations system (16.91 ). in the vacuum 
case, results to be 



< ^fj? - r[f,rR^'^ - f^g^l - f,gi% + ^f2Rf + ^f2rRfr] = (6.11) 
/iri?(2) + Qf^[2Rf + rRfr] = 

^ 2gi? + r[2gSl - rR^'^ + 2gfl + rg^^^„\ = 

The last equation of the system (16.1 II) is the definition of Ricci scalar (|1.5I) at 0(2) - order. The 
trace equation (the fourth line in the (16.1 II) ). in particular, provides a differential equation with 
respect to the Ricci scalar which allows to solve, if sign[/i] = — sign[/2], the system (16.111) at 0(2) 
- order : 



91? — ^0 — + 



3A Ar 



6A2 Ar 



, „{2) _ _5i S2{t) Ar+l ^-Ar , ^ajt) Ar-l ^Ar 

\ yrr — 3A Ar ^ 6A2 Ar ^ 



(6.12) 



2A r 



where 
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(6.13) 



with the dimension of length~^. Let us notice that the integration constant 5q has to be dimension- 
less, 5i has the dimension of length, while the time - dependent functions 82 and ^3, respectively, 
have the dimensions of length^^ and length'^. The functions 5i{t) (i = 2, 3) are completely ar- 
bitrary since the differential equation system (16.231) contains only spatial derivatives. Besides, the 
integration constant 5o can be set to zero, as in the theory of the potential, since it represents an 
unessential additive quantity. When we consider the limit / — )• _R, in the case of a point-like source 
of mass M, we recover the perturbed version of standard Schwarzschild solution (|2.8I) at 0(2) - 
order with 5i = Vg. In order to match at infinity the Minkowskian prescription of the metric, we 
discard the Yukawa growing mode present in (16.121) . we have : 



R 



^ fir ^ 



S2{t) 



3A Ar 



]_ _|_ JIs. _|_ ^ijt) Ar+l ^-Ar 
fir 3A Ar 



dr"^ — r'^dVL 



(6.14) 



At this point one can provide the solution in term of the gravitational potential. In such a case, 
we have an explicit Newtonian-like term into the definition, according to previous results obtained 
with less rigorous methods Il32lll32ll . The first of (16.121) provides the second order solution in term 
of the metric expansion (see the definition (16.11 )). but, this term coincides with the gravitational 
potential at the Newtonian order (12.651) . In particular, since we have = 1 + 2$, the gravitational 
potential of a HOG theory, analytic in the Ricci scalar R, is 



^_GM^S^t^ (5.15, 

jir OA Ar 

As first remark, one has to notice that the structure of the potential (|6.15l) . for a given / - theory, 
is determined by the parameter A, (16.131 ). which depends on the first and the second derivative 
of /, once developed around a vanishing value of Ricci scalar. The potential (16.151 ) is coherent 
with respect to the results in [|132[ 1 17211 . obtained for higher order Lagrangians as F = R + 
Sfc=o ^kR □'^-R- In this last case, it is demonstrated that the number of Yukawa corrections to the 
gravitational potential was strictly related to the order of the theory. 
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From (16.141) one can notice that the Newtonian limit of any analytic / - theory is related only 
to the first and second term of the Taylor expansion (16.71) of the given theory. In other words, the 
gravitational potential is always characterized by two Yukawa corrections and only the first two 
terms of the Taylor expansion of a generic f Lagrangian turn out to be relevant. This is indeed a 
general result. 

Let us now consider the system (16.91) at third order contributions. The first important issue is 
that, at this order, one has to consider even the off-diagonal equation: 

f,gl% + 2f2rR^} = 0, (6.16) 

(2) 

which relates the time derivative of the Ricci scalar to the time derivative of grr . From this relation, 
it is possible to draw a relevant consideration. One can deduce that, if the Ricci scalar depends on 
time, so it is the same for the metric components and even the gravitational potential turns out to be 
time-dependent. This result agrees with the analysis provided in the Chapter [3] where a complete 
description of the weak field limit of HOG has been provided in term of the dynamical evolution 
of the Ricci scalar. Moreover it was been demonstrated that if one supposes the time independent 
Ricci scalar (13.211) . static spherically symmetric solutions result of the form (13.241) . 

Eq. (16.161) confirms this result and provides the formal theoretical explanation of such a behav- 
ior. In particular, together with Eqs. (16.141) . it suggests that if one considers the problem at a lower 
level of approximation (i.e. the second order) the background spacetime metric can have static 
solutions according to the Birkhoff theorem; this is no more verified when the problem is faced 
with approximations of higher order. In other words, the debated issue to prove the validity of the 
Birkhoff theorem in higher order theories of gravity, finds here its physical answer. In the Chapter 
[3] and here, the validity of this theorem is demonstrated for / - theories only when the Ricci scalar 
is time independent or, in addition, when the solutions of Eqs. (|1.25l) are investigated up to the 
second order of approximation in the metric coefficients (16. 1|) . Therefore, the Birkhoff theorem is 
not a general result for HOG but, on the other hand, in the limit of small velocities and weak fields 
(which is enough to deal with the Solar System gravitational experiments), one can assume that 
the gravitational potential is effectively time independent according to A6.14\) and A6.15\) . 

The above results fix a fundamental difference between GR and HOG theories. While, in GR, a 
spherically symmetric solution represents a stationary and static configuration difficult to be related 
to a cosmological background evolution, this is no more true in the case of HOG. In the latter case, 
a spherically symmetric background can have time-dependent evolution together with the radial 
dependence. In this sense, a relation between a spherical solution and the cosmological Hubble 
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flow can be easily achieved. 

From the system (16.91) . one can notice that the general solution is characterized only by the first 
three orders of the / expansion (16.71) . Such a result is in agreement with the / reconstruction which 
can be induced by the post-Newtonian parameters adopting a scalar-tensor analogy as discussed in 
Chapter [51 



6.1.1 Newtonian and post — Newtonian limit in the harmonic gauge 

Up to now, the discussion has been developed without any gauge choice. In order to overcome 
the difficulties related to the nonlinearities of calculations, we can work considering some gauge 
choice obtaining less general solutions for the metric entries. If we consider the gauge (12.431) we 
can use the Ricci tensor components (12.531 ) and the Ricci scalar expression (12.541 ). The gauge 
choice does not affect the Christoffel symbols. Thus, by solving the field equations (|1.25l) one 
obtains 



9u,^Jt,r) = 1 + + - + ^^^^^^ - 2/tf)/tf ^ 



(6.17) 



where the constants fcj are relative to approximation level 0(2), while k\ to 0(4). The Ricci 
scalar is zero both at 0(2) and at 0(4) approximation orders. 

Eqs. (16.171) suggest some interesting remarks. It is easy to check that the GR prescriptions 
are immediately recovered for a particularly choice of integration constants. The grr component 
displays only the second order term, as required by a GR - like behavior, while the component 
shows also the fourth order corrections which determine the second post-Newtonian parameter (3 
[I129II . It has to be stressed here that a full post-Newtonian formalism requires to take into account 
matter in the system (16.91 ): the presence of matter links the second and fourth order contributions 
in the metric coefficients [I129I . 



6.2 The Newtonian limit of quadratic gravity 

Since terms resulting from _R" with n > 3 do not contribute in the Newtonian limit, as seen 
previously, the most general choice for the Lagrangian is 
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f = aiR + a2R^ + a^Ra^R"^ (6.18) 

where oi, a2, ajll are constants. We have to note that the field equations (|1.25|) (with / = R^), 
(ll.74|) and (11.781) satisfy, in four dimensions, the condition 



Ric 



jRie 



0, 



(6.19) 



then only two of the three expressions are independent [l35ll . Such a quantity is related to the Gauss- 
Bonnet topological invariant ll36l . We can consider Eq. (16.181) as the most general quadratic theory 
of gravity. The field equations of (16.181) are a linear combination of (11.101) . (11.251) and (11.741) . that 



IS 



(6.20) 



If we introduce the generalization of the gravitational potentials in the isotropic metric (11.881) by 
the quantities $ and ^ linked 
(16.201) in the Newtonian limit: 



(2) (2") 

the quantities $ and ^ linked to gl^ and g], , we can investigate the solution of field equations 



1 + 2$ 



de 



1-2^ 



dijdx^dx^ 



(6.21) 



Up to the Newtonian order the left-hand side of the field equations (I1.10|) . (|1.25l) . (11.741) and (11.781) 

are 



1. for the GR - theory: 



G.,~Glf = A($-vl/)5,^.-(<|.-xl/),, 



(6.22) 



2. for R^ - theory: 



'Note that [a2] — [03] — length^ and [ai] ~ length^ 
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4A2($-2^) 

A2(2^ - + (A$ - 2A^),y 



(6.23) 



3. for RauR"^ - theory: 



[Ric 
Ht 



Ht 



TTRic ^ rrRic{2) 



2A2($ - ^) 



A2(3^ - ^)6ij + (A$ - 3A^), 



4. for Rc^p^sR"^^^ - theory: 



(6.24) 



f HS'' ~ i/f = 4A2$ 



(6.25) 



If we take into account the results (16.221) . (|6.23l) and (16.241) for the geometric side and the results 
(12.571) for the matter side, the explicit form of the field equations (16.201) up to the Newtonian order 
is 



2ai A^ - 2(4a2 + ag) A^^ + 2(2a2 + ag) A^^ = . 



A 



ai(^ - $) + (4a2 + a3)A<l> - (802 + 3a3)A^ 



ai(^ - $) + (4a2 + a3)A$ - (802 + 3a3)A^ = 
By introducing two new auxiliary functions ($ and the equations (16.261) become 



(6.26) 
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_2a2_^2^ — ^"3(3a2+a3) _ Aa2+ci3 

2a2+a-i 01(202+03) 202+03 



202+03 r- 



(6.27) 



A 



$ + A^ 



where $ and \I/ are linked to $ and ^ via 



_ (8o2+3o3)#+oi1' 
201(202+03) 



^ _ _ (4o2+a3)*+a]_£ 
201(202+03) 



(6.28) 



Obviously we must require ai(2a2 + 03) 7^ 0, which is the determinant of the transformations 
(16.281) . Let us introduce the new function S defined as follows: 



S := $ + A^. (6.29) 

A this point, we can use the new function S to uncouple the system (16.261) . With the choice 
$ = S — A^, it is possible to rewrite equations (|6.26|) as follows 



203(302+03) ^3^ ^ |02±a3 ^2^ _ = Xp+T, 

01(202+03) 202+03 202+03 1^ ' 

(6.30) 

A S 5ij - E^ij = 

where r = p2±B3/^^ + 203(302+03) ^2g_ interested in the solution of (l627l) in terms of the 

202+03 01(202+03) "I " 

Green function ^(x, x') defined by 



^(x) = Y J rf3x^(x,xXx') , (6.31) 



where 
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a(x) = A'p(x) + r(x) , (6.32) 

and Y being a constant, which we introduce for dimensional reasons. Then Eqs. (|6.26l) are equiv- 
alent to 



^^^#^^A3^(x,x') + |^A2^(x,x') - T^A^(x,x') = r-M(x-x') 

ai(2a2+a-i) ' / ' 2a2+a3 ' / 2a2+a-i ^\ ' J V / 

(6.33) 

AE{x)5i,-E{x),ij = 

where 5(x — x') is the 3-dimensional Dirac 5-function. The general solution of (16.271) for $(x) 
and ^(x), in terms of the Green function ^(x, x') and the function S(x), are 



A'p(x') 



-|^A,.S(x') + H^fe±^i2lA2,H(x') 

8a2+3a3 



(6.34) 



A'p(x') 



-g±MA,,S(x') + ':l^°^+"^M '.H(xO 



ai(2a2+a3) 



4a2+a3 ~;(x') 



2ai (202+13) ' 



Eqs. (16.271) represent a coupled set of fourth order differential equations. The total number of 
integration constant is eight. With the substitution (16.291) . it has been possible to decouple the set 
of equations, but now the differential order is changed. The total differential order is the same, 
indeed we have one equation of sixth order and another equation of second order, while previously 
we had two equations of fourth order. The number of integration constants is conserved. We can 
conclude that, with our approach, also introducing the new quantities $, ^ does not contradict the 
paradigm of a metric theory of HOG. The price is that now the r.h.s. of tt - component of field 
equation has been modified: there is an additional matter term r coming from the ij - component, 
(see the redefinition of the matter density (16.321) ). In Table [611, we show particular cases of Eqs. 
(16.271) for different choices of coupling constants of the theory with vanishing the determinant of 
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transformations (16.281) . 



6.3 Green functions for systems with spherical symmetry 

We are interested in the solutions of field Eqs. (16.271) at Newtonian order by using the method of 
Green functions and assuming a system with spherical symmetry: ^(x, x') = ^(|x — x'|). Let us 
introduce the radial coordinate r = |x — x'|; with this choice, the first equation of (16.331) for r 7^ 
becomes 

203(802 + as) A3^(x, x') + (6a2 + 03) A2^(x, x') - a? A,^(x, x') = , (6.35) 

where = r~'^dr{r~'^dr) is the radial component of the Laplacian in polar coordinates. The 
solution of (16.351) is: 



r 



K2 + -[ Kse 
ai 



2(3a2 + as) 



I ^1 , 

2(302+03) 



(6.36) 



where Ki are constants. We note that, if 02 = as = 0, the Green function of the Newtonian 
mechanics is found. It is the same of the Electromagnetism. The integration constants have to 
be fixed by imposing the boundary conditions at infinity and in the origin. In fact Eqs. (16.361) is 
a solution of (|6.35l) and not of the first equation in (16.331) . A physically acceptable solution has 
to satisfy the condition ^(x, x') — )• if |x — x'| — t- 00, then the constants Ki, K4, Kq in (|6.36l) 
have to vanish. To obtain the conditions on the constants K2, K^, we consider the Fourier 
transformation of ^(x, x'), that is 



^(x,x ) = ^(x - x') 



(27r)3/2 



(6.37) 



Eq. (16.331 ). in terms of Fourier transform, becomes 
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Vw-llvJlCCa Ui CUllslu-llls 


Corresponding field equations 


A 


02 = 

ag = 


U(x)+ ^ [dv,''^"'!,! =0 

^ ' ai J |x— x'l 


B 


ai = 

«3 = 


A 


(2^ - $) = 
A(2^ - $) 




A(2^-<l>) 


= 


C 


ai = 
02 = 


A^ 
A 


($ = ^ 
A($-3^) 




A($ - 3^) 


= 

,17 


D 


txg -^'J'2 


2a2A2$ - aiA^ = -f p, 

V2 ai$(x)-2a2A$(x) + G/rfV^ 5,,- 

- ai$(x) - 2a2A<l>(x) + G / rf^x'^ = 


E 


cti — U 
03 = — 4a2 


A2 
A 


$ = 
A^ 


Sij - 


A^ 


= 


F 


ai = 
as = — 2a2 


A A(^-$) 


S^j - 


A(^-<l>) 


= 

,«.7 


G 


ai = 

„ — St^z 

03 - — ^ 


A^ 
A 


(2^ 
A$ 


+ $) = 

6ij - 


_ 3r] 

4a 

A<l> 


= 

,j7 



Table 6.1: ExpUcit form of the field equations for different choices of the coupling constants for which the 
determinant of the transformations (16.281) vanishes. 
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(27r)3/2 



3ik-(x— x') 



203(802 + 03) g 602 + 03,4 



0i(202 + 03^ 

Oi 



-k^ + 



202 + O3 



202 + O3 



k^ + 



0. 



(6.38) 



The Green function can be expressed as follows: 



6?(x-x') = -Y 



-1 



gik-(x-x') 



ai (202+03) 202+03 202+03 



(6.39) 



Since we are investigating systems with spherical symmetry, it is better to introduce polar coordi- 
nates in the k-space. Eq. (16.391) becomes 



6^(x-x') 



Y 



-1 



d\k\ 



(27r)2/2 Jn 203(302+03) ^^4 _ 602+03 ^^2 _ 01 
01(202+03) 202+03 202+03 



01(202 + 03) 1 



dn e'^<''-''''^ 



(6.40) 



27r 03(302 + 03) |x-x'| 



d\k 


sin 


k| 


X — x' 




\k\ 


k2 - 


03 




lc-2 1 "1 

^ 2(802+03) 



The analytic expression of Q(x — x') depends on the nature of the poles of |k| and on the values 
of the arbitrary constants oi, 02, 03. If we define two new quantities Ai, A2 G TZ: 



03 



A 



Oi 



2(3o2 + 03) 



(6.41) 



we obtain a particular expression of (16.391) : 



^(x-x') 




TT Y 



-1 



18 X- x' 



Xj _ Xj e-Ai|x-x'| g-A2|x-x'| 



2 \2 
2 



A?A 



A? 



+ 



Ai 



(6.42) 



This Green function corresponds to the one in (16.361) but now we have also the conditions in the 
origin. Obviously, we have three possibilities to introduce Ai and A2. In Table 16.21 we provide the 
complete set of Green functions ^(x — x'), depending on the choices of the coefficients 02 and 03 
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(with a fixed sign of ai). 



Case 


Choices of constants 


Green function 


A 


3a2 + a3 > 


y ^) - V 18 |x-x'| A^Ai 1 A^ Xi 


B 


ag > 

3a2 + a3 < 




A|-A^ cos(Ai|x-x'|) 1 cos(A2|x-x'|) 








C 


as < 

3a2 + as < 




A2+a2 , e-^il---'l , cos(A2|x-x'|) 




^^(---') = v^^ 


AjA2 Aj A2 



Table 6.2: The complete set of Green functions for equations (I6.39I ). It is possible to have a further choice 
for the scale lengths which turns out to be dependent on the two knows length scales. In fact, if we perform 
the substitution Ai ^ A2, we obtain a fourth choice. In addition, for a correct Newtonian component, we 
assumed ai > 0. In fact when 02 = as = the field equations (16.261 ) give us the Newtonian theory if 
ai = 1. 



When one considers a point-like source, p oc 5(x), and by setting S(x) = 0, the potentials 
(16.341) are proportional to ^(x — x'). Without losing of generality, we have: 

^(x) ~ ^ + —— + —— , (6.43) 
|x| |x| |x| 

an analogous behavior is obtained for the potential ^^(x). We note that, in the vacuum case, we 
found Yukawa-like corrections to Newtonian mechanics but with two scale lengths related to the 
quadratic corrections in the Lagrangian (|6.18l) . See also the above expressions (|6.41l) . This behav- 
ior is strictly linked to the sixth order of (16.331) . which depends on the coupled form of the system 
of equations (16.261) . In fact if we consider the Fourier transform of the potentials $ and ^: 

$(k) e^^- , v,(x) = y v,(k) e^^- , (6.44) 



the solutions are 
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v^(x) = -f/ 



X_ r d^k [ai+(8a2+3a3)k^]p(k)e'''-^ 
' 2 J (27r)3/2 k2(ai-a3k2)[ai+2{3a2+a3)k2] 



X r rf^k [ai + (4a2+a3)k2]p(k)e'''-" 



(6.45) 



(27r)3/2 k2(ai-a3k2)[ai+2(3a2+a3)k2] 



where p(k) is the Fourier transform of the matter density. We can see that the solutions have the 
same poles as (16.391 ). Finally, if p(k) = - — --tt (the Fourier transform of a point-like source) the 
solutions (|6.45l) are similar to (|6.43l) . In fact, if we suppose that as 7^ and 3a2 + 03 7^ 0, the 
solutions (|6.45l) are 



<|.(x) 
M/(x) 

Then the metric (16.211) becomes 



, GM ( _ 4 -Ai|x| I lg-A2|x| 
ai|x| \ 3 3 



. GM / ^ _ 2g-Ai|x| _ lg-A2|x| 
ai|x| \ 3 3 



(6.46) 



ai X 



1 - ^e-^^l'^l + ^e- 



A2|x| 



1 + 



ai X 



2 

— ( 

3 



1 

— t 

3 



(6.47) 



It is interesting to note that, if 03 = (Ai — )• 00), we have the missing of a scale length (a pole 
is missed) with only a Yukawa-like term as for the Electrodynamics. The Green function, in this 



case, is: 



and the Lagrangian becomes: / = aiR + a2-R^- Since at the level of the Newtonian limit, as 
discussed, the powers of Ricci scalar higher then two do not contribute, we can conclude that 
(16.481) is the Green function for any /-gravity at Newtonian order, if / is an analytic function of 
the Ricci scalar. We found the same situation in the newtonian limit of / - theory in standard 
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coordinates (§ 16.11) . In fact if we consider the analogy 



oi = /i, 02 = /2 (6.49) 

and if sign [02] = — sign[ai] we found the same characteristic scale - length: A = A2. Finally the 
presence of the pole is achieved considering a particular choice of the constants in the theory, e.g. 
03 = — 2a2. In Table 1611 (Case D), we provide the field equations for this choice and the relative 
Green function is: 



W-a.v^)(k)oc^^^^^^^:^. (6.50) 

The spatial behavior of (|6.48l) - (16.501) is the same but the coefficients are different since the theories 
are different. In conclusion we need the Green function for the differential operator A^. The only 
possible physical choice for the squared Laplacian is: 



^(A2)(x - X ) oc -, (6.51) 

^ ' |x — x'l 

since the other choice is proportional to |x — x'| and cannot to be accepted. Considering the 
last possibility, we will end up with a force law increasing with distance [|3TI . Summary, we 
have shown the general approach to find out solutions of the field equations by using the Green 
functions. In particular, the vacuum solutions with point-like source have been used to find out 
directly the potentials, however it remains the most important issue to find out solutions when we 
consider systems with extended matter distribution. 



6.4 Solutions using the Green function 

Unlike the §|6T|we are going to find solutions with Green functions method. Before to investigate 
the general solution of Eqs. (16.261) we want discuss, in the first subparagraph, all cases shown in 
the Table 16. 1[ While in the next subsections we will analyze the solution in presence of matter 
using the Green functions shown in Table 16.21 
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6.4.1 Particular solutions 

In Table (63] we provide solutions, in terms of the Green function of the corresponding differential 
operator, of the field equations shown in Table 16.11 Case A corresponds to the Newtonian theory 
and the arbitrary constant ai can be absorbed in the definition of matter Lagrangian. The solutions 
are: 



^$(x) = A^(x) = -G / d'x'j^^^ . (6.52) 

For Case D, instead, we have the field equations of a sort of modified electrodynamic-like repre- 
sentation. The solution can be expressed as follows: 



D 



D 



G / d^x.' 



P(x') 



(6.53) 



The solutions make sense only if ai/a2 > 0, then we can introduce a new scale-length. A partic- 
ular expression of (|6.53l) . for a fixed matter density p(x), will be found in a more general context 
in the next section. Nevertheless these two cases are the only ones which exhibit the Newtonian 
limit (obviously the first one!), while for the remaining cases there are serious problems with di- 
vergences and incompatibilities. In fact. Case B presents an incompatibility between the solution 
obtained from the tt - component and the one from the ij - component. The incompatibility can 
be removed if we consider, as the Green function for the differential operator V*, the trivial solu- 
tion: G(a2)\b = const. With this choice, the arbitrary integration constant $o can be interpreted as 
—GM. However another problem remains: namely the divergence at the origin. The interpretation 
of the constant $o as a total mass and not as a generic integral / (i^x'p(x') does not avoid the 
singularity. We can conclude, then, the solution 



2b^(x) -b $(x) = -^^^ (6-54) 

holds only in vacuum. Before continuing our analysis of the various cases, the term J (i^x'^(^2) (x— 
x')/)(x') has to be discussed for the choice (16.511) . A generic field equation with (from Table 
[61]) is 
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A^$(x) (xAlf c^V ^^^'\ oc A,p(x) ^ 47rp(x); (6.55) 
J |x X I 

we conclude that the only consistent possibility is to set p(x) = 0. In the remaining cases, we can 
only consider vacuum solutions. 



Case 


Solutions 


Newtonian behavior 


A 


$(x)=vl/(x) = -f /rfV^ 


yes 


B 


2v^(x) - $(x) = ll 

2M/(x) - <|.(x) = -ME Jrf3^'^ (x-xOp(x') 


no 


C 


W = S + ^ / ^'^'^(v^)(^ - ^')P(X') 
^(^) = R + ^ / ^'x'^(v^)(x - x')p(x') 


no 


D 


$(x) = -AnG J rfV^(2a2V4-aiV2)(x-x')p(x') 
^(x) = -AttG J c/V^(2a2V4-aiV2)(x-x')p(x') 


yes 


E 


<^'(x) = -^/rfV^(v.)(x-x')p(x') 


no 


F 


^W = -^/^'x'^(v^)(x-x')p(x') 


no 


G 


= §1 

= -Ir - ^/rfV^(v^)(x-x')p(x') 


no 



Table 6.3: Here we provide the solutions of the field equations in Table \6A] The solutions are found by 
setting H = in the ij - component of the field equation (16.33b . The solutions are displayed in terms of the 
Green functions. <I>o is a generic integration constant. 
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6.4.2 The general solution by Green function ^ (x — x') 

In this section, we explicitly determine the gravitational potential in the inner and in the outer re- 
gion of a spherically symmetric matter distribution. The first consequence of the extended gravity 
theories which we are considering is the non-validity of the Gauss theorem. In fact, in the Newto- 
nian limit of GR, the equation for the gravitational potential, generated by a point-like source 



^^QNew.mechX^ - = -47r5(x - x') (6.56) 

is not satisfied by the new Green functions developed above. If we consider the flux of gravitational 
field gNew.mech. defined as 



',New.mech. — i — ^ 7iT~^ — —G MV xQ New.mech.{^ ~ y (6.57) 



|x — x'P 



we obtain, as standard, the Gauss theorem: 



gNew.mech. ■ fl (X M , (6.58) 

where S is a generic two-dimensional surface and n its surface normal vector. The flux of field 
gNew.mech. ou the surfacc S is proportional to the matter content M, inside the surface indepen- 
dently of the particular shape of surface (Gauss theorem, or Newton theorem for the gravitational 
field I1205II ). On the other hand, if we consider the flux defined by the new Green function, its value 
is not proportional to the enclosed mass but depends on the particular choice of the surface: 



C/S gNew.mech. ■ fl My, ■ (6.59) 



s 



Hence My, is a mass-function depending on the surface S. Then we have to find the solution 
inside/outside the matter distribution by evaluating the quantity 



(i3x'^?^(x - x')p(x ) , 



(6.60) 
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and by imposing the boundary condition on the separation surface. By considering solutions (|6.34l) 
with the Green function — x') from Table 16.21 and by assuming S(x) = 0, we have 



= {^^^ + ^l2A^) G/rfV^^(x-x')p(x') 

A^W = rA'(^gg^/rfV^^(x,x')p(x') 

= (/ii + /i3AJ G/rfW(x-x')p(x') 



where /ii := 



47ry 



2a2+a3 



12ttY AfA| 



47ry(8a2+3a3) _ 4ttY 4Af-A 



(6.61) 



47ry(4a2+a3) 
ai(2a2+a3) 



ai(2a2+a3) 

x2^x2 ■ We have to note that our working hypothesis, H(x) = 0, is not particular, since when 
we considered the Hilbert-Einstein Lagrangian in § l6.2l to give the Newtonian solution, we imposed 
an analogous condition. For the potential $(x), one obtains: 



^$(x) = (/ii + /i2Aj Gy"rfVp(x') Y^g^^ 



3-Ai|x-x'| 



X — X 



(6.62) 



i=0 



i=0 



where 



r 



A, Ai , 



G / rf'^xV(x' 



g-Ai|x-x'| 



(6.63) 



<l>^'^'(x) = (/ii + /i2Ax)r^'"»(x), (6.64) 

and gf are the coefficients. Here the numbers Aj assume the above values 0, Ai, A2. Supposing 
a matter density p(x) = p(|x|) and denoting the radius of the sphere with total mass M by ^, we 
have 
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G Trfix'iixfpdx'i) r, 

JO -'0 



/•TT -A,^|xP+|x'|2-2|x||x'|cosQ 

Jo a/|xP + |x'p — 2|x||x'| coso; 



(6.65) 



where cos a = cos 6^ cos 6*' + sin 6* sin 6^' cos(0 — (p') and a is the angle between two vectors x, 
x'. In the spherically symmetric case, we can choose = without losing generality. The (16.651 ) 
becomes 



r^'^'(x) 



2nG 

Jo 



d|x'||x'|p(|x'|) 



3M 



g-Ai||x|-|x'|| _ g-A,(|x| + |x'|) 



For a constant radial profile of density, p(|x|) = ^ we have 



(6.66) 



3GM 



1 - (1 + A,,Oe-^'«^i^^^^ 



f|||[A,ecosh(A,0 - sinh(A,0]^ 



in the inner and in the outer region respectively. The limit 



(6.67) 



lim r^'^^(x) 

A,-s>0 



3GM _ GM\ \2 
2( 2^3 l-^l 



GM 
|x| 



|X| > 



(6.68) 



gives us the internal and the external Newtonian behavior. The internal and the external potential 
for a given Aj is 



3GM 



A«|x| 



(6.69) 



^t^'i^) = f|f [A.ecosh(A.O - sinh(A.O](/^i + 
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The boundary condition on the surface |x| = ,^ is 



A$m(0 -A ^outiO 



3GM 



A 
i 1 



(6.70) 



1=0 



but the last relation is identically vanishing (see Table I6.2I ). The internal and external potential is 
given by: 



l$m(x) 



V2 7r3/^ GM 

ai e 



Af(2+3Aig^^)-8Ai _ ^ g^_,,^^^ ^ ^^^^ 



sinh(Ai|x|) 
A?|x| 



e--^ll^l 

1^ 



■?4f^§|[A2ecosh(A20 - sinh(A20] 



(6.71) 



The relations (16.711 ) give the solutions for the gravitational potential $ inside and outside the con- 
stant spherically symmetric matter distribution. A similar relation is found for ^ by substituting 
1^2 — ^ At3- We note that the corrections to the Newtonian terms are ruled by G^{x — x'). If we 
perform a Taylor expansion for A|x| <^ 1, we have: 



smh(A X ) , ,r, 

' ~ constant + xr + . . . . 

A|x| 



(6.72) 



For fixed value of the distance |x|, the external potential $out(x) depends on the value of the radius 
^, then we have that the Gauss theorem does not work. In this case, the potential depends on the 
total mass and on the distribution of matter in the space. In particular, if the matter distribution 
takes a bigger volume, the potential |$o«t(x) | increases and viceversa. We can write 



(6.73) 
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obviously the limit of ^ has to be interpreted up to the maximal value where the generic position 
|x| in the space is fixed. If we consider the limit ^ — )■ (the point-like source limit), we obtain 



lim 3 



A^cosh(AO -sinh(AO 



(6.74) 



For A^nutM, we have 



, , 2y/2 71^/^ GM 8^2 7r3/2 GMe-^iW 

hm A^outi^) = r-r H 5 n 

C-s-o ai |x| Sai |x| 

3ai |x| 

The last expression is compatible with the discussion in 



(6.75) 



6.4.3 Further solutions by the Green functions ^^(x — x') and (x — x' 



For the sake of completeness, let us derive solutions for the other Green functions. Starting from 
Case B in Table 16.21 we have: 



r 



3«^-^rsin(A.O-A.ecos(A.O]^^^ 



|x|<e 



|x| > ^ 



(6.76) 



in the inner and outer region. Also in this case, if we consider the limit of Aj — )■ 0, one obtains the 
Newtonian limit (|6.68l) . The internal and external potential for given Aj is 



m \ J 



3GM^.. \ 2 



A»|x| 



J + _ A,V2)[cos(A,0 + A,esin(A,0]^^^ 

cos(Aj|x|) 



(6.77) 



^(^1 - Ai /i2)[sin(AiO - cos (AiO 1^,7^ 
The boundary condition on the surface |x| = ^ is 
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3GM 

B^^niO -B^outiO = 7^/^2$^6^f = 

^ i=0 

(see Table l6.2l) . The internal and external potential are given by 



(6.78) 



V2 7r3/2 GM J Af(3A|g2_2)+8A| 

ai e ] AfA^ 1^1 A^ 



^ - |Xp - ^[cos(AiO + 



+Aiesin(AiO]^^ + 4[cos(A20 + sinlA^O]^^ 



^ + — Afi^[sm(Ai4j - Ai4 cos(Ai4jJ — ^ 

2Vlzr!^ GM r • / \ f\ \ (^Nl cos(A2|x|) 



N). 



(6.79) 



The above considerations hold also for the first of (16.791) . The correction term to the Newtonian 
potential in the external solution, second line of (|6.79|) . can be interpreted as the Fourier transform 
of the matter density. In fact, we have: 



and in the limit 



dh 



(27r)3/2 



p(x')e 



_,k.x' ^ 3M sin(|k|0-|k|ecos(|k|0 

(27r)2/3 |k|3^3 



(6.80) 



lim 



rP(x)e 



/\ „-ik-x' 



M 



>o J (27r)3/2^^ ^ (27r)2/3 ' 

we obtain again the external solution for point-like source as limit of external solution (|6.79l) : 



(6.81) 



lim B<l>o„t(|x|) 



2V2 7r3/2 GM ^ 4\/27r3/2 GMcos(Ai|x|) 



ai |x| 3ai 
x/27r=^/2 G'Mcos(A2|x|) 
3ai Ixl 



(6.82) 



Finally for Case C in Table l6.2l we have 
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V2jr^GM J Aj(3Aig2_2)_8A2 , ,2 g ^-Ai£/i , \ sinh(Ai|x|) 



+ 4[cos(A20 + A2esin(A20]^^ 



(6.83) 

■^4f^^[sin(A,0 - A2ecos(A,0]^^ 



The limit of point-like source is valid also in this last case: 



hm c^outi^) = r-r H 5 n 

C->o ai |x| Sai |x| 

2^2 7r3/2 GA'/cos(A2|x|) 

3ai Ixl 



(6.84) 



These results means that for suitable distance scales, the Gauss theorem is recovered and the theory 
agrees with the standard Newtonian limit of GR. 



6.5 Post — Newtonian scheme of / — gravity 

In this last paragraph of sixth Chapter we want to trace a methodological approach to "perturbed" 
Eqs. (11.251) when we consider / ~ /i_R+ /2-R^ + + . . . (see the Taylor expanse (16.71) ) and the 
metric tensor completed at post-Newtonian order (I2.35|) . Obviously the matter tensor is Eq. (I2.57|) . 
Eqs. (fL25]) - (fOgl) at 0(2) - order become 



rf-fi?(2)_2/^Ai?(2) = A'T,? 



(0) 



< m 



(2) 



(6.85) 



= _6/2Ai?(2) _ f^R{2) = ;f 2^(0) 
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at 0(3) - order 



(1) 

ti 



(6.86) 



and by remember the expressions for the Christoffel symbols (12.391 ) and In ~ | [gu^ — gmL] + 
. . ., finally, at 0(4) - order, 



-2/2 



-6/3 



|Vi?(2)|2 + i?(2)Ai?( 



2) 



XT 



(2) 



= -6/2Ai?(4) - - I8/3 



|Vi?(2)|2 + i?(2)Ai?(2) 



+6/2 



p(2) _ ^(2) p(2) _ lv7/^(2) _ ^(2) X „p(2)_ (2) p(2) 



(6.87) 



The solution for the Ricci scalar i?^^), in the last line of (|6.85l) is similar to (|5.24l) . In fact we have 



i?(2)(x) = X 



(27r)3/2 6/2k2 - A 



/i y (27r)3/2 k2 + A2 
/i J |x-x'| 



(6.88) 



where, again, A is one defined previously (16.131) or (16.411) . Considering a generalization of the 
metric (12.711) 



_ 1 + 2$ + 2 6 



(6.89) 
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and requiring the harmonic gauge (12.431) the solution for the metric tensor is found. Now, from 
the expression of time-time component of Ricci tensor at Newtonian level (|2.53l) we obtain the 
modified newtonian potential as solution of the first line of (|6.85l) 



We can check immediately that when / — )■ i? we find $(x) — t- —G J d^x' |^^J| . This outcome is a 
generalization of solution (16.151) . From the (16.861) we find the "vectorial" solution 



(6.9.) 

fi J \x- x'\ QttX^ J - 
and from second line of (16.851) the spatial tensorial solution 



(2) 

9ij 



fi 

2 

3A2 



I ,1 + . 
X — x' 4tt 



X — x' 



3A2 

|x| 



rf|x'||x|'2i?(2)(x') 



(6.92) 



We can affirm that it is possible to have solution non-Ricci-flat in vacuum: HOG mimics a matter 
source. 



From the fourth order of field equation, we note also the Ricci scalar (i?^'*)) propagates with 
the same A (the second line of (16.871 )) and the solutions at second order originates a supplementary 
matter source in r.h.s. of (|1.25l) : 



J_ A/?(^) - = l^T, + ^T, + ^ r(2) (6.93) 
A^ Ji A^ /i 

where the functions T^, T5 are known. The Ricci scalar at fourth order is 
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1 r A|x— x'l 

Also in this case we can have a non-vanishing curvature in absence of matter. At last the tt - 
component at fourth order can be reformulated as follows 



AO — — + contributions from previously order (6.95) 

/i 

then it is possible to find a generale solution for tt - component at fourth order of metric tensor. 

With this last paragraph we wanted to resume a methodological approach to Post-Newtonian 
limit of /-gravity, if / is an analytical function of Ricci scalar. The development of / is performed 
in = and the Ricci tensor components are expressed in the harmonic gauge. 



Chapter 7 



The post-Minkowskian approximation in / 
— gravity: Gravitational Waves in higher 
order gravity 

In this chapter, we develop the post-Minkowskian limit of HOG theories [F]. It is well known 
that when dealing with GR such an approach provides massless spin-two waves as propagating 
degree of freedom of the gravitational field while ETGs imply other additional propagating modes 
in the gravity spectra. We show that a general analytic HOG model, together with a standard 
massless graviton, is characterized by a massive scalar particle with a finite-distance interaction. 
We briefly discuss how such massive gravitational mode can have relevant consequences both on 
cosmological and small scales distances affecting the stochastic background of gravitational waves 
and representing a valid alternative to Dark Matter on galactic scales. Furthermore we develop an 
analytic definition of the energy-momentum tensor of the gravitational field in such a scheme. 
Such a tensor represents a basic quantity in order to calculate the gravitational time delay in Pulsar 
timing. 

7.1 The Post — Minkowskian approximation in spherically sym- 
metric solution 

In Chapter|6l we have found the spherically symmetric solution of /-gravity in the Newtonian limit 
with the metric (|1.89l) . Here we want to discuss a different limit of these theories, pursued when 
the small velocity assumption is relaxed and only the weak field approximation is retained. This 
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situation is related to the Minkowski limit of the underlying gravity theory as well as the discussion 
of the Chapter [6] was related to the Newtonian one. In order to develop such an analysis, we can 
reasonably resort to the metric (|1.89l) . considering the gravitational potentials guit, r) and grr{t, r) 
in the suitable form 



gu{t,r) = l + gil\t,r) 

(7.1) 

grr{t,r) = l + gil\t,r) 

with gu\gil'' ^ 1. Let us now perturb the field equations (11.251) . with respect to approach (13.311) 
- (13.331 ). considering, again, the Taylor expansion (16.71) for a generic / - theory. For the vacuum 
case (Tfj^y = 0), at the first order with respect to g^l^ e gil\ it is 



where 



/o = 



,(1) 



(7.2) 



(0) 



(7.3) 
(7.4) 



In this approximation, the Ricci scalar turns out to be zero while the derivatives, in the previ- 
ous relations, are calculated at i? = 0. Again we find /-Lagrangians without the cosmological 
contribution as in (|6.10l) . 



Let us now consider the limit for large r, i.e. we study the problem far from the source of the 
gravitational field. In such a case the (17.21 ) become 
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ytt yrr 



8/2 



q2„(1) #0*^' r)2„(l) 
O Srr I " 9tt Q O ffrr 

^ 9*2 



(7.5) 



any function of time t 



The (17.51) are two coupled wave equations in term of the two functions g^f* and gr^r ■ Therefore, we 
can ask for a wave-like solutions for the gravitational potentials 



(7.6) 



9 



(1) _ [■ dujdk~il) 



and substituting these into the (17.51 ). We find the condition 



9tt\^, k) 



grr{oo, k) UJ = ±k 



1 + 



32d 

4 fc2 



(7.7) 



gn\uj, k) uj = ± ■sj', 



where A is defined in the (16.131) . In particular, for /i = or /2 = one obtains solutions with a 
dispersion relation u = ±k. In other words, for fi ^ (i = 1, 2), that is in the case of non-linear 
/, the above dispersion relation suggests that massive modes are in order. In particular the mass 
of the graviton is mgrav = and, coherently, it is obtained for a modified real gravitational 
potential. As matter of fact, a gravitational potential deviating from the Newtonian regime induces 
a massive degree of freedom into the particle spectrum of the gravity sector with interesting per- 
spective for the detection and the production of gravitational waves ^20611 . It has to be remarked 
that the presence of massive gravitons in the wave spectrum of HOG is a well known result since 
the paper of [32] . Nevertheless it is our opinion that this issue has been always considered under a 
negative perspective and has been not sufficiently investigated. 

In the post-Minkowskian approximation, as expected, the gravitational field propagates by 
means of wave-like solutions. This result suggests that investigating the gravitational waves behav- 
ior of HOG can represent an interesting issue where a new phenomenology (massive gravitons) has 
to be seriously taken into account. Besides, such massive degrees of freedom could be a realistic 
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and testable candidate for cold dark matter, as discussed in [I207II . 



7.2 The post — Minkowskian limit of / — gravity toward grav- 
itational waves 

Let us formally develop, in this section, the post-Minkowskian limit of HOG models. Such inves- 
tigation completes the analysis of the weak field regime of /-gravity and it has to be considered in 
this present work thesis. The post-Minkowskian limit of whatever gravity theory arises when the 
regime of small field is considered without any prescription in term of the propagation velocity of 
the field. This case has to be clearly distinguished with respect to the Newtonian limit which, differ- 
ently, requires both the small velocity and the weak field approximations. Often in literature such 
a distinction is not clearly remarked and several cases of pathological analysis can be accounted. 
The post-Minkowskian limit of GR naturally furnishes massless waves as the propagating behavior 
of gravity in this regime. We can now develop an analogous study (see paragraph 12.41 ) consider- 
ing in place of the Hilbert-Einstein Lagrangian a general function of the Ricci scalar. Actually, 
in order to perform a post-Minkowskian development of field equations one has to implement the 
field equations (|1.25l) with a small perturbation on the Minkowski background rj^y (|2.74l) . It is 
reasonable to assume that the / - Lagrangian is an analytic expression in term of the Ricci scalar 
(16.71) {i.e. Taylor expandable around the Ricci scalar value R = Rq = 0). In such a case field 
equations (11.251 ). at the first order of approximation in term of the perturbation become : 



/i 



-2/2 



(7.8) 



From zero-order of (|1.25l) one gets again /(O) = while T^,^ is fixed at zero-order in (17.81) as in 
the paragraph 12. 4[ The explicit expressions of Ricci tensor and scalar are the same of (12.771) . The 
(|7.8I) can be rewritten in a more suitable form by introducing the constant A (|6.13l) : 



+ - V,un,){h^r^-^ - D,h) = (7.9) 
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and by choosing the harmonic gauge (12.431 ): h^^ = — ^r]^„ with the condition h^^" ^ = 0, one 
obtains that field equations and the trace equation respectively read 



(7.10) 

In order to deduce the analytic solutions of (17.101) . we can now adopt a dual space (momentum 
space) description, this approach can simplify the equations system and, above all, allows to di- 
rectly observe what are the physical properties of our problem. In such a scheme we have : 



(7.11) 

where 



(7.12) 

are the Fourier transforms of the perturbation h^p{x) and of the matter tensor t'i^J . We have 
defined, as usual, kx = — k ■ x and k"^ = u"^ — k^. On the other side h{k) and T^^\k) are the 
traces of h^v{k) and TjiJ{k). In the momentum space one can easily recognize the solutions of 
(17.111) : the expression for h^,y{k) turns out to be 



M^-) ^ f ^ f ^^^T<«.W, (7.13, 

which fulfils the condition h^^"^ ^ = {h^^^ik) k^ = 0). The true perturbation variable h^y{k) 
can ne obtained inverting the relation with the tilded variables, in particular inserting the matter 
functions Sf}{k) = T^°^(fc) - lr]^^T^^^{k) and5W(A;) = r]"^ S^^) (k) , one ohtaim: 
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H,Ak) - ^ (7.14) 

which represents a wavelike solution, in the momentum space, with a massless and a massive 
contribute since the pole in the denominator of the second term, whose mass is directly related 
with the pole itself. The explicit wavelike solution can be obtained returning the the configuration 
space inverting the Fourier transform of h^^. 

Let us remark that field equations (11.251 ), for a generical / - Lagrangian, can be rewritten 
isolating the Einstein tensor in the l.h.s. as usual for Curvature Quintessence [l56l[57l[58l . In such 
a case higher than second order differential contributes, in term of the metric tensor, are considered 
in the r.h.s. as a source component of the space-time dynamics as well as the energy momentum 
tensor of ordinary matter does : 



G,. = T(-) + r(-™), (7.15) 

where 



(7.16) 



rp{curv) _ 1 f-f'R I f'^y-Qi^yOf 



Actually if we consider the perturbed metric (|2.74l) and develop the Einstein tensor up to the first 
order of perturbation we have 



1 11 

G^u ~ G^l^ = - 2^'?V - 2^,f^'^ - ^ih^r"^ - Dr,h)r]^^ (7.17) 

while the curvature tensor will give the other contributes 



Tjfr^ - ^(^M^n. - dl,)ih.r''^ - Or,h) (7.18) 
This expression easily allow to recognize that, in the dual space of momentum, such a quantity 
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will be responsible of the pole-like term which implies the introduction of a massive degree of 
freedom into the particle spectrum of gravity. In fact, inserting these two expressions into the the 
field equations (|7.15l) and considering the (12.771) we obtain the solution : 



2X 

1\ 



(7.19) 



where ^^^{x) is related with the Curvature tensor and is defined as 



^-(^) = -5 / ^%^|^S<»)(*) . (7.20) 

The general solution for the metric perturbation h^y{x), when equation are given as in (17.151) . can 
be rewritten as 



V(^) = ^ / 7^^^ e^'^ I ^'^P^^^^^S'^'^k) e^'^ , (7.21) 



/i J (27r)2 P 3/i7 (27r)2 k^{k^ - X^) 

which displays in the second term a pole whose properties can be easily evaluated in vacuum. In 
fact, in such a case (i.e. T^,^ = 0), the (17.101) becomes 



(7.22) 

showing that allowed solutions are of two types along the relations : 



±|k| 



(7.23) 



(2^) 



Akx 



with h{k) = 



and 
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< (7.24) 

h{k) e'^'' with h{k) ^ 

It is evident, that the first solution represents a massless graviton according with standard pre- 
scriptions of GR while the second one gives a massive degree of freedom with m? = X'^. In this 
sense we can furtherly rewrite the (|7.10l) introducing (p = Oh so that the general system can be 
rearranged in the following way 



(2^)2 



'im? 



(7.25) 



/: 



-m 



which suggests that the higher order contributes act in the post-Minkowskian limit as a massive 
scalar field whose mass depends on the degree of deviation (/', /"), calculated on the background 
unperturbed metric, of the initial / - Lagrangian with respect to the standard Hilbert-Einstein 
expression. 

It is important to remark that the peculiarity of a massive contribute in the wave spectrum of 
HOG is strictly related with the peculiar behavior of the trace equation with respect to GR. In fact 
in the case of HOG theories the trace equation establishes a constraint for the Ricci scalar under the 
form of a dynamical equation. This relation allows a more complex evolution of the system since 
the Ricci scalar is not univocally fixed by the trace equation as in GR. In fact, while in the frame- 
work of GR vacuum solutions imply i? = {i.e. it holds the Poisson equation in the Newtonian 
Limit), in the HOG models R can assume a generical dynamical evolution according with (|1.26l) . 
which assumes the zero value only under certain hypotheses on the nonlinear Lagrangian. This 
behaviour, as widely displayed in chapter[6l implies, as natural consequence, a modify of the Pois- 
son equation in the Newtonian limit to a form which allows a modified gravitational potential in 
such a regime. This characteristic is directly related with the massive degree of freedom obtained 
in the post-Minkowskian limit of these theories. In other words, such peculiarity is a different 
representation, at a different scale (or energy range), of the same effect. One can easily notice that 
the characteristic length of the modified gravitational potential enters in the wave solution exactly 
as the mass parameter = — ^ of the additive component in the gravitational wave spectrum. 
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7.3 Strong gravitational waves in a general / — gravity 

We are interesting to study the field equation for a small perturbation h^i, on the background metric 
g'^^J {0{hY ^ 1). Where g^l^} is a solution of GR with Ricci scalar R = R'^^^ = (this solution is 
in the vacuum). Then the relativistic invariant is described as 



ds^ = g^rdx'^dx^ = {g^^^ + Kr)dx''dx^ (7.26) 

Obviously the lowering and rising of the index have been made with the metric background gj^J . 
The field equation (ll.25|) at zero order, if we consider the development shown in the paragraph 13 .41 



/'^°^4°^-^^7i°V^°^ = '^T(°) (7.27) 

from the which the trace equation states /*^°^ = —^T. But in the vacuum the trace is vanishing 
and we have to impose the condition f^^^ = (we neglect the cosmological contribute). Then the 
field equation (17.271) becomes 



4°) = (7.28) 

and the metric gj^J is solution for the field equation (11.251) at zero order in the vacuum (Schwarzschild 
solution). At the first order we have 

- Ig^^jR^''] - r^''[Ril - e^n.(o)i?(^)} = XTj^'J (7.29) 

where the derivatives covariant have been calculated with respect to metric gjS- Now the expres- 
sions for Ricci tensor and scalar are 



p(l) Ucr 1| — I /, 1 !, 



(7.30) 
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and the general field equation perturbed (17.291) is 



1^ (v,V. - ^7i°jn,,o,) (^Kr''''^ - ng,,h}j = (7.31) 



2 

+ 
3A2 



where A = — . If f R one has /" — 0, A — )• oo and the (17.311) becomes 



3/"(0) 



/.e. the first order for Einstein equation in GR with the same conditions for the metric g^J . The 



trace of (17.311) is 



□^(0) (h^r'"^ - Dgmh^ + P (^Kr"'^ - n^fo/i^ = r(i) . (7.33) 

By using the harmonic gauge condition (12.431 ) from the equations (17.311 ) and (17.331 ) one has 



(7.34) 



^ (□,(o)+A2)n^(o)/i = --p^A2T« 

7.4 Energy-momentum tensor of / — gravity 

In order to detect gravitational waves the construction of a number of sensitive detectors for grav- 
itational waves (GWs) is underway today. At the moment there are several laser interferometers 
already built like the VIRGO detector (Italy), the GEO 600 detector (Germany), the two LIGO 
detectors (United States), the TAMA 300 detector (Japan) and many bar detectors are currently 
in operation too. Since very soon there will be a huge amount of experimental data the results 
of these detectors will have a fundamental impact on astrophysics and gravitation physics. GW 
detectors will be of fundamental importance in order to probe GR and, above all, to check every 
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alternative theory of gravitation 11208 [ 12091 12101 . A possible target of these experiments is the 
so called stochastic background of gravitational waves [|220l 12211 IITTI l2T2l l2T3l l2T4ll which can 
be related with the inflationary scenario settled in the early universe evolution. Actually there is 
another very challenging test dealing with gravitational waves phenomenology: the gravitational 
time delay in Pulsar timing. This experiment is one of the most important evidence of GR validity, 
since allows to verify the correction to the orbital period of pulsars as predicted by Einstein gravity 
theory [121 511 . Therefore, this experiment represents an unescapable test in order to check a what- 
ever viable gravity theory. An analytic calculation of this problem has been performed in the case 
of the Brans - Dicke theory with positive results since there are not significant constraint on the 
Brans - Dicke parameter uj Il216[|217l . Actually, in order to calculate what is the physical effect of 
HOG model on a pulsar system one has calculate the energy-momentum tensor of the gravitational 
field. This quantity will characterize the energy loss due to the gravitational irradiation. Although 
the procedure to calculate the energy-momentum tensor of the gravitational field in GR is often 
debated, one can extend the formalism developed for a generical field theory and obtain this quan- 
tity varying functionally on the gravity Lagrangian in term of the Lagrange operator obtaining a so 
called pseudo-tensor whose properties does not completely fulfils dippheomorphisms invarianceQ. 
Such calculation need to be extended when dealing with an HOG model since field equations are 
of order higher than two. 

In standard field theory, given a generical Lagrangian £ = C{g^v, 9fj.u,p, 9fj.u,pa) depending even 
on accelerations, field equations are obtained considering a variational principle which considers 
all the explicit functions. Thus, in the case of a HOG Lagrangian which depends on the metric and 
its derivatives up to the second order one has 



5 j (i^x^f^f = 5 j (fxC{g^y,g^y^p,g^y, 



pa) 



^AT^ + d%- 5g^ 



(fx^^RP^Sp^ = (7.35) 



where ~ means we neglected a pure divergence. Then we can set 



<fxd\ 



dg, 



dC ^ dC 



pcr,X 



■dg 



pa,XS, 



^9 pa + — C'5'p<7,C 



dg, 



pa,Xi 



0, 



(7.36) 



'This quantity is tipically referred as the Landau-Lifshitz energy-momentum tensor, nevertheless other kinds of 
energy-momentum tensor of the gravitational field can be defined 
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As matter of facts, one can write generalized Euler-Lagrange equations for this framework : 



-9 



- dx- + d^^ 



dC 



9g, 



per 



dgpa, 



A? 







(7.37) 



which coincide with the field equations (ll.25|) in the vacuum. Actually, even in the case of such 
general model it is possible to define an energy momentum tensor of the field, in particular t'^ ^ 
turns out to be defined as follows : 



-9 



dC ^ dC 



pcr,X 



pa,X^ 



dC 



dg, 



pa,\^ 



(7.38) 



The (17.381) quantity together the energy-momentum tensor of matter T^u satisfies a conservation 
according with standard requirements. In fact since in presence of matter H^^ = X T^^, one has 



(7.39) 



and as a consequence 



(7.40) 



which demonstrates the conservation law. We can now write down the expression of the energy- 
momentum tensor ^ of the gravitational field in term of the / - gravity action and the respective 
derivatives, in such a way to have a completely general expression : 



-f'R 



dR _ 

dgpa,\ 

OR 



-A 



-g- 



dR 



-g ' V dgp^^xii 



dR 

gpa,a ~l~ 7^ gp(J,^o 



dg, 



pa,X^ 



dg, 



pa,\^ 



'gpcT,a ^af ■ 



(7.41) 



Let us notice that while in GR t'^ ^ is non-covariant quantity, the relative generalization in 
HOG models turns out to satisfy the covariance prescription behaving as an ordinary tensor. One 
can easily verify that such an expression reduces to the usual definition of the Landau-Lifshitz 
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energy-momentum tensor of GR when f = R. 



where the GR Lagrangian has been considered in its effective form, i.e. the symmetric part of the 
Ricci tensor, which effectively characterizes the variation principle leading to the motion equations 

It is important to stress that GR definition of the energy-momentum tensor and HOG definition 
are quite different. These discrepancies are due to the presence, in the second case, of higher 
than second order differential term in the gravity action, which cannot be discarded by means of a 
boundary integration as it is done in GR. We have remarked above that the effective Lagrangian of 
GR turns out to be the symmetric part of the Ricci scalar since the second order terms present in 
the definition of R can be discarded by means of integration by part. 

A generic analytic / - Lagrangian is characterized from the dynamical point of view only by 
the first two terms of its Taylor expansion once the perturbation is implemented at the linear level, 
i.e. f ~ /'(O) R + J^{R), where the function satisfies the condition: hm^j^o-^ R^- As a 
consequence one can rewrite the explicit expression of (|7.41l) as : 



dR 



-.3, 



pcr,X 



dR 



'd v-9- 



dR 

9pcr.a ~l~ ~^ 9pcr,^a 



OgpaM J J pa, 



pcr,X^ 



—J^"R^^— (7po-,a — S^J^ 



' 99 pa, 



(7.43) 



Let us recall the general expression of the Ricci scalar (11.51) splitting its linear (i?*) and quadratic 
{R) dependence once a perturbed metric is considered 



R = R* + R (7.44) 

(notice that £gr = —V—gR)- Actually, in the case of GR the Landau - Lifshitz tensor shows 
as a first non vanishing term a h'^ contribute. A similar result can be obtained in the case of HOG 
models. In fact considering the expression (17.431 ) one obtains that at the lower expansion order ^ 
reads : 
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OR* \ OR* 

C"o I 9 per, a ~l~ 9pcr,^a 



ogpaM ^ 



dR* 

-d^lR*- ]gp^^, 



9gpa,\s, 

Now, since for a perturbed metric (I2.74h R* ~ R^^\ where R^^^ is defined as in (12.771) . one has 



(7.45) 



dR* ^ dR^^^I _ r^pXr^uS, _ ^X^^pcr 



(7.46) 



and the first significant term in (17.451) is of second order in the perturbation. We can now write 
down the expression of the energy - momentum tensor explicitly in term of the perturbation h: 



in term of the tilded perturbation hnu the new contribution reads : 



(7.47) 



(7.48) 



As matter of facts the energy-momentum tensor of the gravitational field, which expresses the 
energy transport of this field during its propagation, can have a natural generalization in the case 
of HOG models. We have adopted in our case the Landau-Lifshitz definition, however some other 
approaches are in order as outlined in [i219il . The general definition of t'^ ^ obtained above consists 
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of a sum considering the GR contribute plus a term which takes into account corrections induced 
by the higher differential order of / - theories : 

t\ = /'(0)t\|,, + r(0)t\|^, (7.49) 

and again when f = Rwe obtains t"^ a = ^^Igr already discussed. 

Quantities obtained along this section represent the basic elements in order to develop an an- 
alytic calculation of the gravitational time delay in the pulsar timing in the framework of HOG 
models. Nevertheless such analysis is beyond the purposes of the current study and will be argu- 
ment of a forthcoming investigation. 



Chapter 8 



Discussions and conclusions 



ETGs are good candidates to solve several shortcomings of modern astrophysics and cosmology 
since they seem, in a natural way, to address the problem of cosmological dynamics without intro- 
ducing unknown forms of dark matter and dark energy (see e.g. lfT0ll56| '). Nevertheless, a "final" 
alternative theory solving all the issues has not been found out up to now and the debate on modi- 
fying gravitational sector or adding new (dark) ingredients is still open. Beside this general remark 
related to the paradigm (extending gravity and/or adding new components), there is the method- 
ological issue to "recover" the standard and well-tested results of GR in the framework of these 
enlarged schemes. The recover of a self-consistent Newtonian limit (or a weak field limit)is the test 
bed of any theory which pretends to enlarge or correct the GR. In fact GR has been consistently 
tested in physical situations implying, essentially, spherical symmetry and weak field limit [I129II . 
One of the fundamental and obvious issue that any theory of gravity should satisfy is the fact that, 
in absence of gravitational field or very far from a given distribution of sources, the spacetime has 
to be asymptotically flat (Minkowski). Any alternative or modified gravitational theory (beside the 
diffeomorphism invariance and the general covariance) should address these physical requirements 
to be consistently compared with GR. This is a crucial point which several times is not considered 
when people is constructing the weak field limit of alternative theories of gravity. 

In our opinion, such a task has to be pursued in the natural frame of the theory otherwise 
the results could be misleading. Specifically, we have to develop the limit in the Jordan frame 
without conformal transformations to the Einstein frame since such transformations could alter the 
interpretation of the results. 

In this PhD thesis, we have considered the Taylor expansion of a generic /-theory, obtaining 
general solutions in term of the metric coefficients up to the second order of approximation when 
matter is neglected. In particular, the solution relative to the gu metric component gives the grav- 
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itational potential which is corrected with respect to the Newtonian one of / = R. The general 
gravitational potential is given by a Yukawa-like terms, combined with the Newtonian potential, 
which is effectively achieved at small distances. Besides also starting from the standard corrections 
to the Hilbert - Einstein Action (the well-known quadratic ones), but now the matter is present, we 
have faced, in same systematic way, the problem to find out solutions. The solutions are found 
using the Green function method and we have derived several solutions where the Newton poten- 
tial is corrected by combinations Yukawa-like terms. We have classified the results considering i) 
the parameters in the Lagrangian, ii) the field equations and Hi) the resulting potential. In rela- 
tion to the sign of the characteristic coefficients entering the gu component, one can obtain real or 
complex solutions. In both cases, the resulting gravitational potential has physical meanings. A 
discussion on the non- validity of the Gauss theorem has been given. Furthermore, for spherically 
symmetric distributions of matter, we discussed the inner and the outer solutions. Furthermore, it 
has been shown that the Birkhojf theorem is not a general result for / - gravity. This is a funda- 
mental difference between GR and HOG. While in GR a spherically symmetric solution is static, 
here time-dependent evolution can be achieved depending on the order of perturbations [C, E, G]. 

From other hand it is possible also to calculate Newtonian limit of such theories with a redefi- 
nition of the degrees of freedom by some scalar field leading to the so called O'Hanlon Lagrangian 
[I203II . In fact, considering this latter approach, we get a scalar- tensor theory with vanishing kinetic 
term and a potential term linked to /-theory. Also in this case we found a Yukawa-like correction 
to classic Newtonian potential. Nevertheless when we turn off the modification of Hilbert-Einstein 
Lagrangian we do not obtain the right Newtonian potential. In fact only in this limit / — )■ i? it has 
sense speaking about the Eddington parameter 7 and its value is 1/2 and not 7 ~ 1 as observed. 
The origin of inconsistency is in the not-well defined field equation when / — )• _R. In fact this 
problem in present also in Brans-Dicke theory and only by requiring ubd — ^ C)0 we obtain the GR 
[H]. 

We have discussed the differences between the post-Newtonian and the post - Minkowskian 
limit in / - gravity. The main result of such an investigation is the presence of massive degrees 
of freedom in the spectrum of gravitational waves which are strictly related to the modifications 
occurring into the gravitational potential. This occurrence could constitute an interesting oppor- 
tunity for the detection and investigation of gravitational waves. To do this it needs to generalize 
the energy -momentum tensor for a generic /-gravity. In the last chapter we tried to find a new 
expression for a HOG [F] . 

Starting from Tensor-multi- scalar theory of gravity ll23l we can show how a polynomial La- 
grangian in the Ricci scalar R, compatible with the PPN-limit, can be recovered in the framework 
of HOG. The approach is based on the formulation of the PPN-limit of such gravity models de- 
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veloped in analogy with scalar-tensor gravity HI 3711 . In particular, considering the local relations 
defining the PPN fourth order parameters as differential expressions, one obtains a third-order 
polynomial in the Ricci scalar which is parameterized by the PPN-quantity 7 and compatible with 
the limit (3 = 1. The order of deviation from the linearity in R is induced by the deviations of 
7 from the GR expectation value 7 = 1. Actually, the PPN parameter 7 may represent the key 
parameter to discriminate among relativistic theories of gravity [A]. 

Besides we investigated also the viability to find spherical solutions is in /-theories with an 
perturbation methodic analysis with respect to standard results of GR when we consider the limits 
r — )■ 00 and / — j- i?. Essentially, spherical solutions can be classified, with respect to the Ricci 
curvature scalar R, as R = 0, R = Rq ^ 0, and R = R{r), where Rq is a constant and R{r) 
is a function of the radial coordinate r. In order to achieve exact spherical solutions, a crucial role 
is played by the relations existing between the metric potentials and between them and the Ricci 
curvature scalar. In particular, the relations between the metric potentials and the Ricci scalar 
can be used as a constraint: this gives a Bernoulli equation. Solving it, in principle, spherically 
symmetric solutions can be obtained for any analytic / function, both for constant curvature scalar 
and for curvature scalar depending on r. Such spherically symmetric solutions can be used as 
background to test how generic / theories of gravity deviate from GR. Particularly interesting are 
theories that imply f R in the weak field limit. In such cases, the experimental comparison 
is straightforward and also experimental results, evading GR constraints, can be framed in a self- 
consistent picture [I128II . Finally, we have constructed a perturbation approach in which we search 
for spherical solutions at the Ot/i-order and then we search for solutions at the first order. The 
scheme is iterative and could be, in principle, extended to any order in perturbations. The crucial 
request is to take into account / - theories which are Taylor expandable about some value R = Rq 
of the curvature scalar. A important remark is in order at this point. Considering interior and 
exterior solutions, the junction conditions are related to the integration constants of the problem and 
strictly depend on the source (e.g. the form of T^u)- We have not considered this aspect here since 
we have, essentially, searched for vacuum solutions. However, such a problem has to be carefully 
faced in order to deal with physically consistent solutions. For example, the Schwarzschild solution 
R = 0, which is one of the exterior solutions which we have considered, always satisfies the 
junction conditions with physically interesting interior metric. This is not the case for several 
spherically symmetric solutions which could give rise to unphysical junction conditions and not be 
in agreement with Newtons law of gravitation, also asymptotically. In these cases, such solutions 
have to be discarded [D]. 

We have discussed a general method to find out exact solutions in ETGs when a spherically 
symmetric background is taken into account. In particular, we have searched for exact spherically 
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symmetric solutions in /-gravity by asking for the existence of Noether symmetries. We have 
developed a general formalism and given some examples of exact solutions. The procedure con- 
sists in: i) considering the point-like / Lagrangian where spherical symmetry has been imposed; 
ii) deriving the Euler-Lagrange equations; Hi) searching for a Noether vector field; iv) reducing 
dynamics and then integrating the equations of motion using conserved quantities. Viceversa, the 
approach allows also to select families of / models where a particular symmetry (in this case the 
spherical one) is present. As examples, we discussed power law models and models with constant 
Ricci curvature scalar. However, the above method can be further generalized. If a symmetry ex- 
ists, the Noether Approach allows transformations of variables where the cyclic ones are evident. 
This fact allows to reduce dynamics and then to get more easily exact solutions. These consid- 
erations show that the Noether Symmetries Approach can be applied to large classes of gravity 
theories. Up to now the Noether symmetries Approach has been worked out in the case of FRW - 
metric. In this PhD work, we have concentrated our attention to the development of the general 
formalism in the case of spherically symmetric spacetimes. Therefore the fact that, even in the case 
of a spherical symmetry, it is possible to achieve exact solutions seems to suggest that this tech- 
nique can represent a paradigmatic approach to work out exact solutions in any theory of gravity. 
A more comprehensive analysis in this sense will be the argument of forthcoming studies [B]. 
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